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Abstract 

This review outlines our present experimental knowledge and theoretical understanding of deep- 
inelastic scattering on nuclear targets. The emphasis is primarily on nuclear coherence phe- 
nomena, such as shadowing, where the key physics issue is the exploration of hadronic and 
quark-gluon fluctuations of a high-energy virtual photon and their passage through the nuclear 
medium. New developments in polarized deep-inelastic scattering on nuclei are also discussed, 
and more conventional binding and Fermi motion effects are summarized. The report closes with 
a brief outlook on vector meson electroproduction, nuclear shadowing at very large Q 2 and the 
physics of high parton densities in QCD. 
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1 Introduction 



This review is written with the intention to summarize and discuss nuclear phenomena 
observed in the deep-inelastic scattering (DIS) of leptons (mostly muons and electrons) 
on nuclear targets. Experimental developments in the last decade have brought several 
such effects into focus (the EMC effect; shadowing, etc.). This first came as a surprise. 
At the high energy and momentum transfers characteristic of DIS one did not expect to 
see sizable nuclear effects which usually occur on length scales of order 1 fm or larger, 
governed by the inverse Fermi momentum of nucleons in nuclei. 

Today such nuclear effects are well established by a large amount of high-quality experi- 
mental data. Also, their theoretical understanding has progressed in recent years, so that 
an updated review of these developments appears justified. Our presentation, however, 
does not aim for completeness in all details. We wish to emphasize those effects in which 
two or more nucleons act coherently to produce significant deviations from the incoherent 
sum of individual nucleon structure functions. The most prominent effect of this kind is 
shadowing. Its close relationship with diffraction in high-energy hadronic processes is now 
quite well understood, which points to the significance of optical analogues when dealing 
with the interaction of high-energy, virtual photons in a nuclear medium. 

Other, less pronounced nuclear effects such as binding and Fermi motion will also be 
discussed. Some overlap with previous reviews |],|]||f||| is not unwanted for reasons 
of continuity. At the same time, this report incorporates plenty of more recent material, 
including polarization observables in DIS on nuclei, a field in which experimental activities 
progress rapidly and forcefully. 

Before turning to our main subject it is necessary and useful to summarize, in the following 
Section 2, our knowledge on free nucleon structure functions. Deep-inelastic scattering 
probes the substructure of the nucleon with very high resolution down to length and time 
scales of order 10~ 2 fm. The QCD analysis of the structure functions gives detailed insights 
into the composition of nucleons in terms of quarks and gluons, and their momentum and 
spin distributions. 

A fundamental question from a nuclear physics perspective is then the following: how 
do the quark and gluon distributions of the nucleon change in a nuclear many-body 
environment? What are the mechanisms responsible for such changes? These issues will 
be addressed starting from Section 3 in which the basic observations and phenomenology of 
DIS from nuclear targets will be described. A particularly instructive way of illustrating 
the physics content of nuclear structure functions is provided by a space-time (rather 
than momentum space) analysis to which we refer in a separate Section 4. Shadowing is 
discussed in detail in Section 5. Binding effects, Fermi motion and pionic contributions are 
dealt with in Section 6. Section 7 turns to a discussion of more recent work on polarized 
DIS from nuclei. A status summary and further perspectives follow in Section 8. 
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We close this introduction with a remark on references to the literature. As usual, aiming 
for completeness is an impossible task. What we hope to be a reasonable compromise is 
a combination of references to previous reviews in which earlier references can be found, 
together with selected original references to data and theory whenever they are of direct 
relevance in the text. 



2 Structure functions of free nucleons 

2. 1 Deep-inelastic scattering: kinematics and structure functions 

Consider the scattering of an electron or muon with four-momentum k^ = (E, k) and 
invariant mass m from a nucleon carrying the four- momentum P M = (E p , P) and mass M. 
Inclusive measurements observe only the scattered lepton with momentum k'^ = (E', k') 
as indicated in Fig.[2.1|. 



Neglecting weak interactions which are relevant at very high energies only, the differential 



k 




P 



Fig. 2.1. Inclusive deep-inelastic lepton- nucleon scattering. 



cross section is given by:Q 



d 2 a a 2 E' 



(2.1) 



dVLdE 1 Q 4 E 



For introductions to deep-inelastic lepton scattering see e.g. Refs.^^Q. 
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to leading order in the electromagnetic coupling constant a = e 2 /4tt ~ 1/137. Here 

qP = jfe" - k"* = (u, q) (2.2) 

is the four-momentum of the exchanged virtual photon, and Q 2 = —q 2 . The lepton-photon 
interaction is described by the lepton tensor L^ v . Let the spin projections of the initial 
and final lepton be s and s'. After summing over s' the lepton tensor can be split into 
pieces which are symmetric and antisymmetric with respect to the Lorentz indices ji and 
v. 

L^(k, s; k') = L%{k; k>) + i L* u (k, s; k'), (2.3) 
with: 

L^(fc; k') = 2 (kX + kX) + 9,u q\ (2.4) 
L^(k,s;k') = 2me flua ps a q f3 , (2.5) 

where the lepton spin vector is defined by 2m s a = u^ a ^u. For unpolarized lepton scat- 
tering the average over the initial lepton polarization is carried out. In this case only the 
symmetric term, L^ u , remains. 

The complete information about the target response is in the hadronic tensor W^ u . 
We denote the nucleon spin by S. Gauge invariance and symmetry properties allow a 
parametrization of the hadronic tensor, 

W^(q; P, S) = W* v (q; P) + i W*(q; P, S), (2.6) 
in terms of four structure functions. The symmetric part is 

W^P)={^-9^W 1 (P-q,q 2 ) 
and the antisymmetric part can be written 

G 2 (P-q,q 2 



W*(q;P,S) = e flua(S q c 



S p MGi{P ■ q,q 2 ) + (P • qS - S ■ q P?) 



M 



(2.8) 
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Here the nucleoli spin vector S^, with 2MS^ = U (P, S) 7^75 U (P, S), has been introduced. 
In the conventions used in this paper nucleon Dirac spinors U are normalized according 
to U (P) , ~f fl U(P) = 2P M . It is evident that W\p can be measured in unpolarized scattering 
processes, whereas the complete investigation of G^ 2 requires both beam and target to 
be polarized. 

It is convenient to introduce dimensionless structure functions 



F 1 (x,Q 2 )=MW 1 (P-q,q 2 



F 2 (x,Q 2 



P-q 



W 2 (P-q,q 2 ), 



which depend on the Bjorken scaling variable, 

Q 2 



X 



2P-q 



(2.9) 
(2.10) 

(2.11) 



In terms of the charged lepton scattering cross section (|2.1|) for an unpolarized lepton 
and nucleon is: 



cPa 4ira 2 



dxdQ 2 Q 4 



Mxy\ F 2 



(2.12) 



with 



P-q 
P-k' 



(2.13) 



Let us recall the behavior of the structure functions in the Bjorken limit, i.e. at large 
momentum and energy transfers, 



Q' 2 



-q 



00, P ■ q — > 00, 



(2.14) 



but fixed ratio Q 2 /P ■ q. Here the unpolarized structure functions 



F 2 (x,Q 2 



F 2 (x) 



(2.15) 
(2.16) 



are observed to depend in good approximation only on the dimensionless Bjorken scaling 
variable x. Variations of the structure functions with Q 2 at fixed x turn out to be small. 



A similar scaling behavior is expected for the spin-dependent structure functions: 



8 



g l (x J Q 2 )=MP-qG 1 (P-q,q 2 ), (2.17) 

In 

M 



92(x,Q 2 )= i -^ T ^G 2 (P- q ,q 2 ), (2.18) 



which likewise reduce to functions of x only when the limit Q — > oo is taken. 



2.2 Parton model 



The approximate Q 2 -independence of nucleon structure functions at large Q 2 has led to 
the conclusion that the virtual photon sees point-like constituents in the nucleon. This 
is the basis of the naive parton model which gives a simple interpretation of nucleon 
structure functions. In this picture the nucleon is composed of free pointlike constituents, 
the partons, identified with quarks and gluons. Introducing distributions qj(x) and qj{x) 
of quarks and antiquarks with flavor / and fractional electric charge e/, one finds: 



F 1 (x) = We 2 f (q f (x) + q f (x)), (2.19) 
A f 

F 2 {x)=2xF 1 {x). (2.20) 

The Bjorken variable x coincides with the fraction of the target light-cone momentum 
carried by the interacting quark with momentum I: 

Q 2 l-q . 

x = — ^ — = — . 2.21 

2P-q P-q y ! 



The Callan-Gross relation ( j2.20|) connecting F\ and F 2 reflects the spin- 1/2 nature of the 
quarks. 

For the spin structure functions the naive parton model gives: 



9x{x) = \ E 4 l A( lf( x ) + A lA x )] . ( 2 - 22 ) 
z / 

g 2 (x) = 0. (2.23) 

The helicity distributions Aqj(x) = ql(x) — q\{x) and Aqf(x) = qj(x) — q~f(x) involve 
the differences of quark or antiquark distributions with helicities parallel and antiparallel 
with respect to the helicity of the target nucleon. 



9 



2.3 Virtual Compton scattering 



The hadronic tensor (2.6) can be expressed as the Fourier transform of a correlation 
function of electromagnetic currents, with its expectation value taken for the nucleon 
ground state \P, S) normalized as (P', S'\P, S) = 2E p (2tt) 3 5 3 (P - P') 5 SS > 

W^(q;P,S) = — ^ J d'ze^i^SlJ^J^l^S). (2.24) 

It is related to the forward virtual Compton scattering amplitude: 

T lu ,faP,S)=i J tfze** (P,S\T '(^(2)^(0)) |P, 5), (2.25) 



where T denotes the time-ordered product. By comparison of Eqs. fl2.24"|) and ( |2.25| ) one 
finds the optical theorem: 

27cMW flu = lmT llu . (2.26) 



Consequently, nucleon structure functions can be represented in terms of virtual photon- 
nucleon helicity amplitudes, 

AHH,h>H> = e 2 e%T^(H, H') e v h . (2.27) 

Here and ey are the polarization vectors of the incoming and scattered photon with 
helicities h and h', respectively. They have values +1,-1,0 (abbreviated as +,— ,0). 
Helicities of the initial and final nucleon are denoted by H and H' . Their values are ±1, 
symbolically denoted by |, j. Choosing the z-axis in space to coincide with q/\q\, the 
direction of the propagating virtual photon, and quantizing the angular momentum of 
the target and photon along this axis yields the following relations: 



Fi = ^ QmA+ u+l + Im^l +T , +T ) , (2.28) 
F2 = QmA+ u+i + Im^ +Ti+T + 21m At,o T ) , (2.29) 

where k — 1 + (2Mx/Q) 2 . For the spin-dependent structure functions one finds: 



3i = ^( [ui^+i.+ i - | + \/2(,; - I) lmA+ ) . (2.30) 
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92 = -r^Y~ ( Im ^+T,+T - Im *4+i,+| + / 2 N : Im^+|,OT | • ( 2 - 31 ) 
^e 2 K y ^J2(k - 1) 

In the scaling limit the structure functions F±, F 2 and g\ are determined by helicity 
conserving amplitudes. It is therefore possible to express them through virtual photon- 
nucleon cross sections defined as: 

° hU = 2MK lmAhH,hH: (y2 ' 32 ^ 

with the virtual photon flux K = (2P ■ q — Q 2 )/2M. For example, the structure function 
F 2 reads: 

where the longitudinal and transverse photon-nucleon cross sections o~l,t(v, Q 2 ) are given 
by: 

a L = ^{a 01 + a ol ), (2.34) 

Ot = ^ (c^+T + <r+i + <M + ■ ( 2 - 35 ) 
An interesting quantity is their ratio: 
^ = F 2 (l + (2M,/Qff_ 1 

In the simple parton model the Callan-Gross relation (|2.20|) implies R = as Q 2 — > 00. 
Due to their interaction with gluons, quarks receive momentum components transverse to 
the photon direction. Then they can absorb also longitudinally polarized photons. This 
leads to R 7^ 0. 



2.4 QCD-improved parton model 



Nucleon structure functions systematically exhibit a weak Q 2 -dependence, even at large 
Q 2 . These scaling violations can be described within the framework of the QCD-improved 
parton model which incorporates the interaction between quarks and gluons in the nucleon 
in a perturbative way (see e.g. |||7|ii)- The scale at which this interaction is resolved is 
determined by the momentum transfer. The Q 2 -dependence of parton distributions, e.g. 
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F 2 (x,Q 2 ) =J2 e } x (lf( x iQ 2 ) + Vf( x iQ 



Sifa Q 2 ) = lH 4 [Ag/(x, Q 2 ) + Ag7(x, Q 2 



[2.37) 
'2.38) 



is described by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equa- 
tions. They are different for flavor non-singlet and singlet distribution functions. Typical 
examples of non-singlet combinations are the difference of quark and antiquark distribu- 
tion functions, or the difference of up and down quark distributions. The difference of 
the proton and neutron structure function, F$ — F£, also behaves as a flavor non-singlet, 
whereas the deuteron structure function F$ is an almost pure flavor singlet combination. 
For the flavor non-singlet quark distribution, g NS , and the flavor-singlet quark and gluon 
distributions, q and g, the DGLAP evolution equations read as follows: 



d 



dlnQ 2 



dq m (x,Q 2 
dhiQ 2 



q s (x,Q 
g(x,Q 



2tt J y H \y) 




f g s (y,Q 2 ' 
g(y,Q' 



(2.39) 
(2.40) 



Here a s (Q 2 ) is the running QCD coupling strength. The splitting function P qq (x/y) de- 
termines the probability for a quark to radiate a gluon such that the quark momentum is 
reduced by a fraction x/y. Similar interpretations hold for the remaining splitting func- 
tions. For further details we refer the reader to one of the many textbooks on applications 
of QCD, e.g. 



2.5 Light-cone dominance of deep-inelastic scattering 



The QCD analysis of deep-inelastic scattering has generated its own terminology and 
specialized jargon. In this section we summarize some of the basic notions. The general 
framework is Wilson's operator product expansion applied to the current-current correla- 
tion function. A detailed investigation reveals that the hadronic tensor 



2n MW, 



Im 



# ze *» ( P \T(J^z)J u (0))\P) 



(2.41) 



at Q 2 

~2 _ 



oo but fixed Bjorken x, is dominated by contributions from near the light-cone, 
— z 2 ~ The operator product expansion makes use of this fact by 
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expanding the time-ordered product of currents around the singularity at z 2 = 0: 

oo 

T(J(z)J(0)) ~£c£(*V) z^-'-z^O^in 2 ), (2.42) 

n=0 



where the Ml ... Mn are local operators involving quark and gluon fields. The coefficient 
functions c n are singular at z 2 = 0. They are grouped according to the order of their 
singularity. Both the operators O and the c-number coefficient functions c n depend on 
the renormalization point fj, 2 . 

The operators O can be organized according to the irreducible representation of the 
Lorentz group to which they belong. Each operator has a characteristic dimensionality, 
d, in powers of mass or momentum. For example, the symmetric traceless Lorentz tensors 
of rank n with minimum dimensionality are the operators 



^x-Mn = {^7w^» • "iD^} s , (2.43) 
£>l-» n = {G^JD, 2 ■ ■ ■ iD IM _ 1 Gl n } s , (2.44) 

local bilinears in the quark field ip and the gluon field tensor G^, with any number of 
gauge-covariant derivatives inserted between them. The brackets {}$ indicate sym- 
metrization with respect to Lorentz indices and subtraction of trace terms. The operators 
O q and O 9 have dimensionality d = 3 + (n — 1) = 4 + (n — 2) = 2 + n. The difference 
t = d — n is called "twist" (r = 2 in our example), a useful bookkeeping device to classify 
the light-cone (z 2 — > 0) singularity of the coefficient function c n . Comparing dimensions 
in Eq.(|2.42|) one finds that, for each given operator on the right-hand side, the coefficient 



behaves as c n ~ {l/z 2 )^ 2dj T ^ 2 when z 2 — > 0, where dj = 3 is the dimensionality of each 
of the currents on the left-hand side of Eq.( p.42j ). 

Matrix elements of the operators O between nucleon states are of genuinely non-perturbative 
origin. For spin-averaged quantities they must be of the form 



(P\Ol^ n {p 2 )\P) = al{p 2 )P„ 1 ...P lln , (2.45) 
(P\0° n ( f i 2 )\P) = ai( f i 2 )P, 1 ...P IM , (2.46) 

since Lorentz-covariant tensorial functions of the nucleon four-momentum P M , with P 2 = 
M 2 fixed, are proportional to the symmetric tensors P w ■ • ■ P Mn . Trace terms have been 
subtracted in Eqs.( p.45| , 2.46|) . The constants and a 9 n are fixed at a given renormalization 



scale \x and represent the non-perturbative quark and gluon dynamics of the nucleon. 

We can now make contact with observables. Since the Fourier transform of 
(P\T(J tl (z)J u (0)\P) is proportional to the forward virtual Compton scattering amplitude 
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and its imaginary part determines the structure functions F\ and F2, it is clear that the a n 
represent moments of those structure functions, with Q 2 -dependent coefficients. Consider 
as an example the structure function Fi(x, Q 2 ) in the flavor singlet channel. One finds, 



1 

J dxx^Fxix, Q 2 ) = C* n (Q 2 ; /i 2 ) <(/i 2 ) + C°(Q 2 ; /i 2 ) <(/i 2 



(2.47) 



where crossing symmetry implies a restriction to even orders n = 2, 4, . . .. The momentum 
space coefficient functions C n (Q 2 ;iJL 2 ) are related to the c-number functions c n (z 2 ; /j?) of 
Eq. fl2.42|) by Fourier transformation. The important point is that the C n (Q 2 ; fi 2 ) can be 
calculated perturbatively at large Q 2 . Their Q 2 -dependence is determined by renormal- 
ization group equations equivalent to the DGLAP equations in ( |2.40| ). 



It is often useful to express the structure functions in a factorized form, by a convolution 
of "hard" effective cross sections a q and a g for the scattering of the virtual photon from 
quarks and gluons in the nucleon, with "soft" quark and gluon distributions of the target. 
For example, 



dy 



{v q (y, Q 2 ; y" 2 ) [qfa/ y, fi 2 ) + q(x/ y, /i 2 ) 

+ °g{y,Q 2 ;n 2 )g{x/y^ 2 )} ■ 



(2.48) 



(Here we have generically used only one quark flavor with unit electric charge.) The 
perturbatively calculable functions C n then find a simple physical interpretation in terms 
of moments of the "hard" cross sections: 



Cr(QV) = / dxx^a^Q 2 ;^ 



(2.49) 



while the quark and gluon distributions are related to the "soft" reduced matrix elements 
(SHI by 



dxx 



n—l 



q(x,fi 2 ) + (-l) n q(x,fi 2 



1 

;(/i 2 )= J dxx n ~ l g{x^ 2 ), 



(2.50) 
(2.51) 



where Eq. (|2.51|) holds only for even n. To lowest (zeroth) order in the running coupling 
strength a s ([i 2 ), the "hard" cross sections are simply a q ~ 5(1 — x) and a g = 0, so that 
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only quarks contribute to F\. Gluons first enter at order a s . We mention that, in general, 
the representation of a given structure function in terms of separate quark and gluon 
contributions is a matter of definition. It is unique only in leading order and depends 
on the renormalization scheme at higher orders in a s H|7]||. The measured structure 
functions are, of course, free of such ambiguities. 



2. 6 Facts about free nucleon structure functions 



In this section we briefly review the present experimental status on free nucleon structure 
functions as measured in deep-inelastic lepton scattering. We focus on those aspects which 
are of direct relevance for our further discussion of nuclear deep-inelastic scattering. 



2. 6. 1 Spin independent structure functions 



Unpolarized deep-inelastic scattering has been explored in recent years over a wide kine- 
matic range in fixed target experiments at CERN, FNAL and SLAC, and at the HERA 
collider at DESY. Reviews can be found e.g. in Refs.[p,10,ll . 



The proton structure function F£ 



Accurate F£ data are available from fixed target measurements at SLAC, at CERN 
(BCDMS, NMC) and at Fermilab (E665). They cover the kinematic range 10 -3 < x < 0.8 
and 0.2 GeV 2 < Q 2 < 260 GeV 2 |J. Due to experimental constraints fixed target studies 
at small x are possible only at low Q 2 . For example, in the E665 measurements at Fer- 
milab the smallest values of the Bjorken variable, x ~ 0.8 ■ 10~ 3 , are measured typically 
at Q 2 ~ 0.2 GeV 2 
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This is different at the HERA collider where the kinematic range 
3 • 10 -6 < x < 0.5 and 0.16 GeV 2 < Q 2 < 5000 GeV 2 is explored. In these experiments 
the region of small x is accessible also at large Q 2 . 



The data summarized in Figures Til and 
erences see e.g. P, |IT||1"T1| ): 



display several important features (for ref- 



At small Bjorken-x (x <C 0.1) but large Q 2 a strong increase of F$ with decreasing 
x has been found at HERA. This behavior is commonly interpreted in terms of the 
dominant role of gluons at small x, the density of which rises strongly with decreasing 
x. This increase becomes weaker at low Q 2 . Here only a minor x-dependence has been 
observed in fixed target experiments, which is nevertheless enhanced at very small 
i < 0.1 as recently explored at HERA ||14,17|. Note that a rise of F 2 P with decreasing 
x reflects a growing virtual photon-proton cross section as the photon-nucleon center- 
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Fig. 2.2. The proton structure function as a function of x for various Q 2 . The data are taken 
from HI ZEUS |OM||, E665 O, NMC M, SLAC @ 5 and BCDMS M. 



of-mass energy W = ^/s increases. For example, at Q 2 ~ 100 GeV 2 one observes a 



characteristic behavior [13 



7*P 



(W 2 ) A , withAw0.3. (2.52) 



For the real photon-nucleon cross section at high energies, on the other hand, one has 
A fa 0.08 f2~lf . The dynamical origin of the observed variation of the energy depen- 
dence of <7 7 *n with Q 2 is an important issue of ongoing investigations (see for example 
Refs.|JT|). 



Hadron-hadron interaction cross sections have an energy dependence similar to that ob- 
served in photon-nucleon scattering. It is often parametrized using Regge phenomenol- 



ogy pl|, P2l. In Regge theory the dependence of cross sections on the center-of-mass 
energy is determined by the t-channel exchange of families of particles permitted by 
the conservation of all relevant quantum numbers. Each group of particles is charac- 
terized by a Regge trajectory, a(t) « a(0) + ta', which relates their spin with their 
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invariant mass. The resulting dependence of hadron-hadron total cross sections on the 
squared center-of-mass energy s is: 



<T m ~s a(0yi . (2.53) 

The rising hadron-hadron cross sections at high energies are well described by the so- 
called pomeron exchange. It corresponds to multi-gluon exchange with vacuum quantum 
numbers and it is characterized by the trajectory [^l| ,p3]] 

a F (t) « a w (0) +ta' F & 1.08 + 1 0.26 GeV~ 2 . (2.54) 



Note that the fast growth of the interaction cross section fl2.53|) with energy as implied 



by Eq.( p.54j ) cannot persist up to arbitrarily high energies because of limitations im- 



posed by unitarity. At asymptotic energies the Froissart bound does not permit total 
hadronic cross sections to rise faster than (lns/so) 2 with some constant scale Sq 



The slow decrease of hadron cross sections at moderate energies is described by an 
exchange made up from a set of reggeons which lie on the approximately degenerate 
trajectory plj2"4"| 



«n(*) -cvW+^m ~0.5 + t0.9GeV" 2 , (2.55) 
and which carry the quantum numbers of the p, u, a2 and f% mesons, respectively. 



F p 
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Fig. 2.3. The Q 2 -dependence of the proton structure function for x < 0.1. The data are taken 
from ZEUS E665 @, and NMC If. 
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At large energies these so-called subleading contributions are exceeded by pomeron 
exchange 



At small values of Q 2 (i.e. Q 2 < 1 GeV 2 ) the structure function FJ 5 drops. This is quite 
natural in view of the fact that F 2 P has to vanish linearly with Q 2 in the limit Q 2 — >• 
as a consequence of current conservation (see e.g. [pTl ). Bjorken scaling must break 




down in this kinematic regime. In particular, at small x < 0.1 or large photon energy, 
v > 5 GeV, vector meson dominance is expected to play an important role. It describes 
(virtual) photon-nucleon scattering via the interaction of vector meson fluctuations of 
the photon. The contribution to -Fj 3 from the three lightest vector mesons reads (see 
e.g. HI): 

F^ VMD \x,Q 2 ) = - X ( — ) ( -r-TF ) -y, <^c> 



The sum is taken over p, lu, and mesons with their invariant masses my. The vector 
meson-proton cross sections are denoted by ay p . The vector meson-photon coupling con- 
stants gy can be deduced from electron-positron annihilation into those vector mesons. 
One observes that p^ VMD ^ ^ Q 2 at small Q 2 . At large Q 2 , however, the vector meson 
contribution ( |2.56| ) vanishes as 1/Q 2 - Then the scattering from parton constituents in 
the target takes over and leads to Bjorken scaling. 

Finally, at large values of x one observes a rapid decrease of the structure function. This 
can be understood within the framework of perturbative QCD. In the limit x — > 1, a 
single valence quark struck by the virtual photon carries all of the nucleon momentum. 
The only way for such a configuration to evolve from a bound state wave function 
which is centered around low parton momenta, is through the exchange of hard gluons. 
A perturbative description of this process leads to F£(x — > 1) ~ (1 — a;) 3 [E6f. 



2.6.2 The ratio of longitudinal and transverse cross sections 

Extracting the structure function F-i from lepton scattering data requires information on 
the ratio of the total cross section for longitudinally and transversely polarized photons, 
R — ^ from Eq. ( |2.36[ ). Previous data from SLAC and CERN cover the region 0.1 < x < 



(Tj 1 

0.9 and 0.6 GeV 2 < Q 2 < 80 GeV 2 gTJ. In this region R is small. New data from the NMC 
collaboration are available for 0.002 < x < 0.12 A rise of R with decreasing x has 
been observed as shown in Fig. |2.4| . This behavior can be understood within the framework 
of perturbative QCD pH] . Helicity conservation implies that a high-Q 2 longitudinally 
polarized photon cannot be absorbed by a quark moving in longitudinal direction: a non- 
zero transverse momentum is necessary for this process to occur. In the QCD-improved 
parton model such transverse quark momenta result from gluon bremsstrahlung which is 
important for low parton momenta, i.e. at small x. Further studies of R in the domain of 
small x are currently performed at HERA. A first analysis gives R ~ 0.5 at x = 2.4 • 10~ 4 
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Fig. 2.4. The ratio R = ol/ctt as a function of x. The data are taken from NMC [18], BCDMS 

mum, and cdhsw 



and Q 2 = 15 GeV 2 p . 



2.6.3 Spin dependent structure functions 

In recent years polarized deep-inelastic scattering experiments have become a major ac- 
tivity at all high-energy lepton beam facilities. They aim primarily at the exploration of 
the spin structure of nucleons. Detailed investigations have been carried out at CERN 
(SMC), SLAC (E142/143/154/155) and DESY (HERMES). For references see [33 - 41]. 

While the proton spin structure functions gf and g% have been measured directly using 
hydrogen targets, neutron structure functions have been extracted from measurements 
using deuterons and 3 He targets with corrections for nuclear effects. In the data analy- 
sis such corrections have commonly been done in terms of effective proton and neutron 
polarizations obtained from realistic deuteron and 3 He wave functions. They account for 
the fact that bound nucleons carry orbital angular momenta. As a consequence their po- 
larization vectors need not be aligned with the total polarization of the target. At the 
present level of accuracy the use of effective polarizations turns out to be a reasonable 
approximation as discussed at length in Section |7|. 



In Fig.gj| we show a collection of data for g\. The behavior of the proton, deuteron and 
neutron structure functions turns out to be quite different, especially in the region of 
small x. This is in contrast to the unpolarized case where proton and neutron structure 
functions show a qualitatively similar behavior. 
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Fig. 2.5. A compilation of data of the proton, deuteron, and neutron spin structure functions g\ 
from Refs.|33|,36l,37,38|, 40|,42,43|. (We thank U. Stoesslein for the preparation of this figure.) 



The moments 



dxgi' n (x, Q 2 



(2.57) 



of the proton and neutron spin structure functions are of fundamental importance. They 
can be decomposed in terms of proton matrix elements of SU(3) axial currents, as follows 



(for a review see e.g. [fl4] , fi5 



1 U 



1 



(2.58) 
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with the axial vector matrix elements: 

MS,Aq a = (P,S\^ l5 ^ij\P,S), (2.59) 



where ip = (ip u ,'4 ) d,' l l>s) is the quark field. Here A a (a = 1, ...,8) denote SU(3) flavor 
matrices and the singlet Ao is the 3x3 unit matrix. In Eq. (|2.58 ) and below we suppress 
QCD corrections which are currently known up to order o? s . Current algebra and isospin 
symmetry equate the non-singlet matrix element Ag 3 = Au — Ad with the axial vec- 
tor coupling constant = 1.26 measured in neutron /3-decay. One thus arrives at the 
fundamental Bjorken sum rule: 



Tl - rj = ~Ag 3 = \g A . (2.60) 
6 6 

Furthermore, assuming SU(3) flavor symmetry, Ag§ = (Au + Ad — 2As)/v / 3 is determined 
by hyperon /5-decays. The non-singlet matrix elements Ag 3j8 involve conserved currents, 
hence they are scale independent. This is different for the singlet term Ago = Am + Ad + As 
which receives a Q 2 -dependence through the QCD axial anomaly. Note that in next-to- 
leading order both quarks and gluons contribute to Ag . However, the detailed separation 
into quark and gluon parts depends on the factorization scheme used to separate pertur- 
bative and non-perturbative parts of the spin-dependent cross section. 

An evaluation of the structure function moments T± ,n (Q 2 ) from Eq. fl2.57p requires knowl- 
edge of gi in the entire interval < x < 1. Since measurements cover only a limited 
kinematic range, data for gi have to be extrapolated to x — > and x — *■ 1. The large- x 
extrapolation is not critical since g\ becomes small and ultimately vanishes as x — > 1. The 
situation at small x is, however, not yet well understood (for a review and references see 
Ref.||10||). The common approach is to assume Regge behavior which implies that g\ ~ x a 
with < a < 0.5 for x — ► 0. 

A status review of the analysis of spin structure functions and their moments can be 
found in Refs. P3L|I| . All current studies arrive at the conclusion that the flavor singlet 
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contribution to the nucleon spin is small. At Q 2 = 1 GeV 2 one finds (in the AB scheme) 
Ago = 0.23 ± 0.07(stat) ± 0.19(sys). (2.61) 



This would imply that only about one third of the nucleon spin is carried by the quark 
spins alone. The missing two thirds probably involve gluon spin contributions and orbital 
angular momentum of quark, antiquark and gluon constituents. Finally we note that the 
Bjorken sum rule (p.60|) , with inclusion of QCD corrections, is fulfilled at the 5% level 

0- 
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Fig. 2.6. Diffractive scattering from a nucleoli. 



Measurements of g 2 have been performed at CERN [46[] and SLAC |35|j47|1 . For the neu- 
tron case 3 He and deuteron j|7|| targets have been used again. Within large ex- 
perimental errors the data for g 2 are consistent with the twist-2 prediction g 2 (x,Q 2 ) = 
Jl % (1 - 5(1 - x/y)) gi (y,QZ) of Ref.@. 



2.6.4 Diffraction 

A subclass of photon-nucleon processes, namely diffractive lepto- and photoproduction, 
plays a prominent role also in the interaction of real and virtual photons with com- 
plex nuclei at high energies. We focus here on so-called single diffractive processes. They 
are characterized by the proton emerging intact and well separated in rapidity from the 
hadronic state X produced in the dissociation of the (virtual) photon (see Fig. j2.6|) : 



7 



» 



+ p -> X + p'. 



(2.62) 



As in diffractive hadron-hadron collisions such processes are important at small momen- 
tum transfer. Their cross sections drop exponentially with the squared four-momentum 
transfered by the colliding particles. Furthermore, they generally exhibit a weak energy 
dependence. 
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Diffractive leptoproduction 



In deep-inelastic scattering experiments at HERA approximately 10% of the (virtual) 
photon-proton cross section result from diffractive events (for a review see e.g. []49]). 
Their cross section is parametrized in terms of two structure functions, analogous to the 
inclusive case. One has: 



da 



Aira 2 \l — y 



dxdQ 2 dx p dt Q 4 



x 



+ 



The diffractive structure functions depend on x and Q 2 , on the squared momentum trans- 
fer t to the proton, t = (P — P') 2 = (q — px) 2 , and on the variable 



x F — 



(P-P')-q Q 2 + M x -t 



Q 2 + M? 



x 



P ■ q " Q 2 + W 2 - M 2 ~ Q 2 + W 2 ' 



(2.64) 



Here Mx is the invariant mass of the diffractively produced system X in the final state. The 
diffractive structure function, conventionally denoted by F®^ indicating its dependence 
on four kinematic variables, is directly related to the diffractive (virtual) photoproduction 
cross section. At small x one finds in analogy with Eq. (|2.33|) : 



,D(4), 



Q 2 da% f 



4n 2 a dx p dt 



(2.65) 



Most of the data have so far been obtained for the t-integrated structure function 



u 

F? 3 \x,Q 2 ;x F )= J dtF? i4 \x,Q 2 ;x p ,t). 



(2.66) 



Measurements by the HI [ 5D| , |5T | and ZEUS [ 52| , |53 ,53,55| groups cover the range 4.5 < 
Q 2 < 140 GeV 2 , 2 • 10~ 4 < x p < 0.04 and 0.02 < x/x F < 0.9. No substantial Q 2 - 
dependence of F^ 3 ^ has been found. Over most of the explored kinematic region, x ]P i 7 ^ 3 ' ) 
is either decreasing or approximately constant as a function of increasing x p . However 
at small x/x p there is a tendency for x p F2^ to increase at the highest values of x w . A 



typical collection of data is shown in Fig.|2/7 



A reasonably successful description of this behavior has been achieved within Regge phe- 
nomenology which assumes that the interaction proceeds in two steps: the emission of a 
pomeron or subleading reggeon from the proton, and the subsequent hard scattering of 
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the virtual photon from the partons in the pomeron or reggeon, respectively. This picture 
leads to a factorization of the diffractive structure function 



F 2 (x, Q ; x F , t) — /jp (Xjp , t) F 2 (x/x F , Q ) + / M (Xjp , t) F 2 (x/x F ,Q ), 



(2.67) 



where F^^ is commonly interpreted as the "structure function" of the pomeron (reggeon) 
and 



r 2ai(t)-l 
np 



(2.68) 



with z = IP, IR, denotes the pomeron (reggeon) distribution in the proton. 



The HI analysis [EDJ gives a F (0) ~ 1.2 and 0^(0) ~ 0.5. The slope parameters B 



ex 

IP(IR)' JR. 



and were taken to reproduce hadron-hadron data. While 0^(0) is found to be slightly 
larger than the value obtained from parametrizations of hadronic cross sections, a K (0) 



agrees well with the Regge phenomenology of hadron-hadron collisions 



In Fig. (|2.8| ) we show recent ZEUS data [5^] on the ratio of diffractive and total photon- 
nucleon cross sections for different regions of Mx- The data show a similar energy de- 
pendence of both total and diffractive cross sections. Furthermore, the observed Q 2 - 
dependence of the cross section ratio for different regions of Mx suggests that, as Q 2 
increases, diffractive states with large mass become important. 



ZEUS measurements |53[ have investigated the t-dependence of the diffractive leptopro- 
duction cross section. In the range 5 < Q 2 < 20 GeV 2 and 50 < W < 270 GeV the 
t-dependence of the diffractive virtual photoproduction cross section is described for 
0.07 < |t| < 0.4 GeV 2 by the exponential form, da^JJ ' jdt ~ e Bt , with B « 7 GeV 2 . 
This value is compatible with results from high-energy hadron-hadron scattering (see e.g. 
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Fig. 2.7. The diffractive structure function x v F 2 for different values of j3 = x F /x and Q 2 « 28 



GeV . Data from HI |5C| (open squares) and ZEUS |52] (solid points). 
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Fig. 2.8. ZEUS data |3^] for the ratio of diffractive and total photon-nucleon cross sections. The 
diffractive production cross section has been integrated over different intervals of Mx- 



Diffractive photoproduction 

Diffractive dissociation of real photons, 7 + N — > X + N, has been explored with fixed 
target and collider experiments. At FNAL [|5£| photon-proton center of mass energies up to 
W ~ 15 GeV were used to produce diffractive states with an invariant mass up to Mx — 5 
GeV. Recent experiments at HERA []59 i|60| , |6l] ) |6^ , |6^| were carried out at W ~ 200 GeV 



and Mx < 30 GeV. The diffractive cross section amounts to approximately 20% of the 
total photon-proton cross section. Around half of these events come from the production 
of the light vector mesons p, to and <p. This is contrary to diffractive leptoproduction at 
large Q 2 where vector meson contributions are suppressed roughly as l/Q 4 



At sufficiently large mass Mx of the diffractively produced system X, the measured cross 
sections drop approximately as 1/M|, as shown in Fig. |2.9| . This is in accordance with 
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Fig. 2.9. Differential cross section for diffractive photoproduction off nucleons from FNAL 



and HI [61] for different center of mass energies W. The curves corresponds to a Regge fit [ plj| . 

oo with M| — > oo only pomeron exchange 



Regge phenomenology. In the limit W 2 /M^ 
is important and leads to 



da*g f ^4( Qjp (o)-i) 



dM\dt 



M^ (0) 



exp 



/ W 2 



(2.69) 



with a slope parameter 5. Equation ( 2.69| ) implies that at energies W = (15 — 30) GeV 
typical for fixed target experiments at CERN and FNAL, the relative amount of diffraction 
in deep- inelastic scattering is reduced to (10 — 15)% ffilf] . Observed deviations from the 
simple behavior ( |2.69 ) have been associated with contributions involving subleading Regge 
trajectories [JB"0l,|B"l]l . 



3 Deep-inelastic scattering from nuclear systems 



3.1 Introduction and motivation 



We now enter into the central topic of this review: an exploration of new phenomena 
specific to deep-inelastic lepton scattering from nuclear (rather than free nucleon) targets. 

Nuclei represent systems with a natural, built-in length scale. The baryon density in the 
center of a typical heavy nucleus is po — 0.15 fm -3 . The average distance between two 
nucleons at this density is 

d ~ 1.9 fm. (3.1) 
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The nucleons have a momentum distribution characterized by their Fermi momentum, 

(o 2 \ !/ 3 
— po I ^1-3 far 1 ~ 0.26 GeV. (3.2) 

A high energy virtual photon which scatters from this system can expect to see two sorts 
of genuine nuclear effects: 

i) Incoherent scattering from A nucleons, but with their structure functions modified 
in the presence of the nuclear medium. Such modifications are expected to arise, for 
example, from the mean field that a nucleon experiences in the presence of other 
nucleons, and from its Fermi motion inside the nucleus; 

ii) Coherent scattering processes involving more than one nucleon at a time. Such effects 
can occur when hadronic excitations (or fluctuations) produced by the high energy 
photon propagate over distances (in the laboratory frame) which are comparable to 
or larger than the characteristic length scale d ~ 2 fm of Eq. ( |3 . 1| ) . A typical example 
of a coherence effect is shadowing. 

It turns out, as we will demonstrate, that incoherent scattering takes place primarily in 
the range 0.1 < x < 1 of the Bjorken variable. Strong coherence effects are observed at 
x < 0.1. Cooperative phenomena in which several nucleons participate can also occur at 
x > 1. (In fact, the Bjorken variable can extend, in principle, up to x < A in a nucleus 
with A nucleons.) 

The aim of this section is to prepare the facts and phenomenology of nuclear DIS. An 
important subtopic in this context deals with the deuteron. While this is not a typical 
nucleus, it serves two purposes: first, as a convenient neutron target, and secondly, as 
the simplest prototype system in which coherence effects, involving proton and neutron 
simultaneously, can be investigated quite accurately. For this purpose we need to introduce 
the hadronic tensor and structure functions for spin-1 targets as well. Once the nuclear 
structure functions are at hand we will present a survey of nuclear DIS data and give 
first, qualitative interpretations. The more detailed understanding is then developed in 
subsequent sections. 



3.2 Nuclear structure functions 



The deep-inelastic scattering cross sections for free nucleons and nuclei have basically the 
same form as given by Eq. ( |2.1|) . All information about the target and its response to the 
interaction is included in the corresponding hadronic tensor. For nuclei with spin 1/2 the 
hadronic tensor formally coincides with the one for free nucleons given in Eqs.( |2.6|J2.71 . p.8| ). 
In this case nuclei are characterized by four structure functions, F± 2 an d gfV F° r spin-0 
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targets, only the symmetric tensor ( |2.7|) with the structure functions F^ 2 is present. In the 
case of spin-1 targets the situation is more complex. Here the hadronic tensor is composed 
of eight independent structure functions |67|j68|] :H 



A | P»Pu pA | , M A 
"" + P.q 2 ' P-q 



Ma W^ v = - 9iat + 



it + 9$) 



S ■ q 
P~q 



with the Lorentz tensors: 



(3.3) 



/J.U " 



( Mj V 



(3.4) 




The tensors ( |3.4p are functions of the photon and target four-momenta q^ and P M , the 
target polarization vector £ , and the spin vector S a = —iSa^i, E 13 * S^P 6 /Ma. Furthermore 
we have used the notation ka = 1 + MlQ 2 / (P ■ q) 2 where Ma denotes the nuclear mass. 



The nuclear structure functions in Eq.( |3.3| ) depend on the Bjorken scaling variable of 
the target, xa = Q 2 /2P ■ q with < x& < 1, and on the momentum transfer Q 2 . Note 
however that these functions are frequently expressed in terms of the Bjorken variable 
of the free nucleon which is x = Q 2 /2Mv = x^M^jM in the lab frame, and which 
can extend over the interval < x < Ma/M ~ A. The first four structure functions in 
Eq. ( |3.3|) are proportional to Lorentz structures already present in the case of free nucleons 
( p.6|) or spin- 1/2 nuclei. The new structure functions can be measured in the scattering 
of unpolarized leptons from polarized targets. By analogy with the Callan-Gross relation 
Q2.20| ) one finds b 2 = 2x A bi in the scaling limit. The deuteron structure function bf is 



subject of investigations at HERMES |69] 



For spin-1/2 nuclei the relations between nuclear structure functions and photon-nucleus 
helicity amplitudes A1 Hh , H i are analogous to the ones for free nucleons in Eqs.( |2.28H2.31p . 
For spin-1 targets with helicity H,H' = +, — , one obtains |67|,|68|,[70 



2 We omit terms proportional to or q v which do not contribute to the cross section ( 2.1 
due to electromagnetic gauge invariance. 
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F A = ^-L (i m ^;a ++ + Im ^;a + _ + Im aX a +o ^ (3 5) 

F " A= (lm^; + A ++ + Im^t A + - + Im^o + 2Im„4£. A 0+ + Im^o.oo), (3-6) 

01=4^ ( Im ^I-V " Im^T* + A ++ + ^I^l(Im^; A 0+ + Im^l A 00 )) , (3.7) 

\mAXt, ++ ~ Im^t A + _ + -^=(Im^oV + ^X-,oo)) , (3-8) 

bf = (lm^* + A ++ + Im^tV - 2Im^ot + o) • (3-9) 

Corresponding relations for the remaining structure functions can be found for example 
in Ref.@[7§. 



3. 3 Data on nuclear structure functions 



In this section we summarize the existing experimental information on nuclear effects in 
structure functions. Their systematic investigation for light and heavy nuclei has been 
carried out so far only in unpolarized scattering experiments. Most of the data come 
from deep-inelastic lepton scattering. Modifications of nuclear parton distributions have 
also been studied in other high-energy processes. We mention, in particular, heavy quark 
production and Drell-Yan experiments. 



3. 3. 1 Nuclear effects m F 2 A 

Experiments on deep-inelastic scattering from nuclei are reviewed in @||. For a discussion 
of the data it is convenient to use structure functions which depend on the Bjorken scaling 
variable for a free nucleon, x = Q 2 / iflMv). In charged lepton scattering from unpolarized 
nuclear targets these structure functions are defined by the differential cross section per 
nucleon: 



2^A 



(Pa 



Ana 2 



dxdQ 2 Q 4 



i-y 



Mxy\ Ff{x,Q 2 



2E 



1 



x 



+ y 2 F 1 A (x,Q 2 ) 



(3.10) 



Some time ago the EMC collaboration discovered that the structure function F% for iron 



differs substantially from the corresponding deuteron structure function [74], far beyond 



trivial Fermi motion corrections. Since then many experiments dedicated to a study of 
nuclear effects in unpolarized deep-inelastic scattering have been carried out at CERN, 
SLAC and FNAL. The primary aim was to explore the difference of nuclear and deuterium 
structure functions. 
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Fig. 3.1. The structure function ratio F^/F^ for 40 Ca and 56 Fe. The data are taken from NMC 
[7|, SLAC p|, and BCDMS WB. 



Figure |3.1| presents a compilation of data for the structure function ratio F^/F^ over 
the range < x < 1. Here F^ is the structure function per nucleon of a nucleus with 
mass number A, and _F 2 d refers to deuterium. In the absence of nuclear effects the ratios 
F2/F2 are thus normalized to one. Neglecting small nuclear effects in the deuteron, F^ can 
approximately stand for the isospin averaged nucleon structure function, F^. However, the 
more detailed analysis must include two-nucleon effects in the deuteron. Several distinct 
regions with characteristic nuclear effects can be identified: at x < 0.1 one observes a 
systematic reduction of F^/F^, the so-called nuclear shadowing. A small enhancement is 
seen at 0.1 < x < 0.2. The dip at 0.3 < x < 0.8 is often referred to as the traditional 
"EMC effect". For x > 0.8 the observed enhancement of the nuclear structure function is 
associated with nuclear Fermi motion. Finally, note again that nuclear structure functions 
can extend beyond x — 1, the kinematic limit for scattering from free nucleons. 

• Shadowing region 

Measurements of E665 [00,0 at Fermilab and NMC |7I]J7|,|7S|,|SD|,0g at CERN 



provide detailed and systematic information about the x- and A-dependence of the 
structure function ratios F^/F^. Nuclear targets ranging from He to Pb have been 
used. A sample of data for several nuclei is shown in Fig. |3.2j . While most experiments 
cover the region x > 10~ 4 , the E665 collaboration provides data for F^/F^ [[76] down 



to x ~ 2 • 10 . Given the kinematic constraints in fixed target experiments, the small 
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x-region has been explored at low Q 2 only. For example, at x ~ 5 ■ 10 3 the typical 
momentum transfers are Q 2 ~ 1 GeV 2 |75|| . At extremely small values, x ~ 6 ■ 10 -5 , 
one has Q 2 ~ 0.03 GeV 2 [[f6]. 

In the region 5 ■ 10~ 3 < x < 0.1 the structure function ratios systematically decrease 
with decreasing x. At still smaller x one enters the range of small momentum transfers, 
Q 2 ~ 0.5 GeV 2 , approaching the limit of high-energy photon-nucleus interactions with 
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Fig. 3.2. (a) NMC data |7| for the structure function ratio F 2 A /F 2 d for 4 He, 12 C, and 4U Ca. (b) 
The ratio F 2 A /F 2 d for 6 Li, 12 C Wf$, and 131 Xe WE\. 
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Fig. 3.3. The shadowing ratio A e jf /A = o" 7 a/ 'Act^n for 63 Cu as a function of the photon energy 
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Fig. 3.5. The structure function ratio F 2 d /F 2 p . Data from E665 Q and NMC p| 



real photons. As an example we show in Fig. |3.3| data on shadowing for real photon 
scattering from 63 Cu. 



Shadowing systematically increases with the nuclear mass number A. For example, at 
x ~ 0.01 one finds F 2 A /F 2 C ~ A a ~ l with a ~ 0.95 [|ST]| . A similar behavior has been 
observed in high-energy photonuclear cross sections |Mj : their A-dependence is roughly 
<7 7 a ~ v4 92 o" 7 n where <t 7 n is the free photon-nucleon cross section averaged over proton 
and neutron. 

The shadowing effect depends only weakly on the momentum transfer Q 2 . The most 
precise investigation of this issue has been performed for the ratio of Sn and carbon 
structure functions presented in Fig.|3"^| p2| . It reveals that shadowing decreases at 
most linearly with \nQ 2 for x < 0.1. The rate of this decrease becomes smaller with 
rising x. At x > 0.1 no significant Q 2 -dependence of F^/F^ is found. 



Shadowing has also been observed in deep-inelastic scattering from deuterium 
lightest and most weakly bound nucleus. In Fig.|3.5| we show data from E665 |9~I 



the 
and 



NMC [0 for the ratio F 2 / F% of the deuteron and proton structure functions. At 



x < 0.1 this ratio is systematically smaller than one. 
Enhancement region 

The NMC data have established a small but statistically significant enhancement of 
the structure function ratio at 0.1 < x < 0.2. The observed enhancement is of the 
order of a few percent. For carbon and calcium it amounts to typically 2% [jg^l . The 
most precise measurement of this enhancement has been obtained for F 2 Sn / 'F 2 C shown 
in Fig j3.4j . Within the accuracy of the data no significant Q 2 -dependence of this effect 



has been found in this region. 
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Region of "EMC effect" 

The region of intermediate 0.2 < x < 0.8 has been explored extensively at CERN 
and SLAC. In the range 2 GeV 2 < Q 2 < 15 GeV 2 , data were taken by the E139 



collaboration [|72] for a large sample of nuclear targets between deuterium and gold. 
The measured structure function ratios decrease with rising x and have a minimum at 
x ~ 0.6. The magnitude of this depletion grows approximately logarithmically with the 
nuclear mass number. The observed effect agrees well with data for the ratios of iron 
and nitrogen to deuterium structure functions from BCDMS taken at large Q 2 values, 



14 GeV < Q 2 < 200 GeV 2 ppl. These data imply that a strong Q 2 -dependence of 



the structure function ratios is excluded. 
Fermi motion region 

At x > 0.8 the structure function ratios rise above unity [72], but experimental infor- 
mation is rather scarce. The free nucleon structure function F 2 N is known to drop as 
(1 — a;) 3 when approaching its kinematic limit at x — 1. Clearly, even minor nuclear 
effects appear artificially enhanced in this kinematic range when presented in the form 
of the ratio F 2 A /F 2 N . 

The region x > 1 

Data at large Bjorken x and large momentum transfer, 0.7 < x < 1.3 and 50 GeV 2 < 



Q 2 < 200 GeV 2 , have been taken for carbon and iron by the BCDMS and CCFR 
|p5| collaborations, respectively. The results disagree with model calculations at x ~ 1 
which account for Fermi motion effects only. For Q 2 < 10 GeV 2 data have been taken 
at SLAC for various nuclei |96| , |97|j98| , |99| , |l 0|| . Both quasielastic scattering from nucleons 
as well as inelastic scattering turns out to be important here. 



3.4 Moments of nuclear structure functions 




represents the difference of the integrated momentum fraction carried by quarks in a 
nucleus relative to that for deuterium. The constant = (AM /Ma — 2M/M^) corrects 
for the different mass defects of bound systems. Note that in Eq. ( |3.11| ) we have omitted 
QCD target mass corrections ||101|| . An analysis based on the NMC J7^[ and SLAC [72 



data has been performed for He, C and Ca ||]. In the kinematic range covered by these 
experiments, 3.5 • 10~ 3 < x < 0.8, the difference of the structure function moments 
M2 — M 2 turns out to be compatible with zero. Together with the well established result 
of the momentum sum rule for the proton H, one can therefore conclude that, within the 
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Fig. 3.6. NMC data |2j for R Sn - R G as a function of x for Q~ « 10 GeV 2 . The average values 
for R Ca - R c pi, and i? Au - # Fe Jl03| are also shown. 



accuracy of present data, quarks carry about half of the total momentum, in nuclei as 
well as in free nucleons. 



3. 5 Ratios of longitudinal and transverse cross sections 



Investigations of the differences between the longitudinal-to-transverse cross section ratios 
R = o'l/ 't Q2.36 ) for different nuclei have been performed at SLAC for moderate and 
large values of x, while the region of small x has been investigated by NMC. The difference 
R d — R p is found to be compatible with zero P?|fJ^JTT?2| . Similar observations have been 
made for heavier targets pj [102iTOF03|JI05| . In Fig.|J we show NMC data [|2j for 
R Sn — R c as a function of x for an average Q 2 of about 10 GeV 2 . In addition we present 
the average values from the NMC measurement for i? Ca — R c ||104|| , and for R Au — R Fe from 
SLAC E140 ||103|| . All measurements are consistent with only marginal nuclear dependence 
of R. This implies that nuclear effects influence both structure functions F\ and F2 in a 
similar way, and that the ratio of nuclear cross sections directly measures the ratio of the 
corresponding structure functions F 2 . 
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3. 6 Other measurements of nuclear parton distributions 



Nuclear deep-inelastic scattering is sensitive only to the sum of valence and sea quark 
distributions (see e.g. Eq.( |2.37j )), weighted by their respective electric charges. In order to 
separate nuclear effects in the valence and sea quark sectors, and directly measure nuclear 
gluon distributions, other types of processes are required which we briefly summarize in 
the following. 



3.6.1 Drell-Yan lepton pair production 



In the Drell-Yan production of lepton pairs (mostly in hadron- nucleus collisions, 

the underlying partonic sub-process is the annihilation of a quark and antiquark from 
beam and target into a time-like high energy photon, which subsequently converts into 
the observed dilepton. The Drell-Yan cross section reads (see e.g. 



106 



d 2 a 
dxrdxs 



Ana 2 
9mf 



K E4 [qf(xB,Q 2 )qJ(xT,Q 2 ) +qf(x B ,Q 2 )qT(x T ,Q 2 )} , (3.12) 



where mi is the invariant mass of the produced lepton pair. The flavor dependent quark 
distributions of the projectile and target are denoted by qf and qj, respectively. Seen 
from the center-of-mass frame the active quarks carry fractions xb and xt of the beam 
and target momenta. They are determined by the momentum component q^ of the pro- 
duced dilepton parallel to the beam, its invariant mass m ; and the squared center-of-mass 
energy s: 

xtXb = — , %f = -~f = x b — Xt- (3.13) 

S \/S 



Higher order QCD corrections to the production cross section ( |3.12|) turn out to be sig- 
nificant. They are absorbed in the so-called "i^-factor" and effectively double the leading 
order cross section. 



The E772 experiment at FNAL 107 ] has investigated Drell-Yan dilepton production in 
proton-nucleus collisions at s = 1600 GeV 2 . At xp > 0.2 the production process is domi- 
nated by the annihilation of projectile quarks with target antiquarks. Outside the domain 
of quarkonium resonances, i.e. for 4 GeV < mi < 9 GeV and mi > 11 GeV, this experi- 



ment explores possible modifications of nuclear sea quark distributions. In Fig. |3.7| we show 
ratios of dimuon yields for nuclear targets and deuterium taken at xf > 0. At iy > 0.1 
no significant nuclear effects have been observed within admittedly large experimental 
errors. This indicates the absence of strong modifications of nuclear sea quark distribu- 
tions, as compared to those of free nucleons. At xt < 0.1, on the other hand, the observed 
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Fig. 3.7. Drell-Yan dimuon yields per nucleon for 40 Ca and 184 W as a function of xt for xp > 



attenuation for heavy nuclei implies a substantial reduction of nuclear sea quarks, in qual- 
itative agreement with the shadowing effects observed in nuclear deep-inelastic scattering 
at x < 0.1. The detailed comparison of shadowing in Drell-Yan versus DIS requires, of 
course, a careful separation of valence and sea quark effects as well as their Q 2 evolution 

im. 



3.6.2 Lepton-induced production of heavy quarks 

The intrinsic heavy-quark (c- or 6-quark) distributions in nucleons or nuclei are expected 
to be very small. Inelastic heavy-quark production is therefore assumed to receive its major 
contributions from photon-gluon fusion, i.e. the coupling of the exchanged virtual photon 
to a heavy quark pair which is attached to a gluon out of the target. This mechanism is a 
basic ingredient of the so-called color-singlet model ||109||. In this model the cross section 



for heavy quark pair production is proportional to the gluon distribution of the target. A 
comparison of these cross sections for nucleons and nuclei can then be directly translated 
into a difference of the corresponding gluon distributions. 



In this context NMC has analyzed J /if) production data from Sn and carbon nuclei ||110| . 
The average ratio of the corresponding inelastic J / if) production cross sections was found 
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slightly larger than one: 

a( 7 * + C ^ J/ip + X) v 1 

Within the color singlet model this implies an enhancement by about 10% of the gluon 
distribution in Sn as compared to carbon in the region x ~ 0.1, though with large errors. 



3.6.3 Neutrino scattering from nuclei 

Deep-inelastic neutrino scattering permits one to separate valence and sea quark dis- 
tributions. It is therefore a promising tool to investigate modifications of the different 
components of quark distributions in nuclei. The observed nuclear effects in neutrino ex- 
periments are qualitatively similar to the results from charged lepton scattering discussed 
previously | |lll|jll2| , |113| , |114| | , although their statistical significance is poor, given the large 
experimental uncertainties. 



4 Space-time description of deep-inelastic scattering 



So far our picture of deep-inelastic scattering has been developed in momentum space. 
The partonic interpretation of structure functions is particularly transparent in the infinite 
momentum frame in which the nucleon (or nucleus) moves with (longitudinal) momentum 
P — > oo. In this frame the Bjorken variable x has a simple meaning as the fraction of the 
nucleon momentum carried by a parton when it is struck by the virtual photon.^] 

For an investigation of nuclear effects in DIS the infinite momentum frame is not always 
optimal. Instead, it is often preferable to describe the scattering process in the laboratory 
frame where the target is at rest. Only in that frame the detailed knowledge about nuclear 
structure in terms of many-body wave functions, meson exchange currents etc. can be used 
efficiently Also, the physical effects implied by characteristic nuclear scales (the nuclear 
radius Ra ~ A 1 / 3 and the average nucleon-nucleon distance d ~ 2 fm) are best discussed 
in the lab frame. 

In this section we elaborate on several aspects relevant to deep-inelastic scattering as 
viewed in coordinate space. We first discuss the coordinate space resolution of the DIS 
probe. Then we introduce coordinate space distribution functions (so-called Ioffe-time 
distributions) of quarks and gluons and summarize results for free protons. A detailed 



3 A simple interpretation is also possible in the laboratory frame using light-front dynamics. 
In this description, the scattering cross section is determined by the square of the target ground 



state wave function (for a review and references see e.g. [115]). 
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discussion of nuclear effects in coordinate space distributions follows next. In the final 
part we comment on the relationship between lab frame and infinite momentum frame 
pictures. 



4-1 Deep-inelastic scattering in coordinate space 



We follow here essentially the discussion in Ref. ||116| | (see also p] J117|J118|J119| and refer- 
ences therein). Consider the scattering from a free nucleon with momentum P M = (M, 0) 
and invariant mass M in the laboratory frame. The four-momentum transfer = (u, q), 
carried by the exchanged virtual photon, is taken to be in the (longitudinal) ^-direction, 
q = (0_l, (fe) with q 3 = y/u 2 + Q 2 and Q 2 = —q 2 . In the Bjorken limit, v 2 > Q 2 > M 2 with 
x = Q 2 / (2Mv) fixed, the light-cone components of the photon momentum (q ± = v ± q 3 ) 
are q + ~ 2v and q~ ~ —Mx. All information about the response of the target to the 
high-energy virtual photon is in the hadronic tensor 



W^(q,P) ~ / d A y^y(P\J,(y)UO)\P), 



(see Eq.flpp). Using 



1 / , _ _ ,x _ Mx , 

q-y = ^{q y +<? y ) - q±-y±-vy — —y - q±- y±, 



(4.2) 



one obtains the following coordinate-space resolutions along the light-cone distances y ± 

t±y 3 - 



5y ~ — and 5y + ~ . 

y v Mi 



(4.3) 



At y~ = the current correlation function in Eq. ([4.1|) is not analytic since it vanishes 
for y + y~ — (y±) 2 < because of causality (see e.g. 0). Indeed in perturbation theory it 
turns out to be singular at y~ = 0. Assuming that the integrand in ( |4.1| ) is an analytic 
function of y~ elsewhere, this implies that is dominated for q + — > oo by contributions 
from y~ = 0. Causality implies that, in the transverse plane, only contributions from 
(y±) 2 ~ l/Q 2 are relevant: deep-inelastic scattering is dominated by contributions from 
the light cone, i.e. y 2 = 0. 



Furthermore, Eq. (|4.3|) suggests that one probes increasing distances along the light cone 



as x is decreased. Such a behavior is consistent with approximate Bjorken scaling ||117 . 
The coordinate space analysis of nucleon structure functions in Section |4.3| confirms this 
conjecture. In the Bjorken limit the dominant contributions to the hadronic tensor at small 
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Fig. 4.1. Two examples of diagrams illustrating the space-time pattern of deep inelastic scat- 
tering. 



x come from light-like separations of order y + ~ 1/ (Mi) between the electromagnetic 
currents in (14.11). 



In the laboratory frame these considerations imply that deep-inelastic scattering involves 
a longitudinal correlation length 

V3~y = 1 ( 4 - 4 ) 

of the virtual photon. Consequently, large longitudinal distances are important in the 
scattering process at small x. This can also be deduced in the framework of time-ordered 
perturbation theory (see Section |4~5|), where I determines the typical propagation length 



of hadronic configurations present in the interacting photon. 

The space-time pattern of deep-inelastic scattering is illustrated in Fig.[0] in terms of 
the imaginary part of the forward Compton amplitude: the virtual photon interacts with 
partons which propagate a distance y + along the light cone. The characteristic laboratory 
frame correlation length I is one half of that distance. 



4-2 Coordinate- space distribution functions 



Especially when it comes to the discussion of the relevant space-time scales which govern 
nuclear effects in deep-inelastic scattering, it is instructive to look at quark and gluon 
distribution functions in coordinate rather than in momentum space. In this section we 
prepare the facts and return to the underlying dynamics at a later stage. 

It is useful to express coordinate-space distributions in terms of a suitable dimensionless 
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variable. For this purpose let us introduce the light-like vector with n 2 = and 
P -n = P — P 3 . The hadronic tensor receives its dominant contributions from the vicinity 
of the light cone, where y is approximately parallel to n. The dimensionless variable 
z = y ■ P then plays the role of a coordinate conjugate to Bjorken x. It is helpful to bear 
in mind that the value z = 5 corresponds to a light-cone distance y + = 2z/M pa 2 fm in 
the laboratory frame or, equivalently, to a longitudinal distance / = y + /2 pa 1 fm. 

In accordance with the charge conjugation (C) properties of momentum-space quark and 
gluon distributions, one defines coordinate-space distributions by ||120|| : 



Q(z,Q 2 ) = / dx q(x,Q 2 ) + q(x,Q 2 



Q v (z,Q 2 ) = J dx 
o 



Q(z,Q 2 ) = / dx x g(x,Q 2 ) cos(zx), 



sm[z x) 



q(x, Q 2 ) - q(x, Q 2 ) cos(z x) 



(4.5) 
(4.6) 
(4.7) 



where q, q and g are the momentum-space quark, antiquark and gluon distributions, 
respectively. Flavor indices are suppressed here for simplicity. 

At leading twist accuracy, the coordinate-space distributions ( f4.5| - |4.7| ) are related to for- 



ward matrix elements of non-local QCD operators on the light cone ||121| , |122 



Q(z,Q 2 
Qv(z,Q 2 



1 



AiP ■ n 



AP -n 



g(z,Q 2 )=n»n" 



(PMy)tr{y)il>{0)\P) Q >-(y++-y), 

(P$(y) iT{y) V>(0)|P)q 2 + (y <- -y), 



:(P\G, x (y)ny)G\(0)\P) Q i- 



2(P ■ n) 



(4.8) 
(4.9) 
(4.10) 



Here tp denotes the quark field and the gluon field strength tensor. The path-ordered 
exponential 



T(y) = Pexp 



igy»J dXA^Xy) 
o 



(4.11) 



where g denotes the strong coupling constant and N 1 the gluon field, ensures gauge 
invariance of the parton distributions. Note that an expansion of the right-hand side 
of Eqs. (p~5|-pi~7D and (pL~8"H4.10P around y = (and hence z = y ■ P = 0) leads to the 
conventional operator product expansion for parton distributions 0J3,§]. 
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Fig. 4.2. Coordinate-space quark and gluon distributions resulting from the CTEQ4L 
parametrization of momentum-space distributions, taken at a momentum transfer Q 2 = 4 GeV 2 . 
A sum over the u and d quarks is implied in the functions Q v and Q fll6| 1 . 



The functions Q(z), Q v (z) and Q(z) describe the mobility of partons in coordinate space. 
Consider, for example, the valence quark distribution Q v (z). The matrix element in (|4.9|) 
has an obvious physical interpretation: as illustrated in Fig. [4.1| a, it measures the overlap 
between the nucleon ground state and a state in which one quark has been displaced along 
the light cone from to y. A different sequence is shown in Fig. [4.1| b. There the photon 
converts into a beam of partons which propagates along the light cone and interacts with 
partons of the target nucleon, probing primarily its sea quark and gluon content. 



4-3 Coordinate- space distributions of free nucleons 



In this section we discuss the properties of coordinate-space distribution functions of free 
nucleons. Examples of the distributions (" |4.8| — 4. 10 ) using the CTEQ4L parametrization 



123|| of momentum-space quark and gluon distributions taken at a momentum scale Q 2 
4 GeV 2 , are shown in Fig. 



Some general features can be observed: the C-even quark distribution Q(z) rises at small 
values of z, develops a plateau at z > 5, and then exhibits a slow rise at very large z. At 
z < 5, the gluon distribution function zQ{z) behaves similarly as Q(z). For z > 5, zQ(z) 
rises somewhat faster than Q(z). The C-odd (or valence) quark distribution Q v (z) starts 
with a finite value at small z, then begins to fall at z ~ 3 and vanishes at large z. Recall 
that in the laboratory frame, the scale z ~ 5 at which a significant change in the behavior 
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of coordinate-space distributions occurs, represents a longitudinal distance comparable to 
the typical size of a nucleon. 



At z < 5 the coordinate-space distributions are determined by average properties of 
the corresponding momentum-space distribution functions as expressed by their first few 
moments ||124jJL25| . For example, the derivative of the C-even quark distribution Q(z) 



taken at z = equals the fraction of the nucleon light-cone momentum carried by quarks. 
The same is true for the gluon distribution zQ{z) (the momentum fractions carried by 
quarks and by gluons are in fact approximately equal, a well-known experimental fact). 
At z > 10 the coordinate-space distributions are determined by the small-x behavior of 
the corresponding momentum space distributions. Assuming, for example, q(x) ~ x@ for 
x < 0.05 implies Q(z) ~ z~^~ x at z > 10. Similarly, the small- x behavior g(x) ~ x 13 leads 
to zQ{z) ~ z~^~ x at large z. For typical values of (3 as suggested by Regge phenomenology 
one obtains Q v ~ z~ 0,5 while Q(z) and zQ(z) become constant at very large z. 



The fact that Q(z) and zQ(z) extend over large distances has a natural interpretation 
in the laboratory frame. At correlation lengths I much larger than the nucleon size, both 
Q{z) and zQ(z) reflect primarily the partonic structure of the photon which behaves 
like a high-energy beam of gluons and quark-antiquark pairs incident on the nucleon. For 
similar reasons, the valence quark distribution Q v (z) defined in Eq.( [4.9|) has a pronounced 
tail which extends to distances beyond the nucleon radius. An antiquark in the "beam" 
can annihilate with a valence quark of the target nucleon, giving rise to long distance 
contributions in Q v . A detailed and instructive discussion of this frequently ignored feature 
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Fig. 4.3. Contributions from different regions in x to the J~2 combination of coordinate-space 
quark and antiquark distributions at Q 2 = 4 GeV 2 [116]. 
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can be found in Ref.[126|. 



Finally we illustrate the relevance of large distances in deep-inelastic scattering at small 
x. In Fig. [4.3| we show contributions to the nucleon structure function F2 in coordinate 
space, 

^2(z,Q 2 ) = ( — F?{x,Q 2 )sm{zx), (4.12) 
J x 




which result from different windows of Bjorken x. This confirms once more that contri- 
butions from large distances ~ l/(Mx) dominate at small x. 



4-4 Coordinate-space distributions of nuclei 



The implications for scattering from nuclear targets, especially for coherence phenomena, 
are now obvious. If one compares, in the laboratory frame, the longitudinal correlation 
length I from Eq. (f4.4j) with the average nucleon-nucleon distance in the nucleus, d ~ 2 
fm, one can clearly distinguish two separate regions: 

(i) At small distances, I < d, the virtual photon scatters incoherently from the individual 
hadronic constituents of the target nucleus. Possible modifications of the coordinate 
distribution functions ( f4"75| - |4.7| ) in this region are caused by bulk nuclear effects 



such as binding and Fermi motion. 

(ii) At larger distances, I > d, it is likely that several nucleons participate collectively 
in the interaction. Modifications of the coordinate distribution functions are now 
expected to come from the coherent scattering on at least two nucleons in the target. 
Using I ~ l/(2Mx), this region corresponds to x < 0.05. 

This suggests that the nuclear modifications seen in coordinate-space distributions will 
be quite different in the regions / > 2 fm and I < 2 fm. This is best demonstrated by 
studying the ratios of nuclear and nucleon coordinate space distribution functions: 



tz F2 (z,q 2 ) 
n v (z,Q 2 ) 
n g (z,Q 2 ) 



f A ^F 2 A (x,Q 2 )nn(zx) 
Jjf F 2 N (x,Q 2 )sin(^x) 

Q^(z,Q 2 V 

zG A (z,Q 2 ) 



E f ej Qj(z } Q 2 ) 



zG N (z } Q 2 Y 



(4.13) 
(4.14) 
(4.15) 
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Fig. 4.4. The coordinate space ratio Kf 2 at Q 2 = 5 GeV 2 for 4 He, 12 C, and 40 Ca from Ref.El. 



The ratios lZp 2 have been obtained for different nuclei from an analysis of the measured 



momentum space structure functions | 119 |. Furthermore, the ratios of valence quark and 



gluon distributions have been calculated in ||116| as sine and cosine Fourier transforms 



7a - [4.71) °f momentum space distribution functions which result from an analysis of 



nuclear DIS and Drell-Yan data ||127| (see also Section 5.6). 



In Fig. |4.4j we show the ratio TZf 2 for Q 2 = 5 GeV 2 taken from ||119|| . The most prominent 



feature is the pronounced depletion of 1Zf 2 at I > 2 fm caused by nuclear shadowing. At 
/ < 1 fm, nuclear modifications of lZp 2 are small, and deep-inelastic scattering proceeds 
incoherently from the hadronic constituents of the target nucleus. The intrinsic structure 
of individual nucleons is evidently not much affected by nuclear mean fields. In momentum 
space, on the other hand, the pronounced nuclear dependence of the structure function F% 
at x > 0.1 evidently results from a superposition of long and short distance contributions 
as seen in Fig.|4.3|. (For a detailed discussion see Ref. [|119|| .) 



In Fig.^TJ] we show the valence quark and gluon ratios TZ V and TZg for 40 Ca from Ref. ||116| . 



They behave similarly as the structure function ratio 1Zf 2 , where the depletion of gluons 
at large distances is most pronounced. It is interesting to observe that in coordinate space, 
shadowing sets in at approximately the same value of I for all sorts of partons. In momen- 
tum space, shadowing is found to start at different values of x for different distributions 
127|1 . Finally note that the shadowing effect continues to increase for distances larger 



than the nuclear diameter. 
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Fig. 4.5. Coordinate-space ratios at Q 2 = 4 GeV 2 for gluon distributions, valence-quark distri- 
butions, and the F 2 structure function in 40 Ca [ 1 16 1 . 



The results shown in Fig.fl5] clearly emphasize the important role of gluons in the shad- 
owing process. Of course the incident virtual photon does not directly "see" the gluons. In 
the primary step the photon converts into a quark-antiquark pair. At small Bjorken-a;, the 
subsequent QCD evolution of this pair rapidly induces a cascade of gluons. This cascade 
propagates along the light cone over distances which can exceed typical nuclear diameters 
by far: the high energy, high Q 2 photon behaves in part like a gluon beam which scatters 
coherently from the nucleus. This offers interesting new physics. The detailed QCD analy- 
sis of nuclear shadowing can in fact give information on the "cross section" cr ff N for gluons 
incident on nucleons, and a simple eikonal estimate using TZ g at asymptotic distances I 
suggests that this <t 9 n is indeed large, comparable to typical hadronic cross sections (see 



also Refs. 128,125 



In summary, a coordinate space representation which selects contributions from different 
longitudinal distances, lucidly demonstrates that nuclear effects of the structure function 
F 2 and parton distributions are by far dominated by shadowing and have a surprisingly 
simple geometric interpretation. 



4-5 Deep-inelastic scattering in standard perturbation theory 



It is instructive to illustrate the previous results by looking at the lab frame space-time 
pattern of the (virtual) photon-nucleon interaction from the point of view of standard 



46 





target target 



(a) time (b) 

Fig. 4.6. The two possible time orderings for the interaction of a (virtual) photon with a nucleon 
or nuclear target: (a) the photon hits a quark in the target, (b) the photon creates a qq pair 
that subsequently interacts with the target. 



time-ordered perturbation theory. The two basic time orderings are shown in Figs. [4.6| a 
and fOlb: 



(a) the photon hits a quark or antiquark in the target which picks up the large energy 
and momentum transfer; 

(b) the photon converts into a quark-antiquark pair which propagates and subsequently 
interacts with the target. 

For small Bjorken-x the pair production process (b) dominates the scattering amplitude, 
as already mentioned. This can also be easily seen in time-ordered perturbation theory as 
follows (see e.g. |25j and references therein): the amplitudes A a and Ab of processes (a) and 



(b) are roughly proportional to the inverse of their corresponding energy denominators 
AE a and AEf,. For large energy transfers v ^> M one finds: 



AE a = E a (t 2 )-E a (t 1 ) w-^; ' + 
AE b = E h {t 2 )-E b (h) w ,r + Q 



n i/2 | \Pl) + Q 2 

2u 



2v 



(4.16) 
(4.17) 



1/2 . 



where \j) q j is the average quark momentum in a nucleon and /i is the invariant mass 
of the quark-antiquark pair. We then obtain for the ratio of these amplitudes: 



A 



A b 



A E h 



AE n 



Mx 



Pt 



1/2 



(4.18) 
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When analyzing the spectral representation of the scattering amplitude one observes that 
the bulk contribution to process (b) results from those hadronic components in the photon 
wave function which have a squared mass fi 2 ~ Q 2 (see Section [5.4.1 ). The ratio in 



Eq. ([4.18 ) is evidently small for x <C 0.1. Hence pair production, Fig. [4.6| b, is the leading 



lab frame process in the small-x region. On the other hand, at x > 0.1, both mechanisms 
(a) and (b) contribute. 

In process (b) the photon couples to a quark pair which can form a complex (hadronic 
or quark-gluon) intermediate state and then scatters from the target. At small x deep- 
inelastic scattering can therefore be described in the laboratory frame in terms of the 
interaction of quark-gluon components present in the wave function of the virtual photon 



Fig. [4.7|) . The longitudinal propagation length A of a specific photon-induced quark-gluon 



fluctuation with mass \x is given by the inverse of the energy denominator fl4.17|) : 

A ~ J- = 2V ^ J- (4 19) 

AE b {i 2 + Q 2 2xM> { ' 



which coincides with the longitudinal correlation length I of Eq.( |4.4| ). For x < 0.05 the 
propagation length A exceeds the average distance between nucleons in nuclei, A > d ~ 
2fm. For a nuclear target, coherent multiple scattering of quark-gluon fluctuations of the 
photon from several nucleons in the nucleus can then occur, and this is clearly seen in the 
coordinate space analysis discussed in the previous section. 

For larger values of the Bjorken variable, x > 0.2, the propagation length of intermediate 
hadronic states is small, A < d. At the same time the process in Fig. |4.6| a becomes promi- 
nent, i.e. the virtual photon is absorbed directly by a quark or antiquark in the target. 
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Now the incoherent scattering from the hadronic constituents of the nucleus dominates. 



4-6 Nuclear deep-inelastic scattering in the infinite momentum frame 



Let us finally view the deep- inelastic scattering process in the so-called infinite momentum 
frame where the target momentum is large. In this frame the standard parton model 
applies in which a snapshot of the target at the short time scale of the interaction reveals 
an ensemble of almost non-interacting partons, i.e. quarks and gluons. 

Consider the scattering from a nucleus which moves with large longitudinal momentum 
Pa ~ ~~ > 00, where Pn is the average longitudinal momentum of the bound nucleons 
| |130| , |131| , |132| |. The average nucleon-nucleon distance in nuclei is now Lorentz contracted 
as compared to the lab frame: d tn f ~ d Ma/Pa ~ 2fmM/PN. On the other hand the 
derealization of a parton with longitudinal momentum fraction x in the nucleon is given 
according to the Heisenberg uncertainty principle by 51 ~ 1/xPn- At small Bjorken-x, 
x < \/d m 'P^ ~ 0.1, the wave functions of partons from different nucleons have a chance 
to overlap, i.e. d in ^ < 51. Therefore, at x < 0.1 we expect an enhanced interaction 
between partons coming from different nucleons. One can anticipate that, at x <C 0.1, 
the parton derealization extends over the whole nucleus. This is where the quark and 
gluon fluctuations of the photon interact simultaneously with the parton content of several 
nucleons. 



5 Shadowing in unpolarized deep-inelastic scattering 



As outlined in Section 373, the most pronounced nuclear effect in lepton-nucleus DIS 



is shadowing. For small values of the Bjorken variable (x < 0.1), the nuclear structure 
functions are significantly reduced as compared to the free nucleon structure function 
P 2 N . Equivalently, the virtual photon-nucleus cross section is less than A times the one 
for free nucleons, ct 7 .a < ^4cr 7 » N . The analogous behavior is observed for real photons at 
large energies [y > 3GeV). 

This reduction of nuclear cross sections is reminiscent of the features seen in high-energy 
hadron-nucleus collisions. For example, total cross sections for nucleon-nucleus scattering 
behave as cxna ~ ^4 ' 8 o"nn at center-of-mass energies y/s ~ (10 - 25) GeV 2 ||133| . A simple 
geometric picture interprets this effect as the hadron projectile interacting mainly with 
nucleons at the nuclear surface, leading to ctna ~ 

The quantum mechanical description of shadowing in DIS explains this phenomenon by 
the destructive interference of single and multiple scattering amplitudes. Multiple scat- 
tering becomes important as soon as the lab frame coherence length for the hadronic 
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fluctuations of the photon propagator exceeds the average distance between two nucleons 
in the nuclear target. We have seen in our space-time discussion of Section ^ that this is 
precisely what happens in the region x < 0.1 of the Bjorken variable. 

The physics issue of nuclear DIS at small x is therefore, roughly speaking, the optics 
of hadronic or quark-gluon fluctuations of the virtual photon in the nuclear medium. 
Diffractive phenomena play an important role in this context, as we shall demonstrate. 

At extremely small x (i.e. for x < 10 -3 ) in combination with large Q 2 , the measured 
free nucleon structure functions indicate a rapidly growing number of partons (mostly 
gluons). This is the domain of "high density QCD" where individual partons interact 
perturbatively, at large Q 2 , but their number increases so strongly that effective cross 
sections can become large (for references see e.g. | ]134| , |135| . |136| , |137| . p^B| . |139| ). It is of great 
interest to investigate the transition of the observed shadowing phenomena into this new 
domain, accessible by collider experiments, but so far unexplored for nuclear systems. 

In this section we first concentrate on the relationship between diffractive photo- and 
leptoproduction from nucleons and shadowing in high-energy photon- and lepton-nucleus 
interactions. Then we investigate perturbative and non-perturbative QCD aspects of shad- 
owing. After that we summarize existing models which successfully describe data. Finally 
we outline implications of shadowing for nuclear parton distributions. 



5. 1 Diffractive production and nuclear shadowing 



In the shadowing region, diffractive photo- and leptoproduction of high energy hadrons 
gives a substantial contribution to the (virtual) photon-nucleon cross section as dis- 



cussed in Section [2.6.4| . This suggests that the diffractive excitation of hadronic states, 
7*N — > XN, and their coherent interaction with several nucleons inside the target plays 
an important role for shadowing in high energy photon-nucleus scattering, in a similar 
way as for hadron-nucleus collisions. For this effect to be significant, the following two 
conditions have to be met in the laboratory frame: 

(i) The longitudinal propagation length, or coherence length, 

A = why (5 - 1) 



of the diffractively produced hadronic state of invariant mass Mx, Eq. (|4.19 ), must 
exceed the average nucleon-nucleon distance in nuclei: 

A>d~2fm. (5.2) 
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(ii) In addition, the mean free path lx = (po"xn)~ of the diffractively produced system 
in the nuclear medium must be sufficiently short, at least smaller than the nuclear 
radius 

Note that the mean free path of photons in a nucleus with density p amounts to Z 7 ~ 
(pct 7 n) _1 ~ 550 fm, which is much larger than any nuclear scale. Consequently "bare" 
photons do not scatter coherently from several nucleons and therefore do not contribute 
to shadowing. 

Shadowing results from the coherent scattering of a hadronic fluctuation from at least 
two nucleons in the target, i.e. for A > d. Since the longitudinal propagation length A of a 
diffractively produced hadronic state X decreases with its mass Mx, low mass excitations 
with Mx < 1 GeV are relevant for the onset of shadowing. Equation ( |5.2| ) tells again 
that shadowing in deep- inelastic scattering at Q 2 ^> 1 GeV 2 should start at x ~ 0.1, 
in accordance with the observed effect and in close correspondence with the space-time 
picture described in Section f|. 

For real photons diffractive processes at low mass are dominated by the excitation of the 
p- and w-meson. Significant contributions to double scattering and hence to shadowing 
are therefore expected if the photon energy v exceeds about 3 GeV, in line with the 
experimental data. 

Consider now the scattering process in the laboratory frame. Realistic nuclear wave func- 
tions are well established only in this frame (with the exception of recent efforts to con- 
struct relativistic nuclear model wave functions on the light front, see e.g. [140 - 146]). 
Later, in Section |5.5| , we comment on nuclear shadowing as seen in the Breit frame. We 
first neglect effects due to nuclear binding, Fermi-motion and non-nucleonic degrees of 
freedom in nuclei. They are relevant at moderate and large values of the Bjorken variable, 
x > 0.1, as discussed in Section^. 

The (virtual) photon-nucleus cross section can be separated into a piece which accounts 
for the incoherent scattering from individual nucleons, and a correction from the coherent 
interaction with several nucleons: 

<7 7 *a = Z cr 7 . p + (A — Z) a r * n + 5ct 7 »a- (5.3) 



The single scattering part is the incoherent sum of photon-nucleon cross sections, where 
Z is the nuclear charge. The multiple scattering correction can be expanded in contribu- 
tions which account for the scattering from n > 2 nucleons. Expressed in terms of the 
corresponding multiple scattering amplitudes A^) A we have: 

^'A = ^X>4"a. (5-4) 
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d d d 

(a) 



(b) 



d 



Fig. 5.1. Single (a) and double (b) scattering contribution to virtual photon-deuteron scattering. 
The corresponding cross sections are obtained from the imaginary part of the forward scattering 
amplitude indicated by the dashed line. 



where the photon flux ( |2.32| ) is taken in the limit a; < 1. The leading contribution to 
nuclear shadowing comes from double scattering. Its mechanism is best illustrated for a 
deuterium target on which we focus next. 



5.1.1 Shadowing in deuterium 

In this section we review the basic mechanism of shadowing in real and virtual photon- 
deuteron scattering at high energies u, or equivalently, small x. The 7*-deuteron cross 
section can be written as the sum of single and double scattering parts as illustrated in 

FigQ 

<7 7 »d = Cr 7 *p + CT 7 *n + &7 7 * d . (5.5) 



The first two terms describe the incoherent scattering of the (virtual) photon from the 
proton or neutron, while 



accounts for the coherent interaction of the projectile with both nucleons. 

For large energies, v > 3 GeV, or small values of the Bjorken variable, x < 0.1, the 
double scattering amplitude A^) d is dominated by the diffractive excitation of hadronic 
intermediate states (Fig.|5TT| b) described by the amplitude T 7 * n ^ X n- At the high energies 
involved it is a good approximation to neglect the real part of this amplitude. In fact, 
we expect ReT 7 .N^XN ^ 0.15 IhiT 7 *n^xn by analogy with high-energy hadron-hadron 
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scattering amplitudes (see e.g. When including such non-zero real parts, the double 
scattering contribution changes by less than 10% 
dependence for unpolarized scattering 



147||. We neglect the spin and isospin 



70]. Of course, these degrees of freedom play a 



crucial role in polarized scattering as we will discuss in Section 7.3 



Treating the deuteron target in the non-relativistic limit gives ||148|J149| , |i50| , |i51|| : 



A% = -^ fd 3 r\Mr)\ 



d 3 k m ,, s e ik - r 



x 



£ / T ^™ {k) (g - k o y -kl- (q 3 - hf - M x + *e T ™^n(/c),(5.7) 



where k^ = (ko,k) with = (k±,ks) is the four-momentum transfered to the nucleon, 
and ipA is the deuteron wave function normalized as / d 3 r \ipd{ r )\ 2 — 1- The sum is taken 
over all diffractively excited hadronic states with invariant mass Mx and four-momentum 
Px = q — k. We write 



EI^«n| 2 = 64 7 tMV / dM x ^f f (5.1 



d 2 ati¥ 



X Ami X 



in terms of the diffractive production cross section, with t = k 2 . The limits of integration 
define the kinematically permitted range of diffractive excitations, with their invariant 
mass Mx above the two-pion production threshold and limited by the center-of-mass 
energy W = y/s of the scattering process. We introduce the spin-averaged deuteron form 
factor, 

S d (k) = Jd 3 re* k - r \Mr)\ 2 , (5.9) 



perform the integration over the longitudinal momentum transfer in Eq. ( |5.7| ) and then 

l(2) 



(2) 

take the imaginary part of the amplitude A\* d . Actually ks is simply fixed by energy- 



momentum conservation: 



Q 2 + M| 



T) (5-10) 
2v A 



which coincides with the inverse of the longitudinal propagation length ( |4.19| ) of the inter- 
mediate hadronic state. Note that the minimal momentum transfer required to produce 
a hadronic state diffractively from a nucleon at rest amounts to t m i n ~ —k 2 (Mx)- 
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When all steps are carried out, the result for the double scattering correction is [ 148 ,146] 

w 2 



5a rd = --Jd 2 k ± J dMZS d (k ± ,k 3 a\-\Mx)) 

4m 2 - 



dM\dt 



(5.11) 



This equation establishes the close relationship between shadowing in deep-inelastic scat- 
tering and diffractive hadron production. It becomes even more transparent for x <C 0.1, 
i.e. large A. In this limit the magnitude of shadowing is determined just by the ratio of 
diffractive and total 7*N cross sections. To verify this let us parametrize the t-dependence 
of the diffractive production cross section entering in Eq.( [5.11| ) as 



d 2 a% f 
dM\dt 



-B\t\ 



d 2 a% f 



dM\dt 



-Bk 2 



d 2 af/J 



t=o 



dMldt 



(5.12) 



t=o 



neglecting the k 3 dependence of t. Data from FNAL and HERA on diffractive photo- and 
leptoproduction of hadrons with mass M| > 3 GeV 2 give B ~ (5 . . . 7) GeV~ 2 |53]J5^|J55| . 
In the diffractive production of low mass vector mesons (p, u and (f) from nucleons, a 
range of values B ~ (4 . . . 10) GeV -2 has been found, depending on Q 2 and on the incident 
photon energy (for a review and references see e.g. |]49| , |64]| ). Clearly, the soft deuteron form 



factor selects momenta such that the double scattering correction in (|5.11| ) is dominated 
by diffractive production in the direction of the incident photon. 



In Fig. |5.2j we show the deuteron form factor (|5.9|) weighted by the exponential t-dependence 
( |5.12| ) and integrated over transverse momentum, 



B / \ -l\ 



d 2 k A 
W) 



S d (k ± ,k 3 = \-')e 



-Bk 2 



(5.13) 



as obtained with the Paris nucleon-nucleon potential ||152|| for a slope parameter B = 8 
GeV -2 . We observe T d ~ constant as long as the longitudinal propagation length A 
exceeds the deuteron size (r 2 )]/ 2 ~ Afm. From Eq.( |5.2| ) one then finds that hadronic 
states with an invariant mass 



Ml < M 2 max 



W 2 + Q 2 
M(r 2 )l /2 



- Q 2 



(5.14) 



contribute dominantly to double scattering. Combining Eqs. fl5.lID and ( |5.14j ) gives the 
following approximate expression for the shadowing correction in the limit of large longi- 
tudinal propagation length A: 
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Fig. 5.2. Integrated deuteron form factor from Eq.( |5.13| ) for an average slope B 
The dotted curve corresponds to B = 0. 



8 GeV 2 . 



M- 



5a, 



7*d ' 



-STr-F^A" 1 ->0) / 



c 7 di// 

ciM 2 



4m2 



x dM^cit 



> 7 *N 



(5.15) 



t=o 



In the last step we have neglected contributions to the integrated diffractive production 
cross section a i j from hadronic states with invariant masses M max < Mx < W . Since 
da^/J I dM\ drops strongly for large M x as discussed in Section |2.6.4j , this approximation 
is justified at large center of mass energies W or, equivalently, at small x. For the ratio 
between deuteron and free nucleon structure functions we then obtain: 



Rd ~Ff 



2 a. 



7*N 



a 



diff 
(T 7 » N 



(5.16) 



We use o"^Ij//ct 7 »n ~ 0.1 for the fraction of diffractive events in deep-inelastic scattering 
from free nucleons, as suggested by experiment (see Section |2.6.4|) . Furthermore we take 



B = 8 GeV . One finds that shadowing at x 0.1 in deuterium amounts to about 2%, 
i.e. i?d ~ 0.98. The effect is small because of the large average proton-neutron distance in 



the deuteron, but the result agrees well with the experimental data shown in Fig|L5. 



5.1.2 Shadowing for heavy nuclei 

The diffractive production of hadrons from single nucleons also controls shadowing in 
heavier nuclei for which this effect is far more pronounced than in the deuteron. It is an 
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empirical fact that nuclear shadowing increases with the nuclear mass number A of the 
target (see Section |3.3|) . For A > 2 the hadronic state which is produced in the interaction 



of the photon with one of the nucleons in the target may scatter coherently from more 
than two nucleons. However, double scattering still dominates since the probability that 
the propagating hadron interacts with several nucleons along its path decreases with the 
number of scatterers. The double scattering contribution to the total photon-nucleus cross 
section cx 7 *a is obtained by straightforward generalization of the deuteron result ( 5.11|) 



149,150 



+00 +00 

(2) 



<V'A 



8n J d 2 b J dz\ J dz 2 P\\b, Z\\ b, z 2 ) ■ 

d 2 n di JJ 

(5.17) 



-00 Zl 
r2 



W diff 

dM^cos[(z 2 -z 1 )/X] 7 * N 



4ml 



dM\dt 



As illustrated in Fig. |5.3| a diffractive state with invariant mass Mx is produced in the 
interaction of the photon with a nucleon located at position (6, z\) in the target. The 
hadronic excitation propagates at fixed impact parameter b and then interacts with a 
second nucleon at z%. The probability to find two nucleons in the target at the same 
impact parameter is described by the two-body density p£ (b, zi,b, z 2 ) normalized as 
/ d 3 r d 3 r' p$ \r,r') = A 2 . The cos[(z 2 — Zi)/\] factor in Eq. (|5.17|) implies that only 
diffractively excited hadrons with a longitudinal propagation length larger than the aver- 
age nucleon-nucleon distance in the target, A > d ~ 2fm, can contribute significantly to 
double scattering. 

Note that nuclear short-range correlations are relevant only if the coherence length of the 
diffractively excited states is comparable to the range of the short-range repulsive part of 
the nucleon-nucleon force, i.e. for A < 0.5 fm. In this case the shadowing effect is negligible. 
Nuclear correlations are therefore not important in the shadowing domain and the target 
can be considered as an ensemble of independent nucleons with p A (r, r ) « Pa(^)pa(^ ), 
where px is the nuclear one-body density ||153| , |154]| . 



With increasing photon energies or decreasing x down to x <ti 0.1, the longitudinal prop- 
agation length of diffractively excited hadrons rises and eventually reaches nuclear di- 
mensions. Then interactions of the excited hadronic states with several nucleons in the 
target become important. A simple way to account for those is a frequently used equation 



derived by Karmanov and Kondratyuk [155]: 



5<7 7 *a = — b7T 



+00 +00 

J d 2 b J dz\ J dz 2 p A (b, z x ) p A (b, z 2 ) 

— OO Zl 
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Fig. 5.3. Double scattering contribution to deep-inelastic scattering from nuclei. 



U' 2 



dM|cos[(>2 - zi)/A] 



j2 diff 



4m 2 



dM\dt 



exp 



tRiO 



0"XN 



•-2 



p A (M) 



(5.1? 



The exponential attenuation factor describes the elastic re-scattering of the diffractively 
produced hadronic states from the remaining nucleons in the target. The hadron-nucleon 
scattering amplitudes are assumed to be purely imaginary and enter in Eq. fl5.18p through 
the cross sections <txn- 



Equation (|5.18 ) has been applied in several investigations of nuclear shadowing using dif- 
ferent models for the diffractive photoproduction cross section. The more detailed results 
are discussed in Section [5^ , but we can get a simple estimate of nuclear shadowing at 
small Bjorken-x already by just looking at the relative amount of diffraction in DIS from 
free nucleons ||66|| . We restrict ourselves to the double scattering correction ( |5.17 ). For 
x <C 0.1, the coherence length A of the hadronic states which dominate diffractive pro- 
duction in Eq. ( |5.17j ), exceed the diameter of the target nucleus. In the limit A — > oo we 
find: 



a 



(2) 

7*A 



+oo +oo 

-S7T B a^l f J j d 2 b J &z\ J dz 2 p A (b, z x ) p A (b, z 2 ) 



(5.19) 



The slope parameter B and the integrated diffractive production cross section a^JJ have 
been introduced as in Eqs.( |5.12 , 5.15 ). 

For the nuclear densities in Eq. (|5.19|) we use Gaussian, 
/ o \ 3/2 



p A (r) = A 



271^), 



exp 



3r 2 



(5.20) 
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and square-well parametrizations, 



Af(^V) 3/2 for r< ,/f(r 2 )i /2 , 

p A ( r ) = { ^ \5(r 2 )J V3\ /A (521) 

i otherwise, 



with the mean square radius (t 2 )a 
Ra 



a 7 *A/^4cT7*N is easily worked out: 



j ^3 r r 2 p A (j.^ J a. For both cases the shadowing ratio 



l-CA 



B 



diff 



(r 2 ) A J oy N 



(5.22) 



For Gaussian nuclear densities one finds C = 3, while C = 2.7 in the square-well case. 
Using again typical values for the ratio of diffractive and total 7*N cross sections, 



a 



diff 1 

7*N / C"7*N 



0.1, and for the slope parameter, B ~ 8 GeV~ 2 , the magnitude of Ra 
comes out in very reasonable agreement with experimental values as shown in Table [3]. 
This estimate may be simple (in fact, higher order multiple scattering must be included 





6 Li 


12 C 


40 Ca 


131 Xe 


Ra 


0.93 


0.84 


0.73 


0.65 




0.94 ±0.07 


0.87 ±0.10 


0.77 ±0.07 


0.67 ±0.09 



Table 1 

The shadowing ratio Ra estimated according to Eq.( 5.22 ) in comparison to experimental data 



for various nuclei. The data are taken from Ref.[75,76,77] at the smallest kinematically accessible 
values of Bjorken-x (namely, x ~ 10 -4 ). 

in a more detailed analysis) but it certainly confirms that shadowing in nuclear DIS is 
governed by the coherent interaction of diffractively produced states with several nucle- 
ons in the target nucleus. A more detailed investigation of the connection between HERA 
data on diffraction and shadowing effects measured at CERN and FNAL can be found in 



Ref.[|T56 



Inelastic transitions between different hadronic states are neglected in Eq. (|5.18j ). They 
cannot be treated in a model-independent way. Estimates of such higher order diffractive 
dissociation contributions have been performed for high-energy hadron-nucleus scattering 
J133| . |157| | . In the case of neutron-nucleus collisions at center-of-mass energies s ~ 200 GeV 2 
they amount to about 5% of the total reaction cross section. For rising energies the relative 



importance of inelastic transitions is expected to grow [|157 



Given the important role of diffractive production, we can now enter into a more detailed 
discussion of the ^-dependence of shadowing. The coherence lengths A of hadronic states 
with small masses become comparable with nuclear dimensions for x < 0.1. As A increases 
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Fig. 5.4. The quantity log (1 — R\) as a function of logx for data taken on lead \77\. The dashed 
line corresponds to the asymtotic energy dependence ( |5.24| ) with e = 0.1. 



with decreasing x, the shadowing effect grows steadily for x < 0.05. At x <C 0.1 it is also 
the energy dependence of the diffractive production cross section and of the hadron- 
nucleon cross section <txn which influences the x-dependence of shadowing. 



Consider the shadowing ratio Ra = a^x/Aa^ = 1 — <5cr 7 » A/^4oyN, parametrized as: 



Ra-i 



X 



(5.23) 



with a constant c at small Q 2 where data are actually taken, and a characteristic exponent 
e. At asymptotically large energies Regge phenomenology suggests e ~ 0.1 (see Section 

HD D 



In Fig.^| we show the quantity 

log (1 — -Ra) = logc — e logx, (5-24) 



plotted versus logx in comparison with data taken on Pb at small Q 2 . This plot confirms 
that, for x < 3 • 10 -3 , the shadowing effect indeed approaches the high-energy behavior 

4 At the typical average center of mass energies W < 25 GeV used at experiments at CERN and 
FNAL a somewhat stronger energy dependence is expected through the kinematic restriction to 
diffractively produced hadronic states with masses Mx < W. 
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expected from the Regge limit of diffractive production. Deviations from this asymptotic 
behavior at larger values of x indicate how shadowing gradually builds up as the coherence 
length A oc x~ l starts to exceed nuclear length scales for low mass diffractively produced 
states. At sufficiently high energy or small x, the coherence length becomes comparable to 
nuclear dimensions even for heavy hadronic intermediate states. Once a major fraction of 
diffractively produced states contribute to shadowing it starts to approach its asymtotic 
high-energy behavior. Note this asymptotic behavior sets in when the coherence lengths 
Ax of low mass hadronic fluctuations of the photon exceed by far the dimension of the 
nucleus. For example at x — 3 • 10~ 3 and Q 2 ~ 0.7 GeV 2 , which corresponds to the onset 
of the asymptotic behavior in Fig. |5.4j , the p meson coherence length becomes X p ~ 36 fm. 



5.1.3 Shadowing for real photons 

Data on the diffractive production of hadrons in high-energy photon-nucleon interactions 
have been summarized in Section |2.6.4|. They are useful to gain insight into the relative 



importance of p, uo and <f> mesons, as compared to heavier states, for nuclear shadowing 
with real photons. 

Diffractive 7N — > XN production with Mx ^ 1 GeV involves primarily the light vector 
mesons p, u and <fi. Nuclear shadowing at photon energies v up to about 200 GeV is largely 
determined by the coherent multiple scattering of those diffractively produced vector 
mesons. Their propagation lengths A ~ 2v/m\ easily exceed nuclear dimensions as soon as 
v > 20 GeV. With rising energies additional contributions to shadowing from diffractively 
produced states with larger masses, Mx > 1 GeV, become increasingly important. 



This behavior is illustrated for DIS from deuterium in Fig. |5.5| , where we show the ratio of 
the total photon-deuteron cross section compared to the free photon-nucleon cross section, 



Rd = cr 7 d/2<7 7 N, from Ref. ||147|l . The empirical photon-proton cross section from ||158|| has 



been used for o" 7 n- The shadowing correction ( |5.11 ) has been calculated using a fit for the 



diffractive photon-nucleon cross section from Ref . [ 147|[ 



The observed energy dependence of shadowing in Fig. |5.5| results from two sources as 
pointed out previously: the dependence of the diffractive and total photon-nucleon cross 
sections on energy implies Rd — 1 ~ v ' 1 for the shadowing ratio. An additional increase of 
shadowing with rising energy v comes from diffractively produced states with large mass, 
Mx > 1 GeV, which become relevant at high energies. 



5.1.4 Shadowing in DIS at small and moderate Q 2 



So far nuclear shadowing has been measured only in fixed target experiments. The kine- 
matic conditions of such experiments imply that the data for x < 0.01 had to be taken 



at small four-momentum transfers, Q 2 < 1 GeV 2 , as discussed in Section |3.3.1|. The cor 
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Fig. 5.5. The shadowing ratio = cr 7 d/2<7 7 N as a function of the photon energy. The dashed 
line shows the vector meson contribution. The experimental data are taken from the E665 
collaboration J91[| . (The energy values of the data have to be understood as average values which 
correspond to different x-bins.) 

responding energy transfers are typically 50 GeV < v < 300 GeV. The conclusions just 
drawn for real photons apply here too: nuclear shadowing as measured by E665 and NMC 
receives major contributions from the diffractive production and multiple scattering of 
vector mesons. 

In the intermediate range 0.01 < x < 0.1, on the other hand, the E665 and NMC mea- 
surements involve momentum transfers up to Q 2 ~ 30 GeV 2 . At Q 2 > 1 GeV 2 vector 
meson contributions to diffraction and shadowing decrease (Section |2.6.4j) and hadronic 
states with masses M\ ~ Q 2 become relevant. The data reveal that the Q 2 -dependence 
of nuclear shadowing is very weak (Section |3.3.1|) . This suggests that high-mass hadronic 
components of the photon which dominate the measured nuclear shadowing at Q 2 > 1 
GeV 2 , interact strongly with the target, just like ordinary hadrons. The following section 
gives a schematic view of the scales involved, as outlined in Ref. | |159| ]. 

5.2 Sizes, scales and shadowing 

Consider DIS at small x in the lab frame. In this frame of reference the important feature 
is the nuclear interaction of hadronic fluctuations of the virtual photon (see Section §). 
Since the photon and its hadronic configurations carry high energy, the transverse separa- 
tions and longitudinal momenta of their quark and gluon constituents are approximately 
conserved during the scattering process. These hadronic configurations can be classified as 
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"small" or "large" , depending on their transverse extension. "Large" configurations have 
hadronic sizes of order Aqq D ~ 1 fm, whereas "small" configurations are characterized by 
sizes which scale as Q^ 1 . 

The contribution of a given hadronic fluctuation, h, to the photon-nucleon interaction 
cross section is determined by its probability weight u> 7 *h in the photon wave function, 
multiplied by its cross section UhN- The virtual photon-nucleon cross section is: 

c 7 *n = X) W 7*h Cr bN- (5.25) 
h 



The coherent interaction of the virtual photon with several nucleons behaves differently. 
For example, the contribution of a hadronic fluctuation to double scattering, which dom- 
inates shadowing, is proportional to its weight in the photon wave function multiplied 
by the square of its interaction cross section. The double scattering correction to virtual 
photon-nucleus scattering is: 

a 7*A ~ Yl w 7*h (^Iin) 2 - (5.26) 

h 



Now, the probability to find a quark-gluon configuration of large size is suppressed (up to 
possible logarithmic terms) by Aq CD /Q 2 as compared to configurations with small trans- 
verse sizes. On the other hand the interaction cross sections of hadronic fluctuations are 
proportional to their squared transverse radii. These properties and their consequences 
for the cross sections in Eqs.( |5.25l ) and (|5.26| ) are summarized in Table 2. For the scatter- 



fluctuation h 




<7hN 






small size 


1 


i/Q 2 


i/Q 2 


i/Q 4 


large size 


^-qcd/Q 2 




i/Q 2 


V(A 2 1CD Q 2 ) 



Table 2 

Relative contributions of small- and large-size hadronic components of a virtual photon to DIS 
and shadowing at large Q 2 |15£]. 



ing from individual nucleons one finds that both, large- and small-size configurations give 
leading contributions ~ 1/Q 2 to the photon-nucleon cross section ( |5.25|) . On the other 



(2) 

hand contributions from small-size components to the shadowing correction cr^J A are sup- 
pressed by an additional power 1/Q 2 as compared to large-size configurations (apart from 
contributions related to diffractive production from the whole nucleus, not considered in 
this schematic picture). 

In view of these scale considerations, we can now understand some of the previously 
mentioned empirical facts which, on first sight, seemed unrelated: 
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• Nuclear shadowing varies only weakly with Q 2 . 

• The energy dependence of nuclear shadowing for x < 0.01, as measured with fixed 
target experiments at CERN and FNAL, is reminiscent of hadron-nucleus collisions. 

These features follow from the fact that, to leading order in Q 2 , shadowing is primarily 
determined by the interaction of large-size hadronic fluctuations of the exchanged photon, 
even at large Q 2 . These hadronic configurations are expected to interact like ordinary 
hadrons. 



Note, those observations can be applied to diffraction as well as to shadowing, given 
that the two phenomena are closely connected as established in the previous sections: 
diffraction is also a scaling effect, i.e. it survives at large Q 2 . Its energy dependence is 
expected to behave similarly as in hadron collisions. (For limitations to this simple picture 
see Sections O] and IO. ) 



5.3 Nuclear shadowing and parton configurations of the photon 



The results of the previous sections are eludicated by making contact with the underlying 
basic QCD and the parton structure of the virtual photon. The photon wave function can 
be decomposed in a Fock space expansion, 

|7) = c |7o) + Cqq\qq) + c qq - g \qqg) +..., (5.27) 



in terms of a "bare" photon state I70) and partonic (quark-antiquark and gluonic) com- 
ponents. At large Q 2 the minimal Fock component \qq) dominates the hadronic part of 
I7), higher Fock states enter with powers of the strong coupling a s . Let us now have a 
closer look at this minimal Fock component. 

Consider a virtual photon of four- momentum q^ = iy,q), with q = (0, 0, q 3 ) and q 3 > 
defining the longitudinal direction, and Q 2 = —q 2 . Let this photon split into a quark- 
antiquark pair as sketched in Fig. |5.6| . The quark has a four-momentum k^ = (ho, k) with 
k = (k±, kz). The fraction of the photon light-cone momentum carried by the quark is 

k + kn + k-x , 
i= — = ■ 5.28 

q + v + qz 



The momentum fraction of the antiquark with k^ = q^ — k^ is obviously 1 — £. 

For the longitudinally polarized photon, the wave function of its minimal qq fluctuation 
in momentum space is proportional to the longitudinal component of the quark pair 
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quark 



-n antiquark 



Fig. 5.6. Decomposition of a virtual photon into a quark-antiquark pair at large 



current, multiplied by its propagator 
The quantity 



2 rr^+kl 



-1 



where m q is the quark mass |160| . 



M 2 _ m\ + fci 



(5.29) 



can be interpreted as the squared effective mass of the propagating qq pair. 



It is useful to perform the two-dimensional Fourier transform with respect to the transverse 
quark momentum conjugate to the transverse separation b of the qq pair. Neglecting 
the quark mass at large Q 2 and using generically one single quark flavor, the squared wave 
function of the qq component coupled to the longitudinally polarized photon becomes 
||161L|162L|163|: 



^Q 2 e{i-eKi[bQJi{i-i) 
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(5.30) 



where K n denotes modified Bessel functions. The 7* — > qq coupling is proportional to the 
fine structure constant a ~ 1/137. The corresponding result for a virtual photon with 
transverse polarization is: 



3a 
2^2 



Q A m - c + a - o 2 )Kf[bQ mi - . 



(5.31) 



Consider now DIS from a free nucleon at small Bjorken-x and large Q 2 . At the high 
energies involved the photon in the laboratory frame acts like a beam of qq pairs, and one 
can write the cross section for the longitudinally or transversely polarized virtual photon 
with the nucleon in the form ||162| , |164| 



L.T 



1 

d 2 b Jdi \^f(b,o 



qqN 



(b,x), 



(5.32) 
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Fig. 5.7. Short distance interaction of a color singlet quark-antiquark pair with a nucleon through 
two-gluon exchange. 



using the wave functions of the leading qq fluctuations. These wave functions as well as 
the gg-nucleon cross section a qq ^ depend on the transverse separation b of the quark pair. 
Since the modified Bessel functions in Eqs. (|5.30| , |5.31|) drop as Ki t o(y) — > e~ y at large 
y, the wave functions i/j q q receive their dominant contributions from configurations with 
transverse size 



Q 2 £(i-0' 



(5.33) 



Consequently, qq configurations at large Q 2 with comparable momenta of the quark and 
antiquark, £ ~ 1 — £ ~ 1/2, have small transverse size, b 2 ~ ~ ^/Q 2 -, or equivalently, 
large transverse momentum. The interaction of these "non-aligned" configurations with 
the nucleon is therefore determined by the short transverse distance behavior of the cross 
section <r g qN which can be calculated using perturbative QCD. The reasoning goes as 
follows. At large Q 2 the leading mechanism responsible for the short distance interaction of 



the qq pair with the nucleon is two-gluon exchange (see Fig. |5.7|) . The (color singlet) qq pair 
acts as a color dipole. Its interaction strength with the nucleon or any other (color singlet) 
hadron is determined by the squared color dipole moment, hence a qq ^ is proportional to 
b 2 for small transverse separations b. In the leading-logarithmic approximation valid at 



large Q one derives [ 164,165 



&qqTSl{b, x) 



7T 



ct s (Q 2 )b 2 xg N (x,Q 2 



(5.34) 



with the strong coupling constant a s . The Q 2 scale in ( |5.34|) is set by Q 



2 ~ 1/b 2 . All 

non-perturbative effects are incorporated in the gluon distribution g^(x, Q 2 ) of the target 
nucleon. 
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Small qq configurations interact only weakly according to Eq.( |5.34 ). This is the case for 
the kinematic conditions realized in fixed target experiments at CERN and FNAL (see 



Sections |2.6.1| and |3.3.1|) . The situation is different at x <C 0.1 and Q 2 ^> lGeV 2 , the 



extreme region accessible at HERA. Here the strong increase of the nucleon structure 
function at very small x is accompanied by a correspondingly strong increase of the 
gluon distribution function. The gluon density becomes so large with decreasing x that, 
even for small b 2 ~ l/Q 2 , the cross section cr g ^ N can eventually reach magnitudes typical 
for ordinary hadrons |fL6l 



If either the quark or the antiquark becomes soft (that is: if either the momentum fraction 
£ or 1 — £ tends to zero), large qq separations contribute to the wave functions ( |5.30| , |5.31|) . 
In this limit perturbative QCD is not applicable. The interaction cross section for such 
large-size configurations with small transverse momentum is supposed to be similar to 



typical hadron-nucleon cross sections [167 



A detailed analysis of the "transverse" wave function Q5.31 ) shows that both "small" 
(non-aligned) and "large" quark- antiquark configurations give leading l/Q 2 contributions 
to the transverse photon-nucleon cross section in accordance with our previous discussion. 
The situation is different for longitudinally polarized photons (see e.g. (|5.30Q ). In this case 
the contributions from "soft" quarks (with £^0orl— £^0) are suppressed as l/Q 4 
so that, to leading order in the strong coupling a s , only small size qq pairs contribute 



to er^, N . At next to leading order in a s , the Fock expansion ( 5.27|) introduces quark- 



ant iquark-gluon states. Large size qqg configurations are now important, and they are not 
suppressed by additional powers of l/Q 2 ||16£ 



At small momentum transfers, Q 2 < 1 GeV 2 , configurations of large size dominate the qq 
wave function. Strong interactions between quark and antiquark now lead to the formation 
of soft hadronic fluctuations including vector mesons and multi-pion states. Consequently, 
photon-nucleon cross sections at small x and small Q 2 receive important contributions 
from the low mass vector mesons. For example, at Q 2 ~ 0.5 GeV 2 almost half of the 
measured nucleon structure function F 2 N comes from p, uj and <fi mesons according to the 



calculation in Ref . [|169| . 



So far we have focused this discussion on free nucleons. Similar considerations apply to 
deep-inelastic scattering from nuclei which involves the interaction of hadronic compo- 
nents of the virtual photon with the nuclear many-body system. To leading order in a s 
the photon- nucleus cross sections are now obtained from Eq.( |5.32| ) replacing o^n by the 
corresponding gg-nucleus cross section ct^a- The cross section a^A of any hadronic fluc- 
tuation h interacting with a nucleus at high energies, can be related to the cross section 
for the scattering from free nucleons by the Glauber-Gribov multiple scattering formalism 
|. For a Gaussian nuclear density Q5.20| ) this leads to: 



o- hA « A a hN 



1 - ChN 



3 A-l 

16n (t 2 )a 



exp 



3 A 2 



+ . 



(5.35) 
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where A is the propagation length associated with the hadronic fluctuation. Double scat- 
tering gives a negative correction proportional to the squared cross section of the hadronic 
fluctuation. Only those hadronic configurations with large interaction cross sections con- 
tribute significantly to shadowing. Furthermore, since the nuclear mean square radius 
behaves approximately as (t 2 )a ~ A 2 / 3 , the magnitude of the double scattering correc- 
tion grows for large nuclei with the radius of the target, i.e. proportional to A 1 / 3 . The 
exponential in ( |5.35| ) ensures that only hadronic fluctuations with propagation lengths A 



larger than the target dimension contribute significantly to shadowing. For small-sized 
fluctuations, interesting effects beyond those covered by Eq. ( |5.35| ) arise from diffractive 
production on the whole nucleus. 

In accordance with our discussion in the previous section we can conclude: 

i) In the fixed target experiments at CERN and FNAL, where small values of x < 0.01 are 
accessible only at small average momentum transfers, Q 2 < 1 GeV 2 , nuclear shadowing 
is governed by interactions of configurations with large transverse sizes. Contributions 
from the vector mesons p,u and turn out to be particularly important. 

ii) At very small x together with very large Q 2 , the growing number of partons in the 
photon-nucleus system makes them interact like ordinary hadrons, even if the parton 
configurations have small transverse sizes inversely proportional to Q 2 . One now expects 
a complex interplay between soft (large-size) and hard (small-size) partonic components 
of the interacting photon which can no longer be classified by simple book-keeping in 
powers of l/Q 2 - 



54 Models 



In the following we sketch several models which have been used quite successfully to 
describe nuclear shadowing as measured in experiments at CERN and FNAL. As before 
we restrict ourselves to lab-frame descriptions. We do not aim for completeness but rather 
emphasize common features of various models and their implications for the underlying 
mechanism of nuclear shadowing. 



5.4- 1 Vector mesons and aligned jets 

As discussed in Section |5.3| , the quark- ant iquark fluctuation of a virtual photon starts out 
with a transverse size b 2 ~ [Q 2 C,(1 — £)] _1 where £ is the fraction of longitudinal photon 
momentum carried by one of the quarks. "Non- Aligned" qq configurations with £ ~ 1/2 
have small transverse size and interact weakly; "aligned" ones with £ ~ or £ ~ 1 have 
large transverse size and are likely, by subsequent strong interactions, to turn into vector 
mesons if the qq invariant mass matches appropriately. 
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Models which combine aspects of vector meson dominance and the aligned-jet picture 
167| | are described in Refs. ||169| , |170| | . Their starting point is the hadronic spectrum of the 



virtual photon exchanged in the deep-inelastic scattering process. The spectral function, 
II(s), is determined by the cross section for electron-positron annihilation into hadrons, 
where s = q 2 is the squared photon or e + e~ center-of-mass energy: 

12 TT Z CJe+e-.+n+n-is) 



with 



n(g 2 ) = ~ E(2vr) 3 5 4 (g - Px) (0| J„(0) |X) (X| J"(0) |0) . (5.37) 



Here J M is the electromagnetic current operator. The sum in Eq.( |5.37| ) is taken over all 
hadronic fluctuations of the photon with four-momenta Px = Q and squared invariant 
masses /z 2 = p\ = q 2 ■ At small center-of-mass energies, s < lGeV 2 , the spectrum ( |5.37| ) 
is dominated by the vector mesons p, u and as shown in Fig. |5.8| . The high energy 
spectrum at s > 1 GeV 2 is characterized by quark- ant iquark continuum plateaus together 
with isolated charmonium and upsilon resonances. 

The lab frame space-time pattern of deep-inelastic scattering (Section [|) suggests that 
the nucleon structure function at small x can be described by the following expression 

if(*,Q 2 ) = ^ / V ,V n 22 [ ^a m (W,n 2 ;0- (5-38) 

^44 



The basic idea behind this ansatz is the following. For x -C 1, or large lab frame 
propagation length A ~ 2u/{Q 2 + jj 2 ) of a given qq fluctuation of mass /x, the vac- 
uum spectrum H(fi 2 ) remains more or less unaffected by the presence of the target 
nucleon. The high-energy virtual photon with v ^> Q 2 /2M behaves like a beam of 
hadrons with masses p < /i max . Their maximum possible mass is determined by the 
condition that A must exceed the size R ~ 5M _1 of the target nucleon, so that (roughly) 
Pmax ~ V vM. The interaction of this beam with the nucleon is described by the cross 
section <ThN which depends on /i 2 and on the hadron/photon-nucleon center-of-mass en- 
ergy W = \flMv + M 2 — Q 2 ~ \JlMv = Q/ \fx. For a qq pair treated to leading order 
in a s , it also depends on the fraction £ of the photon light-cone momentum carried by the 
quark. The sum in Eq.( |5.38p is taken over hadronic fluctuations of the photon with fixed 
invariant mass. The ansatz neglects contributions to the forward virtual photon scattering 
amplitude in which the mass \i can change during the interaction. 
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Fig. 5.8. Cross section a e + e -^ hadrons /cr e+e -^ +/ ,- . 



The small- a; structure function Ff(x,Q 2 ) as given in Eq. (|5.38 ) is governed by contri- 
butions from intermediate hadronic states with an invariant mass p 2 ~ Q 2 . For small 
momentum transfers, Q 2 < lGeV 2 , low mass vector mesons p, u and are of major 
importance. Their dominance leads to the scale breaking behavior F 2 N (x,<5 2 ) ~ Q 2 for 
Q 2 — > at small x. 
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Fig. 5.9. The nucleoli structure function F^ for small x plotted against Q 2 . The full line has 
been obtained in Ref . 1 169 1 from Eq.( |5.38| ), The dashed line indicates the contribution of the 
vector mesons p, lo and 0. The data are from the NMC [171]. 



For larger momentum transfers, Q 2 > m 2 « 1 GeV 2 , the structure function F 2 N is de- 
termined primarily by the interactions of quark-antiquark pairs from the qq continuum. 
The color singlet nature of hadronic fluctuations of the virtual photon implies that their 
interaction cross section is proportional to their transverse size. Quark pairs with mo- 
menta "aligned" along the direction of the virtual photon have a large transverse size. 
Their cross sections should be comparable to typical hadronic cross sections. On the other 
hand "non-aligned" quarks are characterized by small transverse size. Their cross sections 
should therefore be small. 



With these ingredients, Eq. fl5.38Q gives a good description of the free nucleon structure 
function Ff for x < 0.1 and moderate Q 2 . A comparison with data from NMC is shown 
in Fig. |5.9| . While the vector meson contribution vanishes as 1/Q 2 for large Q 2 , the qq 
continuum pairs are responsible for scaling, F^^x, Q 2 ) ~ ln(Q 2 ), at large Q 2 . Note however 
the importance of vector mesons at small Q 2 . One finds that at Q 2 = 1 GeV 2 almost half 
of F2 at x = 0.01 is due to vector mesons. At Q 2 = 10 GeV 2 they still contribute about 
15%. In these calculations the vector meson part of the spectrum IT(/z 2 ) is p5f| : 



n( v V) 




(5.39) 



with V = p, oj, 0, the empirical vector meson masses my and the coupling constants g p = 
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Fig. 5.10. Results from Ref. |169| 1 for the shadowing in He, Li, C, and Ca compared to experimental 
data from NMC (dots and squares) and FNAL (t riangles) |77jj . The dashed curves show 

the shadowing caused by the vector mesons p, oj and <p only, the solid curves are the results 
including the qq continuum. 



5-0, g u = 17.0 and = 12.9. The vector meson-nucleon cross sections are <t p n 
mb, ct^n ~ 10 mb. 



(T^N ~ 25 



Nuclear structure functions F% for x < 0.1 are expressed in an analogous way as in 
Eq.( |5.38|) , with the hadron-nucleon cross sections cr hN replaced by the corresponding 
hadron- nucleus cross sections OhA- The relation between <7hA and <JhN is given by Glauber- 
Gribov multiple scattering theory, see Eq. (|5.35| ). In Fig. |5.10| we present typical results for 
the shadowing ratiof^] Ra = F^/F^ from Ref. j |169|l . 



Finally we comment on the observed weak Q 2 -dependence of the shadowing effect. In 



5 It is common practice to normalize such that it represents the nuclear structure function 
per nucleon. 
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Fig. 5.11. The slope b = ^(i^ 11 / 'F^) / 'd\aQ 2 indicating the Q 2 dependence of the shadowing 



ratio Sn/C. The calculation is described in [169]. Data are taken from 



the spectral ansatz ( |5.38| ) the given value of Q 2 selects that part of the hadron mass 
spectrum around /z 2 ~ Q 2 which dominates the interaction, and hence determines which 
cross sections crhN(/-t 2 ) contribute significantly to the multiple scattering series. While the 
interaction cross sections decrease as 1/fJ 2 with increasing mass as required by Bjorken 
scaling, pairs which are aligned with the photon momentum interact with large cross 
sections, even for large /i, and therefore produce strong shadowing. This is the reason for 
the very weak overall Q 2 -dependence of shadowing in this framework. A comparison of 
results from Ref.[|T69| with NMC data for the slope b of the ratio -Ff n /F 2 c « a + b\nQ 2 
is presented in Fig. |5.11| . For a more detailed discussion of these issues including QCD 
corrections, see Ref.|l[]. 



5.4-2 Vector meson dominance and pomeron exchange 



As indicated in Eqs. ( [5.1 1],|5. 18| ) , nuclear shadowing is directly related to the diffractive pro- 



duction cross section da^f/dM^ dt or, equivalently, to the diffractive structure function 

Diffractive production at Q 2 < 1 GeV 2 is dominated by the excitation of the vector 
mesons p, uj and <fi. Their contributions can be described within the framework of vector 
meson dominance (see e.g. p5|). Neglecting transitions between different vector mesons 
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Fig. 5.12. Shadowing in Xe. Details of the calculation are given in Ref.[154|. The dashed curve 
shows the contribution of vector mesons p, oj and (f>, while the solid curve includes pomeron 
exchange. The data are from the E665 collaboration 



and omitting contributions from longitudinally polarized virtual photons one finds: 



da diff(Y) 



dM\dt 



a 



n( v )(MD 



Ml 



{Q 2 + Miy 



a 



XN 



(5.40) 



Here the vector meson part fl5.39|) of the photon spectral function enters. Combining 
Eqs. (|5.18| , [5.40| ) shows that the contribution of vector mesons to nuclear shadowing van- 
ishes indeed as l/Q 2 . 



The diffractive excitation of heavy mass states is commonly parametrized according to the 
Regge ansatz as in Eq.( p)7D . Most descript ions concentrate on the dominant contribution 
from pomeron exchange. The pomeron structure function is modeled in agreement 
with available data on diffraction. At large Q 2 it is supposed to scale, i.e. it depends 
at most logarithmically on Q 2 . This leads to scaling for nuclear shadowing at large Q 2 . 

< lGeV 2 



On the other hand, at 



one assumes 



Q 2 [172 



which ensures that 
vanishes at Q 2 = 0, 



the shadowing correction to the nuclear structure function, 5F^ , 
just like F£ itself. Investigations of shadowing effects along theses lines can be found in 
]^g72}|T^g7g , g7g , [T7§ . In Fig.|o~T2l a typical result for shadowing in Xe from |Tol 
shown. 



is 
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5.4-3 Generalized vector meson dominance 



Generalized vector meson dominance models describe deep-inelastic lepton scattering at 
small x purely in terms of hadronic degrees of freedom [|1 77| . |1 78| . |l 79| . |1 80| . |T8T| . |1 82| . |i~83l . |l 8 4|l85|] 
For a free nucleon this leads to the following picture: prior to its interaction with the target 
the virtual photon fluctuates into a hadronic state with invariant mass p. This fluctuation 
scatters from the target and converts into a hadronic state with mass p! . For transversely 
polarized photons this translates into a forward Compton amplitude of the form: 



A T 



(W,Q 2 



d p 2 dp 



p(p,p';W) 



(/i 2 + g 2 )(^ + g 



(5.41) 



with a double spectral distribution p depending on the photon-nucleon center-of-mass 
energy W. Here the integrals over initial and final hadronic fluctuations and their propa- 
gators are made explicit. 

The continuum of hadronic intermediate states which determines the double spectral 
function p(p, p!\ W) is commonly approximated by a discrete set of narrow vector meson 
states V n {n = 1,2,...). The resulting transverse photon-nucleon cross section is: 



M 2 M 



a 



7*N 



£ 9m Ml lq^ mn{W) Ml + Q 2 g n ■ (5 ' 42) 



In Refs.| |18U| , |184j| the vector mesons are assumed to be equally spaced in mass, starting 



with the p-meson. The photon- vector meson couplings g n are chosen to reproduce average 
scaling in e + e~ annihilation into hadrons (see Fig. [5.8|) . S mn denotes the imaginary part 
of the vector meson-nucleon transition amplitude, V m N — > V n N, in the forward direction. 
For diagonal terms it is equal to the total V n -nucleon cross section, S n „ = cr(V n N), which 
is taken to be constant. 

The next step in simplification is to consider only diagonal (m = n) and nearest off- 
diagonal (m = n ± 1) contributions. A fine-tuned cancelation between the corresponding 
amplitudes S mn leads to a reasonable description of the nucleon structure function at 
moderate momentum transfers Q 2 . 

An extension of this approach to nuclear targets involves multiple scattering of hadronic 
fluctuations from several nucleons. The multiple scattering process is described by a cou- 
pled channel optical model | |179|jl84| which accounts for the shadowing criteria in Eq.( p.2| ), 
i.e. only those hadronic fluctuations with longitudinal interaction lengths larger than their 
mean free path in the nuclear medium contribute significantly to multiple scattering and 
thus to shadowing. 



GVMD calculations applied to nuclear DIS data can be found in Refs. [ |186| , |187| . 
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5.4-4 Vector mesons and quark scattering 



We add a few remarks and references about approaches dealing with DIS in terms of 
quark dynamics. 



The starting point in Ref.|188| is a description of DIS from nucleons at large Q and small 



x in terms of quark-nucleon scattering | 189 |. The quark- nucleon scattering amplitude is 
formulated using Regge phenomenology and constrained by the quark distributions of 
free nucleons. The interaction strength of quark-nucleon scattering is determined by the 
quark-nucleon cross section, taken to be about 1/3 of the nucleon-nucleon cross section. 
At center of mass energies s ~ 200 GeV 2 one finds a q ^ ~ 13 mb ||190 . 



An extension to DIS from nuclei at small x involves the quark-nucleus scattering am- 
plitude. Its connection with the amplitude for the scattering from free nucleons is given 
through the Glauber-Gribov multiple scattering series. 



In Ref. ||190|| the interactions of strongly correlated quark-antiquark pairs, i.e. vector mesons, 
have been added. One finds that vector mesons carry more than half of the shadowing 
effect measured at E665 and NMC. On the other hand, the interaction of uncorrelated 
quarks is also important to ensure a weak Q 2 -dependence of shadowing. 



5-4-5 Green function methods 



The previously mentioned models have outlined in different ways the ingredients needed 
in order to understand the physics of shadowing: the mass spectrum of quark-gluon fluctu- 
ations of the virtual photon, and the dynamics of the expanding and strongly interacting 
quark-gluon configurations in the surrounding nuclear system. In most of the models 
the longitudinal propagation of hadronic fluctuations of the photon is treated by multiple 
scattering theory, while the transverse degrees of freedom are more or less "frozen" during 
the passage through the nucleus. Several questions are faced in this context. The trans- 
verse size of quark-gluon fluctuations needs to be connected with their effective mass; 
the relationship with diffractive production must be elucidated; higher order terms in the 
multiple scattering series must be systematically incorporated, at least for heavy nuclei. 



A coordinate space Green function method which permits to unify all those aspects has 



been developed in Refs. ||191| , |192|| . This work considers only quark-antiquark fluctuations of 
the photon. It turns out that some previous approximations can be recovered as limiting 



cases. We follow Ref. [ |191|| and give here a brief summary of the essentials. 



Consider the scattering of a virtual photon with high energy v and large squared four- 
momentum, Q 2 > 1 GeV 2 , through a nucleus as illustrated in Fig. |5.13| . The longitudinal 
(z-) direction is defined by the photon three-momentum, as usual. At point z\ the photon 
produces a quark-antiquark pair with transverse separation b±. Along its passage to point 
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Fig. 5.13. Propagation of a quark-antiquark fluctuation of the virtual photon 7* between points 
z\ and Z2 where the pair has transverse separation b\ and b 2 . The Green function G(b2, Z2; &i, z±) 
sums all possible paths of the pair through the nucleus. 

z 2 where it has a transverse separation b 2 , the qq fluctuation experiences multiple inter- 
actions with nucleons in the nuclear target. We are interested in the full Green function 
G(b2, Z2', bi, z\) which describes the propagation of the qq pair from z\ to Z2, including its 
dynamics in the transverse space coordinate. 

This Green function enters in the shadowing part of the 7*-nucleus cross section, as follows. 
One writes: 

5a rA (x, Q 2 ) = cr rA (x, Q 2 ) - A a rN (x, Q 2 ) 



The quantity W(z 2 ,Zi) has dimension (length) 4 and describes the production of a qq 
fluctuation in the process 7*N — > ggN, its propagation from z\ to z 2 , and its subsequent 
conversion back to a virtual photon: 



oc 



OC' 




(5.43) 



1 



W(z 2 ,z 1 ) = Re J d% J d 2 b 2 J dZF f (b 2 ,Z)G(b 2 ,z 2 ;b 1 ,z 1 )F(b 1 ,S) e 








It involves the amplitude 



Hb,0 = l^Y^(b,t)v(b) 



(5.45) 
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for the 7*N — > ggN process leading to a qq pair of transverse separation b in which the 
quark carries the fraction £ of the photon light cone momentum (see also Section |5.3| ). 
The color dipole cross section cr(b) has the characteristic color screening behavior, i.e. it 
vanishes as b 2 at b — > (see also Eq. ( |5.34|) ), and the distribution of transverse separations 



is determined by the amplitude ■0 7 *->gg- The normalization of is such that its Fourier 
transform gives the 7*N — > qqN diffractive dissociation amplitude, 



f(k ± )= Jd 2 bf(b,Oe 



ik±b 



(5.46) 



in plane wave impulse approximation. The phase factor e 1 ^ 2 involves the characteristic 
wave number of the qq fluctuation: 



K 



Q2^ (1 _ 0+m 2 

2^(1-0 



(5.47) 



where m q is the (constituent) quark mass. For £ = 1/2 the resulting k 



= \- 



is just the inverse coherence length of a quark and antiquark which travel side by side. 
(For arbitrary £ this coherence length includes the transverse momentum, A -1 = k + 
kj[2u 



Let us now return to the propagation function G. It satisfies a wave equation |191|| which 
can be made plausible by the following considerations. The longitudinal motion along the 
z-axis is equivalent to the time evolution of the qq fluctuation, represented by the operator 
%-§^-- The transverse dynamics has a kinetic term, 



Vl 2 



2^(1-0 



(5.48) 



with the 2-dimensional Laplacian acting on the transverse separation coordinate b 2 , and 
the denominator reflecting the effective mass of the pair. Interactions of the qq pair with 
the nuclear medium at an impact parameter b are introduced by an absorptive term, 



v(b 2 ,b) 



v(b 2 )p A (b,z 2 ) 



(5.49) 



The wave equation for G is then of the generic form idG/dz = {tkin + v) G or, more 
precisely, 



dz 2 2z/£(l-0 



+ 2°"( & 2)Pa(&, z 2 ) 



G(b 2 , z 2 ]b^ zx) = 



(5.50) 
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with the initial condition G(b 2 , Z\\ bi, z\) = 5 2 (b2 — t>i). One can now discuss several 
interesting limits: 

i) The "frozen" limit 



Take the energy v — > oo so that the kinetic term (|5.48|) vanishes (with some extra care 
required at the kinematic corners, £ = and £ = 1). Then 



G(b 2 , z 2 \ b u zi) = S 2 (b 2 - bi) exp 



a(b 2 



22 

J dz p A (b h z) 



(5.51) 



Inserting this expression into Eqs.( |5.44] , |5.43|) with A — > oo one finds: 



a rA = 2 I d b I d b q q J d£ \ip r ^ qq (b qq , f)| \ 1 - exp 

o 



a(b, 



mi 



dz p A (b,z) 



(5.52) 



and recovers the shadowing correction 5ct 7 *a as in Glauber-Gribov multiple scattering the- 
ory by expanding the exponential. Note the difference compared to the standard Glauber 
eikonal approximation where the cross section (j{b q q) is averaged in the exponent. 

ii) No absorption 



Take the limit a — > in the wave equation ( |5.50|) . Then G reduces to the free Green 
function of the qq pair, 



G(b 2 ,z 2 ;bt,Zi) = -j- / d 2 k± exp 



ik±-(b 2 - &i) + 



ik\(z 2 - gi) 
2^(1-0 



(5.53) 



Inserting this into Eq. ( |5.44|) and using the diffractive dissociation amplitude ( |5.46|) one 
finds 



W( Zl ,z 2 ) = 7^fd£ I d 2 k ± \f(k ± )\ z cos 
o 



r Q 2 za-o+m 2 q + ki 

— [z 2 - Zi) 



2^(1-0 



(5.54) 



We identify the squared effective mass, M\ = (m 2 + fc^_)/£(l — ^), of the qq pair as in 
Eq.( |5.29D and introduce its coherence length A = 2i>/(Q 2 + M|). Inserting Eq.( [5.54j ) into 
Eq.( [5.43| ) one then recovers the double scattering result, Eq.( |5.17| ), with the factorized 

(2) 

two-body density p A (b, z\\ b, z 2 ) = p A (b, z\) p A (b, z 2 ). It is now also apparent how the 
additional absorption factor in Eq. (|5.18|) is obtained, introducing an average cross section 
<7xn in the exponent. 



78 



iii) Propagation in uniform nuclear matter 



Assume that the qq pair moves in a nuclear medium of uniform density px(b,z) = Po = 
const. (p = 0.17 fm~ 3 for normal nuclear matter). Suppose that the color dipole cross 
section is approximated by 



a(b q q) = cb 2 q g 



(5.55) 



with a constant parameter c. In this case the wave equation ( |5.50| ) reduces to 
d 



+ 



v 6 2 2 



IC 



dz 2 2*^(1-0 2 



+ 7T PO b 2 



G(b 2 ,z 2 ;bx,zi) = 0. 



(5.56) 



This is formally reminiscent of the Schrodinger equation for a harmonic oscillator with 



complex frequency. One finds ||191 



G(b 2 ,z 2 ;bx,0) 



a a 
ex P \ - « 



2n smh(ujz 



bf + bt,) coth(a;z) 



with 



2 fcpoY ■ cp 
to = = %- 



2bi ■ b 2 

sinh(ujz) 



(5.57) 



(5.5J 



This is a convenient approximation to account for multiple scattering and absorption of 
the qq fluctuation, still keeping track of its transverse dynamics during its passage through 
the nuclear medium. Instructive results are discussed in Ref. ||191 . 



5.4-6 Meson exchange and shadowing 

Up to now we have concentrated on diffractive contributions to nuclear shadowing, in 
which the nucleons interacting with the virtual photon are left unchanged. The coherent 
interaction of the photon with several nucleons in the target nucleus can also involve 
non-diffractive processes, in particular, reactions in which nucleons change their charge. 
These are commonly described by the exchange of mesons and sub-leading Reggeons. 

Modifications to nuclear structure functions at small x through meson exchange have been 
investigated in Refs.[ 1931 , 194] for deuterium. In this work significant effects come from the 
interaction of the virtual photon with pions emitted from the target proton or neutron. 
Here, as in diffraction, a hadronic state X is produced which subsequently re-scatters from 
the second nucleon. Contributions from the exchange of other mesons, e.g. p and u, turn 
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out to be negligible. For the double scattering contribution through pion exchange one 
finds in analogy with Eq.( |5.11| ): 



w 



Km = \ I d " k ± \ dM x SUk) ^g- (5.59) 



4mJ 



dM\dt ' 



Here d 2 oZ f *^j 'dM\dt is the cross section for the semi-inclusive production of a hadronic 
state with invariant mass Mx from a proton or neutron via pion exchange. The form 
factor in Eq. (|5.59|) accounts for the spin- dependent response of the deuteron: 



S l ( fc ) = \ E/ d 3 P^f(P) <T P k<Tn- kW(P - fc), (5.60) 



where k is the pion momentum and k = k/\k\. The momentum-space wave function of 
the deuteron with polarization m is denoted by ip™. Furthermore the non-relativistic form 



of the pion-nucleon coupling is used [ |195|| . Note that the energy of the exchanged pion is 



determined by fco = — a/ M 2 + P 2 — y M 2 — (P — A;) 2 , where P is the momentum of the 
parent nucleon. We denote the pion four-momentum by fc = (ko, k). For the longitudinal 
pion momentum one has A; 3 rs yM with the pion light-cone momentum fraction y = 
k ■ q/P ■ q, and we introduce t = k 2 along with the usual Bjorken-x. 

It is common to factorize the semi-inclusive differential cross section: 
d 2 a w * 

j^ {x,Q 2 ;y,t) = U/ N {y,t) a r7r (x/y, Q 2 ). (5.61) 

Here the photon-pion cross section is related to the structure function of the pion by 
F^x.Q 2 ) = (Q 2 /4:7i 2 a) a~,* n (x, Q 2 ), and the pion distribution function in the nucleon is 
given by: 

f („ +\ 3 ^NN I-^ttNnWI 2 (~t) fKKO\ 

U/N{y,t) = ——^ T- ^ — y, 5.62 

167T^ (t — m^Y 

with the pion-nucleon coupling constant ^nn and the 7rNN form factor Jvnn normalized 
to unity for on-mass-shell pions, i.e. .F^nn = 1 for t = k 2 = m 2 ,. 

For practical calculations the pion structure function has been approximated by that of 



the free pion, as parametrized in [|196j|197l1 in accordance with Drell-Yan leptoproduction 
data. Here, however, only the region x > 0.1 has been measured. An extraction of the pion 
structure function at small x from semi-exclusive reactions at HERA has been discussed 
recently in gggggg . 
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The resulting pionic correction 5a^* d to double scattering turns out to be positive, i.e. 
it causes "anti-shadowing" . The relative weight of <5cr^, d decreases with decreasing x. At 
typical values Q 2 = 4GeV 2 and 0.001 < x < 0.1 it amounts to around 30% of the overall 
shadowing correction | 193 |. In Ref. ||194|| the pion correction 5cr^* d has been found to be 
negligible at large Q 2 > 10 GeV 2 . Note however that the quoted results depend sensitively 
on the yet unknown pion structure function at small x, the deuteron wave function and 
the choice for the pion-nucleon form factor. 



5.4-7 Discussion 



The models sketched above give quite reasonable descriptions of the data on nuclear shad- 
owing measured at CERN and FNAL. All of them support the general observation that 
nuclear shadowing as measured by NMC ff^^f^!^,!!,!!! and E665 |76|,^|7|,|T1 at 
small x < 0.01 receives major contributions from the coherent interaction of the vector 
mesons p, u and (p. In fact those experiments are performed at small average momentum 
transfers Q < 1 GeV 2 . On the other hand, the observed weak Q 2 -dependence of the 
shadowing effect originates from the coherent interaction of strongly interacting quark- 
antiquark fluctuations with large masses, Mx > 1 GeV. 



5. 5 Interpretation of nuclear shadowing in the infinite momentum frame 



In this section we briefly discuss how nuclear shadowing develops in the infinite momentum 
frame where the parton model for deep-inelastic scattering can be applied. We found in 



Section ^LQ that, in this frame, the wave functions of partons from different nucleons in 
the nucleus start to overlap for x < 0.1. One then expects that the interaction of partons 
belonging to different nucleons increases. Shadowing at small x < 0.1 is supposed to be 
due to the fusion or recombination of partons from different nucleons, thereby effectively 
reducing the quark distributions of each nucleon. At the same time parton fusion leads to 



an enhancement of partons at x > 0.1. In Ref. ||201|l modifications of parton distributions 
due to parton fusion have been derived and found to be proportional to l/Q 2 - Therefore 
parton fusion processes seem to be suppressed at large momentum transfers but can be 
significant at low Q 2 . 

Procedures for modeling nuclear parton distributions at small x have been proposed in 
Refs. ||132|J201| , |202|| . Recombination effects modify these distributions dominantly at a low 
momentum scale Qq where parton fusion is calculated and incorporated in the initial 
quark and gluon distribution functions. Parton distributions at Q 2 > are then derived 
through the calculation of radiative QCD corrections using DGLAP evolution (see Section 



2.4). 



To describe the measured shadowing of the NMC and E665 collaboration a typical scale 
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Qq ~ 0.8 GeV 2 has been used in Refs.[ 132 , 202 |. As a result the empirical shadowing for 
F- can be described. It should be mentioned that the calculation of the recombination 
effect within perturbation theory is certainly questionable at a low momentum scale Qq. 
The results are strongly sensitive to model parameters, such as the initial scale Ql and 
the input parton distributions. 

Note that the recombination effects discussed here involve parton distributions at a low 
momentum scale. This "initial-state recombination" is different from the "radiative recom- 
bination" , discussed in Section |S.3|, which modifies the parton evolution by recombination 
of radiatively produced partons. 



5. 6 Nuclear parton distributions at small x 



Any quantitative QCD analysis of high energy processes involving nuclei requires a de- 
tailed knowledge of nuclear parton distributions. In this section we outline the empirical 
information on their difference with respect to quark and gluon distribution functions of 
free nucleons. 

Let us first focus on the nuclear gluon distribution. The Q 2 -dependence of deep-inelastic 
structure functions at small x is dominated by gluon radiation. One can therefore extract 
nuclear gluon distribution functions from a precise analysis of scaling violations of the 
structure functions _F 2 A . In leading order perturbation theory and in the limit i < 0.1 the 
DGLAP equations ( |f39| , |2T10| ) reduce to the simple form pTJ3| : 



dF 2 (x,Q 2 ) _ as_^ 2 <o /v2 
d In Q 2 3n . 



e 2 f xg(2x,Q 2 ). (5.63) 



This relation, with further inclusion of small corrections from quark contributions, has 
been used in an analysis ||204j| of high statistics NMC data on the Q 2 -dependence of the 



structure function ratio F 2 n / F 2 shown in Fig. |3.4| . The result for the corresponding ratio of 



nuclear gluon distributions, gsn/ gc, is shown in Fig. |5.14| . At x < 0.1 the gluon distribution 
is shadowed, i.e. gsn/gc < 1, in a similar way as the structure function F 2 . This observation 
is quite natural since F 2 at x < 0.1 is dominated by contributions from sea quarks. The 
intimate relation between sea quarks and gluons through DGLAP evolution then also 
suggests shadowing for gluons. 

At 0.05 < x < 0.15 an approximate 8% enhancement of nuclear gluons has been found. 
This observation is in agreement with an analysis of NMC data for inelastic J/ip-\epto- 
production ||110|| as indicated in Fig J5.14 The enhancement of nuclear gluon distributions 
around x ~ 0.1 is consistent with the fact that the total momentum of hadrons is given by 
the sum of the momenta of its parton constituents ||108|,|127||. The empirical information 
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Fig. 5.14. Results from Ref.|204] for the ratio of the Sn and carbon gluon densities, gsn(x)/gc( x ), 
together with the measured ratio of structure functions n (x)/F2(x) [82|- The box represents 
the extraction of the ratio of gluon distributions from J/ip electroproduction data [110]. 



on this sum rule applied to quarks has been presented in Section |3.4| . It implies that the 
momentum carried by gluons is, within error bars, equal in nucleons and nuclei, i.e. 



dxx g^(x, Q 2 



dx x qa(x, Q z 



(5.64) 



Consequently, shadowing of nuclear gluon distributions at small x has to be compensated 
by an enhancement at larger values of x. Assuming the latter to be located in the region 
0.05 < x < 0.15 leads to results similar to the ones shown in Fig. f5.14| [|108|| . 

Note that the close relation between shadowing and diffraction allows to estimate gluon 
shadowing using data on diffractive charm and dijet production from free nucleons. A 
corresponding analysis of HERA data has been carried out in Refs. [|128| , |129|| . It suggests 
significantly larger shadowing for gluons than for quarks. 



Nuclear effects in valence and sea quark distributions can be further disentangled using 



Drell-Yan dilepton production data ||108| , |127|| . The E772 collaboration at FNAL has found 
shadowing for nuclear antiquark distributions at x < 0.1 but no enhancement as discussed 
Combining this with the fact that the nuclear structure function ratio 
for 0.05 < x < 0.2, one concludes that nuclear valence quarks have to be 
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Fig. 5.15. Momentum-space ratios from Ref.[127] relative to the free nucleoli, for gluon distri- 
butions, valence quark distributions, and the F% structure function in 40 Ca at Q 2 = 4 GeV 2 . 



enhanced around x ~ 0.1. From the baryon number sum rule 



Jdxq*(x,Q 2 ) = Jdxq^(x,Q 2 ) 



(5.65) 



one then concludes that nuclear valence quark distributions, must be shadowed at 



x < 0.05. Typical results from Ref. [|108|1 are shown in Fig. |5.15 



To summarize, present data on nuclear shadowing imply that all parton distributions 
are shadowed at x <C 0.1, while only valence quarks and gluons are enhanced around 
x ~ 0.1. The kinematic range where enhancement takes place is related to processes 
which involve typical longitudinal distances of 1 fm in the laboratory. This is the region 
where components of the nuclear wave function with overlapping parton distributions 
should be relevant. In Ref. ||108j| it was suggested that at such distances inter-nucleon 
forces are a result of quark and gluon exchange leading to the observed enhancement. In 
such a picture the enhancement of gluons should increase with the density of the nuclear 
target. A 10% enhancement of glue at x ~ 0.1 in Sn as compared to C would then imply 
a 20% increase of the gluon density in Pb as compared to free nucleons 



TUHI. This would 



imply a dramatic change of gluon fields in nuclear matter at distances of 1 fm between 
nucleons. 



More detailed information on nuclear parton distributions is certainly needed. The shad- 
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owing region - where nuclear effects are large - is of particular interest. Further constraints 
on gluon shadowing from deep-inelastic scattering require data on the Q 2 -dependence of 
nuclear structure functions at smaller values of x as indicated by Eq. (|5.63|) . A more 
quantitative separation of nuclear effects in valence and sea quark distributions could 
be obtained from Drell-Yan dilepton production or neutrino scattering experiments with 
high statistics. On the other hand, an extraction of nuclear parton distributions in hadron 
production processes from nuclei, e.g. lepto- or hadropro duct ion of charmonium or open 
charm (see e.g. |205| , |206| ), is complicated by possible final state interactions and higher 
twist corrections. 



6 Nuclear structure functions at large Bjorken-x 



Deep-inelastic scattering from nuclei probes the nuclear parton distributions. On the other 
hand conventional nuclear physics works well with the concept that nuclei are composed 
of interacting hadronic constituents, primarily nucleons and pions. For x > 0.2 DIS probes 
longitudinal distances smaller than lfm (see Section |]), less than the size of individual 
hadrons in nuclei. In this kinematic region, incoherent scattering from hadronic con- 
stituents of the target nucleus dominates. Such processes explore the quark distributions 
of nucleons bound in the nucleus. 

To gain first insights suppose that the nucleus is described by nucleons moving in a mean 
field. The quark substructure of bound nucleons may differ in several respects from the 
quark distributions of free nucleons. First, there is a purely kinematical effect due to the 
momentum distribution and binding energy of the bound nucleons. This effect rescales the 
energy and momentum of the partonic constituents. To illustrate this recall that for a free 
nucleon the light-cone momentum fraction of partons cannot exceed x = 1. A nucleon 
bound in a nucleus carries a non-vanishing momentum which adds to the momenta of 
individual partons in that nucleon. As a consequence light-cone momentum fractions up 
to x = A are possible in principle, although the extreme situation in which a single 
parton carries all of the nuclear momentum will of course be very highly improbable. On 
the other hand intrinsic properties of bound nucleons, e.g. their size, could also change in 
the nuclear environment. This may lead to additional, dynamical modifications of their 
partonic structure. 



6.1 Impulse approximation 



Nuclei are, in many respects, dilute systems. For example, in elastic proton-nucleus scat- 
tering the proton mean free path is of the order of 5 — 10 fm |[207|| , large compared to the 
average distance between nucleons in nuclei. This observation has motivated the impulse 
approximation which reduces the nuclear scattering process to incoherent scatterings from 
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the individual nucleons (for reviews and references see [p].^[4|fl|2| J208|| ). Final state inter- 
actions of the scattered hadron with the residual nuclear system are neglected at high 
energy. (One should note, however, that the validity of this approximation, illustrated in 
FigJOl, is still under debate as discussed in |Jg Jl^ , gO| , g09],gg , gTl|j 2g and references 
therein.) 

Given the small average momenta of the weakly bound nucleons, their quark sub-structure 
is described by structure functions similar to those of free nucleons [ |19C1| , |213|| . For a nucleon 
with momentum p these structure functions depend on the scaling variable x = Q 2 /2p ■ q 
and on the squared momentum transfer Q 2 . However since the energies and momenta of 
bound nucleons do not satisfy the energy-momentum relation of free nucleons, additional 
freedom arises. 



This becomes immediately obvious from the following simple kinematic consideration. In 
the laboratory frame deep-inelastic scattering from a nucleon bound in a nucleus involves 
the removal energy, — e n , of the struck nucleon: 



e n = M A - M 



(n) 
A-l 



M. 



(6.1) 



Here Ma and M A "^ X denote the invariant masses of the initial nuclear ground state and of 
the nuclear system, with a nucleon-hole state characterized by its quantum numbers n. 
The energy of the interacting nucleon is then: 



Po 



M 



(6.2) 



where Tr = p 2 /2M^} l is the recoil energy of the residual nuclear system. We finally 




(A-l)" 



Fig. 6.1. Impulse approximation for deep- inelastic scattering from nuclei at large Q z 
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obtain for the squared four-momentum of this interacting nucleon: 
p 2 = p 2 - p 2 ~ M 2 + 2M(e n - T R - T) ^ M 2 , 



(6.3) 



with T = p 2 /2M. The squared four- momentum of the active nucleon is obviously not 
restricted by its free invariant mass. It is determined by the nuclear wave function which 
describes the momentum distributions of bound target nucleons as well as the mass spec- 
trum of the residual nuclear system. Consequently, the structure function of a bound 
nucleon can in general depend also on p 2 , not just on x and Q 2 . 



6.2 Corrections from binding and Fermi motion 



In the impulse approximation deep- inelastic scattering from a nucleus at large Q 2 proceeds 
in two steps as shown in Fig. |6.1[ the exchanged virtual photon scatters from a quark with 
momentum k. This quark belongs to a nucleon with momentum p which is removed from 
the target nucleus. Treating the nucleus as a non-relativistic bound state, the nuclear 
structure functions factorize into two terms [ 190 1 : the information about the quark and 



gluon substructure of the nucleons is included in the bound nucleon structure functions 
F^(x/y,p 2 ) and F|^(x/|/,p 2 ). They depend on the fraction x/y = k-q/p-q « k + /p + of the 
light-cone momentum^] of the interacting nucleon carried by the quark, and reduce to the 
corresponding free nucleon structure functions at p 2 = M 2 . Details about nuclear structure 
are incorporated in the distribution function of nucleons with squared four-momentum p 2 
and a fraction y of the nuclear light-cone momentum: 



2 „/2 



P )■ 



(6.4) 



Here 



S(p) = 2tt Y, HPo -M-e n + T R ) \^ n (p) 



(6.5) 



is the spectral function of a nucleon in the nucleus. It is determined by the momen- 
tum space amplitude ^ n (p) — ((A — 1)„, — p|^(0)|A), with ^(0) representing the non- 
relativistic nucleon field operator at the origin r = 0. The squared amplitude |\l/ n (p)| 2 
describes the probability of finding a nucleon with momentum p in the nuclear ground 
state |A), and the remaining A — 1 nucleons in a state n with total momentum —p. In 
Eq. ( |6.5|) the sum over a complete set of states with A — 1 nucleons is taken. Note that the 
spectral function is normalized to A, the total number of nucleons in the nucleus. This 

6 Here the photon momentum is chosen as = (qo, 0j_, (73) with 53 < 0. 
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leads to the proper normalization of the nucleon distribution function in Eq. (|6.4j) (see e.g. 
Refs.@,gggT§). 



The nuclear structure functions are then obtained by a convolution over the squared four- 
momentum of the interacting nucleons and their light-cone momentum fraction. For the 



structure function F% per nucleon this gives ||190 



AF 2 A (x) = Jdy Jdp 2 D N/A (y,p 2 )F^(x/y,p 2 ) } (6.6) 



where we have suppressed the dependence on Q for convenience. 



In the following we examine the convolution integral flS.6| ) in more detail. The momentum 
distribution of nucleons in the nuclear target, 

/f^(p)=£i*n(p)| 2 , (6.7) 



falls rapidly with increasing \p\. This implies that the nucleon light-cone distribution ( 6.4]) 
is strongly peaked around y « 1 andp 2 ~ M 2 , with a typical width Ay ~ Pf/M controlled 
by the Fermi momentum pp. Expanding the bound nucleon structure function in Eq. (|6.6|) 
in a Taylor series around y = 1 and p 2 = M 2 , and integrating term by term, leads to the 
following expression for the nuclear structure function per nucleon [fL90|J214 |, valid in the 
range 0.2 < x < 0.7: 



+2 M#(« , (6 . 8) 



M V dp 



p 2 =M 2 



where F% (x) and _F 2 N "(a;) are derivatives of the structure function with respect to x. The 
mean value of the single particle energy e = po — M is 



and 
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represents the mean kinetic energy of bound nucleons. Except for light nuclei the recoil 
energy Tr in Eq.( |6.2[ ) can be neglected. Then (e) coincides with the separation energy. 
Corrections to Eq. (|6.8| ) are of higher order in (e)/M and (T)/M. Note that the approx- 
imate result for in Eq. (|6.8|) is well justified in the region 0.2 < x < 0.7. Here the 
kinematic condition x/y < 1 in Eq.( |6.6| ) can be ignored in accordance with the underly- 
ing expansion. 



Let us briefly discuss the physical meaning of the different terms in Eq.( |ti.8j ) and their 
implications. The second term on the right hand side of Eq. (|6.8|) involves the average 
separation energy of nucleons from the target. As such it is determined by nuclear binding. 
In the range 0.2 < x < 0.7 it leads to a depletion of the nuclear structure function 
compared to the structure function of a free nucleon. The third term accounts for the 
Fermi motion of bound nucleons and yields a strong rise of the structure function ratio 
F<t /F2 a t large x. Finally, the fourth term in ( |6.8|) reflects the dependence of the structure 
function of a bound nucleon on its squared four-momentum. Note that this contribution 
enters at the same order as binding and Fermi-motion corrections. Information about the 
p 2 -dependence of bound nucleon structure functions is rare. Nevertheless such effects may 
lead to significant modifications of the EMC ratio F^/F^ at moderate and large values of 
x. This has been shown for example in the framework of a simple quark-diquark picture 



for the nucleon 190 . 



An important and not yet completely solved problem with respect to the binding and 
Fermi-motion corrections in Eq.( |6.8|) is a reliable calculation of (e) and (T). In a simple 
nuclear shell model the separation energy is averaged over all occupied levels. One finds 
typical values (e) « -(20 - 25) MeV and (T) « 18-20 MeV. Correlations between 
nucleons change the simple mean field picture substantially and lead to high momentum 
components with \p\ > pp in the nuclear spectral function (|6.5| ). This in turn causes an 
increase of the average separation energy (e) [ pl6| , |217| |. In order to see this let us examine 
the Koltun sum rule ||218| 



(e) + (T) = -2n B , (6.11) 



where /ie is the binding energy per nucleon. This sum rule is exact if only two-body 
forces are present in the nuclear Hamiltonian. With fixed //£ ~ 8 MeV, this sum rule 
tells that an increase of (T) due to high momentum components is accompanied by an 
increase of We refer in this context to a calculation ||219|| of the spectral function 
of nuclear matter based on a variational method. This calculation shows that there is a 
significant probability to find nucleons with high momentum and large separation energies. 
An integration of the spectral function of Ref. [219] gives (T) « 38 MeV and (e) ~ 



-70 MeV [ 217|| . In order to estimate these quantities for heavy nuclei one usually assumes 



216| that the high momentum components of the nucleon momentum distribution are 
about the same as in nuclear matter. Together with Eq.( |6.1lD this leads to (s) ~ —50 MeV. 
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Fig. 6.2. The ratio of nuclear and nucleon structure functions, F^/F^, for iron and gold 
taken from Refs. pTa , 190]. (a) solid curve: calculation in Ref.[215], dotted curve: calculation 
in Ref. |2l4 1. (b) results from Ref.[19C]: without p 2 -dependence of the bound nucleon structure 
function (dashed), and including this p 2 -dependence as obtained from a simple quark-diquark 
picture (full). 



In Fig.^72| we show typical results from Refs. [|190| , pi5|| for iron and gold. We observe that a 
qualitative understanding of the EMC effect can indeed be reached, but at x > 0.5 a more 
quantitative description is still lacking. One should note, of course, that the presentation 
of nuclear effects in terms of the ratio F^/F^ 1 magnifies such effects in a misleading 
manner because F2 itself is small in this region (see also the discussion in Section 
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The impulse approximation picture of nuclear deep-inelastic scattering can also be main- 
tained in a relativistically covariant way [ |212j ]. Here, however, a simple factorization of 
nuclear structure functions into nuclear and nucleon parts as in Eq.( |6.6|) is not possible 
any longer. A relativistic calculation of nuclear structure functions requires relativistic 
nuclear wave functions as well as a more detailed knowledge about the structure of bound 
nucleons. Nevertheless, relativistic effects seem to be small: in an explicit model calcula- 
tion of the deuteron structure function F£ relativistic corrections to the non-relativistic 



convolution ( |6.6D are less than 2% for x < 0.9 [ 220|1 . In the region x > 1, where nuclear 



structure functions are very small however, larger deviations are expected. 

In this context a word of caution is in order. A description of nuclear structure functions 
based on nucleons alone is necessarily incomplete since it violates the momentum sum 



rule M. Non-nucleonic degrees of freedom are briefly discussed in Section |Q. 



6.3 Beyond the impulse approximation 



The quality of the impulse approximation has frequently been questioned (see e.g. 
Ref s. jT|,P|, |l 90| . p 1 0| . |2 1 1| , |2 1 2fl and references therein). Here we give a brief summary of possi- 
ble shortcomings in terms of models for nuclear deep-inelastic scattering which go beyond 
this approximation. 



6.3.1 Quark exchange in nuclei 

The impulse approximation includes only incoherent scattering processes from hadronic 
constituents of the target nucleus. On the other hand, contributions involving several 
bound nucleons could also be important, and their role needs to be examined. 

One such possibility, namely quark exchange between different nucleons, has been investi- 
gated in Refs. |210]J221|| . The nuclear quark wave function which is probed in deep-inelastic 
scattering must be antisymmetric with respect to permutations of quarks. This however is 
not realized in the impulse approximation ( |6.6| ). Antisymmetrization introduces additional 
quark exchange terms between different nucleons in the target. 

Under several simplifying assumptions a softening of the nuclear quark momentum dis- 
tribution due to quark exchange has been found. For small nuclei the effect turned out 
to be significant. For 4 He approximately 30% of the observed depletion of the structure 
function ratio F^ e /F2 at x ~ 0.6 has been associated with quark exchange. Only minor 



modifications have been found for heavier nuclei 221 



While the estimates based on a simple quark exchange model may not be reliable at 
a quantitative level, they certainly point to the fact that the impulse approximation is 
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incomplete as soon as correlations between quarks in several nucleons come into play. 



6.3.2 Final state interactions in a mean field approach 

One of the basic assumptions of the impulse approximation is that interactions of the 
struck, highly excited nucleon with the residual nuclear system can be ignored. In general 
there is no solid basis for this assumption since the debris of the struck nucleon includes 
also low momentum fragments as seen from the target rest frame. A proper treatment 
of their final state interaction requires however a description of the nucleus in terms of 
quark and gluon degrees of freedom. 

Investigations in this direction have been made starting out from a quark model for nuclear 
matter, with nucleons modeled as non-overlapping MIT bags |[211U222|,[223|| . The nucleons 
interact via the exchange of scalar and vector mesons which couple directly to quarks. 
Within the mean field approximation for the meson fields it is possible to describe several 
basic properties of nuclear matter, such as its compressibility and the binding energy per 
nucleon at saturation density. 

This model has been applied to deep-inelastic scattering from finite nuclei using a local 
density approximation [|211|| . The debris of the struck nucleon is represented by a pair 
of spectator quarks bound in a diquark bag. Its interaction with the remaining nuclear 
system in the final state leads to a non-negligible effect on nuclear structure functions: 
while the full calculation including final state interactions allows to reproduce the structure 
function ratio F^/F^ 1 , the impulse approximation overestimates nuclear effects at x ~ 0.6 
by about a factor two. In the framework of this model such a difference is expected since 
the binding of the nucleus is the result of the attractive scalar mean field experienced by 
all three constituent quarks of the interacting nucleon. When ignoring the binding of the 
spectator quark pair, as done in the impulse approximation, one assumes that the active 
quark which takes part in the deep-inelastic scattering process behaves as if it carries the 
binding of all three quarks, a feature which obviously needs to be corrected. 

The mean field approach to nuclear quark distributions is based on several simplifying as- 
sumptions, but it nevertheless points to the possible importance of final state interactions 
and, in more general terms, to the relevance of quark degrees of freedom in nuclei. 



6.4 Modifications of bound nucleon structure functions 



The intrinsic properties of nucleons bound in nuclei can be modified as compared to free 
nucleons. We summarize below two examples of models which deal with such possible 
changes in bound nucleon structure functions. 
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6-4-1 Rescaling 



At intermediate values of the Bjorken variable, 0.2 < x < 0.7, the modification of nuclear 
structure functions F 2 A as compared to the free nucleon structure function F 2 can be 
described by a shift in the momentum scale which enters the structure functions. We 
briefly outline here the basic arguments ||224| , p25| , |226| , p271J228|| . Consider the moments^ 



M n A (Q 2 



dx x r ' 



2 F 2 A ( 



with n even. 



(6.12) 



Assume now that the moments of nuclear and nucleon structure functions are related by 
a shift of their momentum scale: 



(6.13) 



At an arbitrary momentum transfer Q 2 the perturbative QCD evolution equations to 
leading order (see Section |2.4j) give: 



m a {Q 2 ) = M n N (UQ 2 )Q 



(6.14) 



with the rescaling parameter 



UQ 2 ) = ^ 



Qs(M A ) 




(6.15) 



Of course Eq. (|6.14 ) can always be satisfied if one allows different £a for different moments 
n. However, when comparing with data it has turned out that the rescaling parameter 
is independent of n to a good approximation. Consequently, the scale change for the 
moments (|6.14j) can be translated directly into a scale change for the structure functions 
themselves: 



F 2 A (x,Q 2 ) = F?(x,UQ 2 )Q' 



(6.16) 



Good agreement with experimental data can be achieved at intermediate x. For example, 
the EMC structure function data on iron suggest £ Fe ~ 2 for Q 2 = 20 GeV 2 ||226| . 



Rescaling gives a reasonable one-parameter description of nuclear structure functions F^ 
at intermediate x, but it does not offer insights into the physical origin of the observed 



7 For simplicity we use the non-singlet part only. 
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change of scale. One possible suggestion to explain the scale change is a modification of 
the effective confinement scale for quarks in the nucleus as compared to free nucleons 



224,225,226 



Scale changes are not simply related to possible "swelling" of nucleons inside nuclei which 
is constrained by inclusive electron-nucleus scattering data in the quasielastic region. The 
experimentally observed y-scaling indicates a rather small increase of the charge radius 
for bound nucleons. For example, the study of ||229| , [230|| comes to the conclusion that any 



increase of the nuclear radius in nuclei should be less than 6% of its free radius. 

Related discussions and a comparison with nuclear DIS data can be found in Refs. [202,216, 
231-234]. 



6.4-2 Color screening in bound nucleons 



The scenario of Refs. [1,235,23(5] assumes that the dominant contribution to the struc- 



ture function if at large x ~ 0.6 is given by small size (pointlike) parton configurations 
in the nucleon. In a nuclear environment such configurations interact only weakly with 
other nucleons due to the screening of their color charge. It is argued that the proba- 
bility for pointlike configurations is reduced in bound nucleons. In fact, the probability 
to find parton configurations of average size in the nucleon should actually be enhanced 
in nuclei since they experience the attraction of the nuclear mean field. Then the varia- 
tional principle with normalization condition on the wave function implies that small-size 
configurations should indeed be suppressed. 



An estimate of such deformations in the wave function of nucleons bound in heavy nuclei 
gives for x ~ 0.5 ||235|| : 



if 



1 + 



4 (£Q 
E* 



0.7-0.8. 



(6.17) 



Here (U) is the average potential energy per nucleon, (U) ~ —40 MeV, and E* ~ 0.5 
GeV is the typical energy scale for excitations of the nucleon. Since (U) scales with the 
nuclear density, the nuclear dependence of the structure function ratio (|6.17|) is roughly 
consistent with data. 

It should be mentioned that the proposed suppression of rather rare pointlike config- 
urations in bound nucleons does not necessarily imply a substantial change of average 
properties of a bound nucleon, such as its electromagnetic radius Q]. 



94 



6. 5 Pion contributions to nuclear structure functions 



In conventional nuclear physics meson exchange is responsible for the binding of nucleons 
in the nucleus. Therefore deep-inelastic scattering from mesons present in the nuclear 
wave function should give additional contributions to nuclear structure functions. Pions, 
which are responsible for most of the intermediate- and long-range nucleon-nucleon force, 
are supposed to play the prominent role (see e.g. ||237| , |238| , ^39|J24CI| , |241|| ). 



The framework is the Sullivan process ||242|| . Its contribution to the nucleon structure 
function F 9 N reads: 



<TF 2 N (x)= dyU /N (y)FZ(x/y), 



(6.18) 



where 



U/n{v) 



167T 2 



dt- 



-t\J~nNN(t)\' 



2\2 



(6.19) 



specifies the distribution of pions with light-cone momentum fraction y in the nucleon, 
while F£ is the pion structure function. Equation ( 6.18Q describes deep- inelastic scattering 
from a pion emitted from its nucleon source. The nucleon receives a momentum transfer 
equal to the pion momentum k^ = (u,k). The minimal squared momentum transfer 
t = k 2 required for pion emission is t m i n = —M 2 y 2 /{\ — y). One finds that S^F™ gets 
its dominant contributions from pions with momenta \k\ ~ 300 — 400 MeV. Pions with 
smaller momenta are suppressed by the explicit factor y in Eq.( p.l9| ), while pions with 
large momenta are suppressed by the pion propagator and the 7rNN form factor Jvnn 
237 . 



The convolution ansatz in Eq. fl6.18Q suffers from similar problems as convolution for nu- 



clear structure functions discussed in Section |6.2| : the interacting pion is not on its mass 
shell, i.e. k 2 « — k 2 ^ m 2 . Therefore the pion structure function depends also on k 2 . Fur- 
thermore, final state interactions of the pion debris with the recoil nucleon are neglected. 



The detailed treatment of pionic effects in nuclei includes the pion propagation in the 
medium with A resonance excitation, Pauli effects and short range spin-isospin corre- 
lations. All these effects are incorporated in the pion-nuclear response function R(k,u) 
which determines the spectrum of pionic excitations in the nuclear medium. The resulting 
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distribution function of pions in a nucleus is [ 237 

oo \k\-My 



U/a{v) 



167T 



-y 



d\k\ 



M 2 y 2 



, fc 2 |^ NN (fc 2 )| 2 

(c — m£)^ 



(6.20) 



k 2 . Using the the Sullivan description ( 6.18|) the contribution of excess 



where t = uo 2 

pions to the nuclear structure function can be calculated according to: 



dy [U/a(v) - U/M) F^x/y) 



(6.21) 



In the original work in Ref . ||23T| , using the empirical pion structure function, a significant 
enhancement of the ratio F^/F^ 1 was found at x < 0.3. This observation was in agree- 
Later data on nuclear structure functions showed 
0.15 (see Section |3.3.1| ). In addition Drell-Yan 



x 



ment with the early EMC data \74 
only a minor enhancement around 
data from E772 [|107|| have demonstrated that the antiquark distribution in nuclei is not 
significantly enhanced as compared to free nucleons, in disagreement with the first pion 
model calculations. However, as already emphasized in Ref. ||237|| , the pion contribution 
to nuclear structure functions is very sensitive to the strength of repulsive short range 
spin-isospin correlations in nuclei. Variations of this strength by 15% can easily lead to 
30% changes in S^F^- While it is not difficult to accommodate the very small observed 
pionic enhancement within such uncertainties, it is still a challange to arrive at a con- 
sistant overall picture of nuclear DIS which rigorously satisfies the requirements of the 
momentum sum rule. 



6. 6 Further notes 



Related studies of pion field effects as well as other nuclear medium corrections and their 
implications on nuclear DIS have been performed in Refs. ||243| , p44| , p45|J246|| . These studies 
include calculations within the delta- hole model | 245 |, the role of NN correlations and 
the energy dependence of nuclear response functions ||246|| , possible effects of "Brown- Rho 
scaling" on nuclear structure functions ||243[| , and implications of low-energy pion-nucleus 



scattering data for nuclear deep-inelastic scattering and Drell-Yan production |[244 



Some further investigations use a relativistic many-body approach to treat mesonic and 
binding corrections to reproduce nuclear effects in the EMC and Drell-Yan measurements 
247 ,1248]. From the point of view of nuclear many-body theory, the best nuclear wave 



functions have been employed in Refs. [[249| , [250|1 , treating both short and long range corre- 
lations in nuclear matter and helium at the most advanced level. Two-nucleon correlations 
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turn out to be important in nuclear DIS. It should be mentioned, however, that the re- 
sults of ||247|J249| , |250| | have met with some debate concerning the proper choice of the "flux 
factor" . Questions of rigorous baryon number conservation M have also been raised. 



It has been suggested to investigate intrinsic properties of bound nucleons in semi-inclusive 
deep- inelastic scattering from nuclei [j235|,[251|,|252|,|253||. Measuring the scattered lepton 



in coincidence with the residual nuclear system should provide detailed information on 
changes in bound nucleon structure functions. Possible experiments are discussed at HER- 
MES |25l. 



7 Deep-inelastic scattering from polarized nuclei 



Understanding the spin structure of the proton and the neutron is a central issue in 
QCD. Both the polarized neutron and proton structure functions, g\ and gf, are needed 



in the investigation of flavor singlet quark spin distributions (see e.g. Ref . | iifl ) , and in the 
experimental test of the fundamental Bjorken sum rule ( |2.60| ). 



Since free neutron targets are not available one must resort to polarized nuclei, such as 
the deuteron and 3 He, where the neutron spin plays a well defined role in building up the 
total polarization of the nuclear target. Polarized deep- inelastic scattering from deuterium 
and 3 He I^,^,l4l[] has been studied with high precision. In order to deduce 



accurate information about the individual nucleon spin structure functions from these 
data, it is essential to correct for genuine nuclear effects. In addition, the presence of 
the tensor interaction between nucleons in nuclei creates specific spin effects which are of 
interest in their own right. 



7. 1 Effective polarization 



As discussed in Section J|, nuclear structure functions at Bjorken-x > 0.2 are dominated 
by the incoherent scattering from bound nucleons. For polarized nuclei, the non-trivial 
spin-orbit structure of the wave functions causes new effects. Bound nucleons can carry 
orbital angular momentum, so their polarization vectors need not be aligned with the 
total polarization of the target nucleus: depolarization effects occur. 

In order to describe such nuclear depolarization phenomena it is useful to introduce effec- 
tive polarizations for nucleons bound in the nucleus. Let \A\) represent a nuclear state 
polarized in the ^-direction. Then the effective polarization of protons or neutrons in that 
nucleus is 
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V^={A]\ £ a z (i)\A1), (7.1) 

i=protons 

^n=^Tl E *.(OI^T>- (7-2) 

i=7ieMiro?i.s 

When the nuclear depolarization effects are described entirely in terms of these effective 

polarizations of bound nucleons, the nucleon spin structure functions (say, g± and g£ ) 
have the following simple additive form: 

gt 2 (x,Q 2 ) =V£gUx,Q 2 )+V^gl 2 (x,Q 2 ). (7.3) 



Nuclear depolarization effects are important over the whole kinematic range of recent 
measurements. Within the impulse approximation these effects exceed by far the influence 
of nuclear binding and Fermi motion at 0.2 < x < 0.7 (see Section TA ). 



7.2 Depolarization in deuterium and 3 He 



In case of the deuteron the effective proton and neutron polarizations are simply deter- 
mined by the D-state admixture in the deuteron wave function, induced by the tensor 
interaction between proton and neutron in the spin triplet state. One finds: 

Vt =Vt = l- \Pv (7.4) 



where Pd is the D-state probability. The numerical values of Pp n range between 0.91 
and 0.94 using deuteron wave functions calculated with the Paris [152] or Bonn [195] 
nucleon-nucleon potential, respectively. 



Apart from the interest in neutron spin structure functions the deuteron with its triplet 
spin structure is of interest all by itself. Its spin-1 property leads to additional structure 
functions as given in Eq.( |3.3| ). In particular, the new spin structure functions b\ and b 2 are 
accessible in deep-inelastic scattering from polarized deuterons and can be investigated 
in forthcoming HERMES measurements [p5| . 

Polarized 3 He can be viewed, to a first approximation, as a polarized neutron target, with 
the proton-neutron subsystem in a spin singlet configuration and the surplus neutron 
carrying the spin of the three-body system. Corrections to this picture come from the 
admixture of S'- and D-wave components to the 3 He wave function. The consequence is 
that the effective neutron polarization is reduced from unity, and the effective proton 
polarization does not vanish: 
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Vl c = 1 - ^P S ' - \Pv = 0.86 ± 0.02, (7.5) 
Vf e = ~ (P D - Ps>) = -0.056 ± 0.008. (7.6) 

These results are obtained from three-body calculations using realistic nucleon-nucleon 
interactions ||255|| , omitting however effects from meson exchange currents (see Section 

7]y]) 



7.3 Nuclear coherence effects in polarized deep-inelastic scattering 



Coherence phenomena such as shadowing at small Bjorken-x are dominated by the inter- 
action of diffractively excited hadronic states with several nucleons over large longitudinal 
distances in the target nucleus. The characteristic space-time properties of DIS are inde- 
pendent of the target or beam polarization. Therefore, nuclear coherence effects are also 
expected in polarized scattering. We explore such effects in the following for deuterium 
and 3 He. 



7.3.1 Polarized single and double scattering in the deuteron 

Consider the deuteron spin structure functions gf and b\ at small values of the Bjorken 
variable, x < 0.1. Following the discussion in Section [3.2| , these structure functions can 
be expressed in terms of virtual photon-deuteron helicity amplitudes. At large Q 2 in 
the Bjorken limit which we keep throughout this section, only the helicity conserving 
amplitudes enter. As usual, we choose a right-handed, transversely polarized (virtual) 
photon (index "+") for reference. We denote the helicity conserving 7*d amplitudes by 
A\^, where H = 0, +, — refers to the helicity state of the polarized deuteron, and we 
choose the direction defined by the photon momentum q as quantization axis. The spin 



structure functions of interest are then expressed as ( |3.7| , |3.9|) 



,f = ^ImUt d -^; + d ), (7-7) 



4ne 

bx = (2^ d - AX? - AX d ) • (7.8) 

Let us now decompose -4.+^ into incoherent, single scattering terms and a coherent double 
scattering contribution. We will use the non-relativistic deuteron wave function, 



Mr) 



u{r) v(r) 1 ^ 

1 ^12M 



Xh, (7.9) 
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1 7*1, and xh denotes the S — 1 spin wave function of the deuteron. The tensor 



where r 

operator Si 2 (r) = 3(<x p • r)(cr n ■ r)/r 2 — cr p • <x n , and u(r) 
radial wave functions normalized as / °°dr [u 2 (r) + v 2 (r)] 



, v(r) are the S- and D-state 
1. The D-state probability is 
P D = / °°drf 2 (r). We have P D ~ 5.8% for the Paris potential [fL5g] and P D ~ 4.3% for 



the Bonn potential ||195 



In the polarized deuteron, the proton or neutron can have their spins either parallel or 
antiparallel with respect to the z-axis defined by q/\q\ - Let the corresponding projection 
operators be P?' n and Pf ,n , respectively. The amplitude for single scattering of the virtual 
photon from a proton in the polarized deuteron is: 



A& = J d 3 r^(r) (pfAtf + PfA£?) Mr), (7.10) 

with the helicity conserving 7*-proton amplitudes A\.^ and A\f . The analogous ampli- 
tudes for single scattering from the neutron are obtained by the replacement [p <-> n]. We 
then have A\^ = ^+if a t the single scattering level. Combining with Eqs.( |2.30"| , |7.7|) 
one finds 

9i=(l-\p»){9l + gl)=2V A g^ (7.11) 

where V d is the effective nucleon polarization ( [7.4|) in the deuteron and gf = (g± +gf)/2. 
Of course, b\ = at the single scattering level since nucleons as spin 1/2 objects do not 
have a structure function b\. 

Next we concentrate on the coherent double scattering amplitude, 

&A\ H = A\ H — (A+H + A+ H ) (7-12) 



which simultaneously involves both the proton and the neutron. At x < 0.1 this ampli- 
tude is dominated, as in the unpolarized case, by diffractive production and rescattering 
of intermediate hadronic states, but now with the polarized target nucleons excited by po- 
larized virtual photons. We introduce the diffractive production amplitudes T™(k) and 
T™(k) which describe the diffractive production process 7*N — ► XN for right handed 
photons on polarized nucleons, with momentum transfer k. Following steps similar to 
those described in Section |5.1.1| , one finds |70| , |256|| : 



A v ' — oo 

4>Ur) (P?T$(k) + PfT^f(k)) (P?T^(k) + P?T$(k)) Mr), (7.13) 
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with the longitudinal propagation length A = 2v (M x + Q 2 ) -1 of the diffractively produced 
intermediate system. We recall from Section [4.5| that a hadronic fluctuation of mass Mx 
contributes to coherent double scattering only if its propagation length A exceeds the 
deuteron diameter, (r 2 ) 1 ^ 2 ~ 4fm. 

In the following we approximate the dependence of the diffractive production amplitudes 
on the momentum transfer t = k 2 « —k\ by: 

T XN (fc)^e- Bfc i/ 2 T XN , (7.14) 



with the forward amplitude T XN = T XN (fc = 0). Various data on diffractive leptoproduc- 
tion at Q 2 < 3 GeV 2 suggest an average slope B ~ (5 . . . 10) GeV -2 (for references see 
e.g. |54],[ljJ). We then define the integrated (longitudinal) form factor 

• F|(A " ) = Jj^Y 2SH{k± > A " } e ~ Bkl > (7 ' 15) 



where 

S H (k)= Jd 3 r\Mr)\ 2 e lkr (7.16) 



is the conventional helicity-dependent deuteron form factor. Next we introduce helicity 
dependent diffractive production cross sections for transversely polarized virtual photons 
by: 



1 
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(7.17) 



t=o 



9 7* N 

and a corresponding expression for da, , with the center-of-mass energy W of the 7*N 
system. The resulting coherent double scattering correction to the spin structure function 
gf becomes [ffO| : 



5g*(x,Q> 



2Q 2 

e 2 x 



w 1 



dM x 



2 TIN' 



d 2 a 



dM^dt dM^dt 



B l \ —In 



t=o 



(7.18) 



Similarly one obtains for b\ from Eqs.( [7.8| , |7.13| ) [ [70| , |194|j256 

w 2 



2Q 2 

e 2 x 



dm 



d 2 a^ 



Am-L 



dM 2 x dt 



t=o 



T*{\- x ) + ^(A- 1 ) - 2F*{\- 1 



(7.19) 



101 



7.3.2 Shadowing in gf 



The difference of polarized diffractive virtual photoproduction cross sections which enters 
in Eq. (|7.18|) has so far not yet been measured. Nevertheless it is possible to estimate 



the shadowing correction 8g\ to an accuracy which is sufficient for a reliable extraction 
of the neutron structure function g\ from current experimental data. With inclusion of 
shadowing and the effective nucleon polarization in the deuteron one finds: 

A ^-^ A-A- (7.20) 



(V 



The measured spin structure functions have the property \gf\ < \g\\, at least for x > 



0.01 p^| . This implies that, at the present level of data accuracy, shadowing effects and 
uncertainties in the deuteron D-state probability do not play a major role in the extraction 

To estimate the amount of shadowing in gf and its influence on the extraction of gf one can 
study the double scattering contribution ( 7.18|) in the framework of a simple model. In the 



laboratory frame at small x < 0.1 the exchanged virtual photon first converts to a hadronic 
state X which then interacts with the target (see Section |4.5|) , dominant contributions 
coming from hadronic states with invariant mass M\ ~ Q 2 . Consider therefore a single 
effective hadronic state with a coherence length A ~ l/(2Mx). Comparing shadowing for 
unpolarized and polarized structure functions gives |f70| : 



2gf 



n 



5F[ 

If 



n 
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5Fj 



with 1Z 



> J r ^(2Mx) 
V B (2Mx) 



(7.21) 



At small x and B = 7 GeV one finds 1Z 



2.2 for both the Paris and Bonn nucleon- 



nucleon potentials |152|J195| . Although shadowing for gf turns out to be approximately 



twice as large as for the unpolarized structure function F 2 , it still leads to negligible effects 
on the extraction of gf, at least at the present level of experimental accuracy. Using the 

2 



experimental data on shadowing for F| 
( [7.201 ) amounts at x ~ 0.01 to less than 
analysis ||257 . 



91 



one finds that the shadowing correction in 
i of the experimental error on gf for the SMC 



7.3.3 The tensor structure function b\ at small x 

The shadowing correction for the unpolarized structure function 8Ff and the deuteron 
tensor structure function b\ are directly related. In order to see this, note again that 

8 Note that shadowing corrections in unpolarized structure functions, 5Fi 2 = F^ 2 ~ {F\ 2 ~ 
i^ 2 )/2, are defined per nucleon, as Ff 2 themselves. 
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Fig. 7.1. Double scattering contribution to the tensor structure function b\ from |70| . The 
dotted and dashed curves correspond to the Bonn |195| 1 and Paris potential [152 1 , respectively. 



the propagation lengths ( 4.19|) of diffractively produced hadrons exceed the deuteron size 
A > (r 2 )d ~ 4fm at small x. The deuteron form factors become approximately constant, 
i.e. .F#(A _1 < 1/4 fm) w JF^(O), and a comparison with the double scattering correction 



for the unpolarized structure function (|5.15 ) gives [70 



b 1 = K bl 5F 1 , with K 



In 



F*(0) - g(0) 



(7.22) 



With B = 7 GeV 2 we find from the Paris nucleon-nucleon potential ||152|| IZ^ = —0.66, 
while the Bonn one-boson-exchange potential ||195|| leads to TZb 1 = —0.58. Using data for 



F2/F2 one can estimate b\ at small x and finds that it reaches about 2% of the 



unpolarized structure function Ff at x 0.1 J7DLESSJ. 



In Fig. |7.1| we present b\ as obtained from Eq. ([7.22 ). The result shown here corresponds 

Estimates of b\ at large Q 2 3> lGeV 2 and small 1 C 0.1 



to the kinematics of E665 
can be found in 
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194 



It should be mentioned that the magnitude of b\ at x < 0.1 exceeds 
estimates from previous model calculations which are applicable at moderate and large 
x, by several orders of magnitude (see e.g. [[259| , |260|| ). Unfortunately, the effect of b\ in 



the observable asymmetry, which is proportional to b\/F^, is only of the order of 10 2 , as 
already mentioned. 
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7.3.4 The 3 He case 



Although 3 He would appear to be an ideal neutron target because of its large effective 
neutron polarization ( |7.5| ), the extraction of neutron spin structure functions nevertheless 
requires dealing with significant nuclear effects. We concentrate again on the structure 
function g\ at small values of a;. As for deuterium two types of corrections are relevant: 
higher angular momentum components of the 3 He wave function lead to effective proton 
and neutron polarizations ( |7.5j , [r.6| ). Furthermore, the coherent interaction of the virtual 
photon with several nucleons causes shadowing. Accounting for both effects the neutron 
structure function can be deduced from gf e as follows: 

-pHer He/ n -pHe 
n ' n °Jh I9l He _ P (n nn^ 

9l (pHe)2 9\ p H e 9l- {' 



Since \V^ e \ <C \V^ e \ the proton contribution in (|7.23|) is indeed suppressed 



Uncertainties in the effective nucleon polarizations which may influence the extraction of 
g^, result also from non-nucleonic degrees of freedom, e.g. mesons and A-isobars, present 
in 3 He. In Ref. [261] this is demonstrated for the non-singlet nucleon and A = 3 structure 
functions: 



Agt=\x,Q*) 
Ag?= 3 (x,Q 2 ) 



9i(x,Q 2 ) 



9i(x,Q 2 ), 
■9?(x,Q 2 ) 



(7.24) 
(7.25) 



where g^ is the triton spin structure function. From the effective nucleon polarizations 
( [7.5| . |7.6| ) and isospin symmetry of the three-body nuclear wave function one obtains: 



Ag( 



A=3 



3 b 



(7.26) 



Applying the Bjorken sum rule ( |2.60| ) to the A 
ratio leads to: 

Jo 1 dx [gf c ( Xl Q 2 ) - gf(x, Q 2 )] Ga (^) 



1 and A = 3 systems and taking their 



(7.27) 



The axial vector coupling constants G^( 3 H) and G^(n) are measured in the /9-decay of 
tritium and the neutron, respectively, with GU(n) = 9a = 1-26 [[262 1. If one considers 
incoherent scattering from individual nucleons and accounts only for the effective nucleon 
polarizations one finds from Eqs.( |7.26| , |7.27|) : 



GU( 3 H) 
9a 



4 2 

i - - 3 p, - - 3 p„ 



0.922 ±0.006. 



(7.28) 
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This result is however at variance with the empirical ratio Ga{ 3 H)/ gA ~ 0.963 ± 0.003 
One concludes that simply using effective nucleon polarizations from realis- 



tic three-nucleon calculations would lead to a violation of the Bjorken sum rule by ap- 
proximately 4% [261 , 264] . On the other hand, it is known from nuclear /5-decay and 



Gamov- Teller transitions that axial coupling constants in nuclei are renormalized by me- 
son exchange currents and A-isobars ||265| . The possible influence of these non-nucleonic 



degrees of freedom on the extraction of gf has to be carefully investigated ||261 |. 



At small values x < 0.1 of the Bjorken variable, coherent multiple scattering from several 
nucleons in the target leads to the shadowing correction 5g±. Assuming again, as in Section 



7.3. 2| , that the photon-nucleus scattering at small x can be represented by the interaction 
of one effective hadronic fluctuation with invariant mass M\ ~ Q 2 , one finds [|261||: 



gf 



F 2 N 



- 1 



(7.29) 



i.e. shadowing in g^ e is about twice as large as for the unpolarized structure function 
F^ e . This is also true for shadowing in the non-singlet structure function Ag^ =3 |[261|| . A 
similar result has been found for the polarized deuterium case, see e.g. ( [7.21[ ). 



In Ref. ||261|| the shadowing correction ( |7.29| ) has been combined with the nuclear Bjorken 



sum rule ( |7.27|) . In the nuclear case, shadowing reduces the small-x contribution to the 
Bjorken integral. Consequently the non-singlet nuclear structure function must be en- 
hanced, i.e. Ag^ =3 > Ag± =1 , somewhere in the region x > 0.1. This is suggested to 
occur around x ~ 0.1, where the projectile may still interact with two nucleons inside 
the target. As a consequence a significant anti-shadowing is obtained in this kinematic 
region. The dynamical origin for such an enhancement is supposed to be independent of 
the flavor channel considered, so that this is also expected to occur for gf e - 

To summarize this discussion we emphasize that a precise extraction of the neutron spin 
structure function gf from 3 He data at x < 0.2 requires a careful analysis of nuclear 
effects due to non-nucleonic (e.g. meson and A-isobar) degrees of freedom in the 3 He 
wave function and a detailed understanding of shadowing and anti-shadowing effects. The 
use of deuterium as a target may have advantages since non-nucleonic admixtures are 
supposed to be smaller due to weaker nuclear binding. Furthermore, especially at small 
x the influence of shadowing on the extraction of gf is less pronounced for deuterium as 
compared to 3 He. 



7.4 Polarized deep-inelastic scattering from nuclei at x > 0.2 



At moderate and large values of x, the distances probed by DIS from nuclei are smaller 
than 2 fm as outlined in Section [|. Incoherent interactions of the virtual photon with 
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hadronic constituents of the nucleus dominate, and the usual starting point is a de- 
scription based on the impulse approximation including nucleon degrees of freedom only 
213y255| . f266| . [267|| . This approach is not without question, but not much progress has so 



far been made beyond that level. 

In the following we first discuss general features of the nuclear spin structure functions 
and g? = 9i + 92 m ^he framework of the impulse approximation. The nuclear target 
is treated non-relativistically. This allows, as in the unpolarized case, to factorize nucleon 
and nuclear degrees of freedom, introducing structure functions of nucleons as bound 
quasi-particles. The nuclear structure functions are then given as two-dimensional convo- 
lutions of the bound nucleon structure functions and the nucleon light-cone momentum 
distributions |]213]: 



9?(x)=£[d P 2 I^D\{y^)g\( X -^ 
75 J J y \y 



i=i 

A 



[dp 1 

i=i j t v 



DUyy)9T\zy)+ D T2(y,p 2 )92^p 2 



(7.30) 
(7.31) 



Here we have suppressed the Q 2 -dependence for simplicity. The expressions in Eqs.( |7.30| , 
7.31|) resemble the result for the unpolarized nuclear structure function in Eq. (|6.6|) : the 
exchanged virtual photon scatters from quarks which carry a fraction x/y of the light- 
cone momentum of their parent nucleons, which in turn have a fraction y = p + /M of the 
nuclear light-cone momentum and a squared four- momentum p 2 . 



The nucleon distribution functions in Eqs.( [7.30| , |7.31|) are given by: 



Di(y, P 2 ) 
D T (y,P 2 ) 
D T 2(y,p 2 ) 



dy 

(2tt)* 
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tr 



tr 



S±( P ') E ± ]6[y-^)6( P 2 -p' 2 



S ± (p')%}5( y -^\8(p 2 -p' 2 



with the polarized nucleon spectral function: 

S aa , (p) = 2tt ]T 5 (po - M - e n + T R ) i, ff (p)< y (p). 



(7.32) 
(7.33) 
(7.34) 

(7.35) 



Here the summation is performed over the complete set of states with A — 1 nucleons. 
The functions i[>n,a(p) — ((A — l)n? — p|"0cr(O) | A) are the probability amplitudes to find a 
nucleon with polarization a in the nuclear ground state and the remaining A — 1 nucleons 
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in a state n with total momentum —p. The separation energy e n and the recoil energy 
of the residual nuclear system T R enter in Eq.( [7.35| ) as for the unpolarized case d6.5|). 
For gi the target nucleus is chosen to be polarized parallel to the photon momentum as 
indicated in Eq. (|7.32| ) by the subscript ||. For the transverse structure function g T the 
target polarization is taken perpendicular to the momentum transfer and denoted by _L. 



The nucleon spin operators which multiply the spectral functions in Eqs. ( |7.32| , |7.33| , [7T3^ ) 
refer to the active nucleon. They read: 



er • p 

M ' 

P±-Sj_ 

M 



<T ■ p 



a ■ p 

M 



+ cr z 1- 




(7.36) 



where j denotes spatial indices. Furthermore, p ± is the transverse component of the nu- 
cleon three-momentum, p = {p±,p z ), and S± determines the transverse spin quantization 
axis relative to the photon momentum which is taken along the z-direction. 



From Eqs. ( |7.30| , |7T31|) we observe that g^ is expressed entirely in terms of the corresponding 



nucleon structure function g^ . This is different for gtf which receives contributions from 
g? as well as from g%. 



If the bound nucleon structure functions in Eqs.( |7.30"l , |7.31|) are replaced by free ones, one 
ends up with the conventional one- dimensional convolution ansatz for nuclear structure 
functions ||213| , |255|| . Relativistic contributions which lead beyond the convolution formula 
( |7.30| , [7.3lD have been investigated in Refs. |[268| , [269|| , and corrections have been estimated 



within a quark-diquark model for the bound nucleon. Deviations from non-relativistic 
convolution were generally found to be small, except at very large x > 0.9. 

Given that the nuclear spectral functions, polarized as well as unpolarized, receive their 
major contributions from small nucleon momenta, systematic expansions of g^ and g^ 
for x < 0.7 can be performed around y = 1 and the mass-shell point p 2 = M 2 , keeping 
terms of order e/M and p 2 /M 2 . Results and applications to spin structure functions of the 
deuteron and 3 He are discussed in Refs. [pl3|J255|J267| , |270|| . As a general rule, the structure 
functions at x < 0.7 are well described using simply the effective nucleon polarizations 



8 Further developments and perspectives 



We close this review with a short summary of the key physics points together with an 
outlook on several selected topics for which investigations are still actively under way. We 
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comment on exclusive vector meson production from nuclei, questions of shadowing at 
large Q 2 and the issue of high parton densities in nuclear systems. 



8. 1 Coherence effects in DIS and in the exclusive electroproduction of vector mesons 



Nuclear shadowing in inclusive deep-inelastic lepton scattering is a prime source of in- 
formation on coherently propagating hadronic or quark-gluon fluctuations of the virtual 
photon in a nuclear medium. By selecting different kinematic cuts in Q 2 and energy trans- 
fer v 1 one can focus on different components of the photon's Fock space wave function. 
An even more stringent selection of such components can be achieved in exclusive photo- 
and electroproduction processes, and in particular in high-energy diffractive vector meson 
production. Data on vector meson production from nuclei have become available in recent 
years at FNAL (E665) f27I], CERN (NMC) PT2] , |2T31 , and DESY (HERMES) JZ7§, and 



further experiments are under discussion at TJNAF [275] and DESY (HERA, HERMES) 
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Depending on energy and momentum transfer, the mechanism of vector meson formation 
can be quite different. In the range v > 3 GeV and Q 2 < 1 GeV 2 the production process 
is well described using the vector meson dominance picture (see e.g. P5|1 ): in the lab 
frame the photon converts into a vector meson prior to scattering from the target. On the 
other hand, at large Q 2 3> 1 GeV 2 perturbative QCD calculations show that the photon- 
nucleon interaction produces an initially small-sized, color singlet quark- ant iquark wave 
packet ||166| , |276|| . At high photon energies the finally observed vector meson is then formed 
at a much later stage. 

The transition from small to large Q 2 interpolates between non-perturbative hadron for- 
mation and perturbative quark- ant iquark-gluon dynamics, a question of central impor- 
tance in QCD. Nuclear targets are particular helpful at this point because they serve 
as analyzers for the coherent interaction of the produced gg-gluon system with several 
nucleons |[277| . The distance between two nucleons provides the "femtometer stick" which 



can be used to measure the relevant coherence lengths (for reviews and references see 



278lg79iP80| ) 



The characteristic scales for this discussion have been encountered several times in previ- 
ous sections. First, there is the typical longitudinal distance (propagation length) 



It represents the distance over which a hadronic fluctuation of invariant mass m propagates 
in the lab frame when induced by a photon of energy v and virtuality Q 2 . At large Q 2 the 
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Fig. 8.1. The rates of produced p 0, s plotted against t' = t — to for hydrogen, deuterium, 3 He, and 
14 N. The solid lines show fits to the data, the dashed lines are inferred incoherent contributions 



(for details see Ref.j27| 



initially produced wave packet is characterized by its transverse size b. For longitudinally 
polarized photons, 



const 



Q 



(8.2) 



In perturbative QCD the minimal Fock space component has const ~ 4-5 at Q 2 > 5 GeV 2 
Thus for Q 2 = 5 GeV 2 and const = 4, the transverse size of the initial wave 



packet is b ~ 0.4 fm, a small fraction of the diameter of a fully developed p meson. 

Recent measurements at HERMES [[274|| have observed effects related to the coherence 
length A in p° electroproduction on hydrogen, deuterium, 3 He, and 14 N. The range of 
energy and momentum transfers covered by the experiment is 9 GeV < v < 20 GeV and 
0.4 GeV 2 < Q 2 < 5 GeV 2 . This implies coherence lengths in the range 0.6 fm < A < 
8 fm covering scales from individual nucleons up to and beyond nuclear dimensions. (The 
interesting upper section of the available Q 2 interval, Q 2 > 4 GeV 2 , has been accessible 
only for small energies with v with A < 1 fm in these measurements.) 

Given the four-momenta q and k of the incoming virtual photon and the produced p 
meson, the t-channel momentum transfer to the nucleon is t — (q — k) 2 . In Fig.^TT] the 
rate of produced p°'s is plotted against t! = t — t , the squared momentum transfer above 
threshold (\to\ — A~ 2 ). At \t'\ < 0.1 GeV 2 coherent production dominates, leaving the 
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Fig. 8.2. The nuclear transparency Ta as a function of the propag ation length A for 14 N. HER- 
MES data p7l are indicated by full dots. The open symbols represent data from previous 
experiments [ 284, 271 1. The dashed curve shows a Glauber calculation from Ref . |282] 



nucleus as a whole in the ground state. Such coherent processes fall off rapidly with the 
nuclear form factor, so that at > 0.1 GeV 2 mostly incoherent p production from 
individual nucleons remains. 

Consider now the incoherent production of vector mesons from nuclei. In the absence of 
coherent rescattering processes the nuclear p production cross section oa would simply be 
A times the production cross section <tn on a free nucleon. Nuclear effects are conveniently 
discussed in terms of the transparency ratio, Ta = oa/(^4ct n ). The measured ratio for 
14 N is plotted as a function of the longitudinal propagation length A in Fig. |S.2| . The 
deviation of Ta from unity for A < 1 fm simply reflects the "trivial" final state rescattering 
of the p meson after being produced on one of the nucleons. More interesting effects 
are visible when A exceeds the average nucleon-nucleon distance of about 2 fm. Now 
the hadronic fluctuations of the photon can scatter coherently on several nucleons also 
prior to the production of the final state vector meson, and the transparency ratio Ta 



systematically decreases until it exceeds the nuclear diameter. The dashed curve in Fig. 3.2 



shows a theoretical prediction calculated within the vector meson dominance model ||282 . 
Its agreement with data suggests that the production process is dominated, given the 
relatively low Q 2 involved, by hadronic fluctuations which interact about as strongly as 
the produced p meson. Further systematic investigations of such coherence length effects, 
especially its detailed dependence on the momentum transfer t, are discussed at TJNAF 

It is interesting to push these observations to more extreme regions of very large Q 2 and 
v. Once the energy transfer exceeds several tens of GeV, a further scale enters. At large 
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Q 2 the qq fluctuation of the photon starts out initially as a small-sized wave packet. The 
time it takes for this wave packet to develop into the final vector meson is called the 
formation time tj ||278| , |279| , |280|| . To be specific, let the observed vector meson again be a 



p . The small initial wave packet is generally not an eigenstate of the strong interaction 
Hamiltonian. Now consider expanding the wave packet in hadronic eigenstates. Clearly, 
for a wave packet with a size small compared to typical hadronic dimensions, several such 
eigenstates are necessary to represent the wave function of the packet. Let one of those 
hadronic eigenstates in the expansion be the p meson itself (mass m p ), let another one be 
a neighboring state with larger mass (say, m p i = m p + 5m). The characteristic propagation 
length of the p component of the gg-gluon wave packet is A ~ 2v/{m 2 + Q 2 ), that of the 
neighboring state is A' ~ 2v j {m 2 + Q 2 + 2m p 8m) . The phase difference between both 
states is determined by their wave numbers 1/A and 1/A'. The time it takes to filter out 
all but the p meson component when passing through the nucleus is then 

* m p 5m ^ 

Thus when v is sufficiently large so that Tf reaches nuclear dimensions and Q 2 is large, 
a small-sized qq wave packet induced by the photon has a chance to travel over large 
distances inside nuclei and interact weakly, its cross section being proportional to b 2 ~ 

i/Q 2 . 

This phenomenon is commonly named color coherence or color (singlet) transparency 
277| . |278| . p79tp80|| . It has been addressed in exclusive vector meson production experiments 



at FNAL (E665) [^7T| and CERN (NMC) pT5j . The interpretation of data in terms of 
color transparency is still under debate (see e.g. ||285| , [286fl ). Possible future developments 



with nuclei at HERA ||254|| could offer an enormous extension of the accessible kinematic 
range. 

8.2 Nuclear shadowing in DIS at large Q 2 



Experimental information on nuclear shadowing in inclusive DIS is available, up to now, 
only from fixed target experiments, with the kinematic range restricted to Q 2 < 1 
GeV 2 at x <C 0.1. Although there are currently no data on nuclear shadowing at large 
Q 2 ^> 1 GeV 2 , it is nevertheless instructive to investigate what one should expect in 
comparison with the previous results at smaller Q 2 . Such a study is possible due to re- 
cent data on diffractive production from the HERA collider and has been performed in 



Refs. 147.287,288 



One finds that contributions from vector mesons are negligible at large momentum trans- 
fers since diffractive vector meson production is strongly suppressed at Q 2 > 10 GeV 2 . 



Furthermore, the ZEUS data [52] on the ratio of diffractive to total photon-nucleon cross 
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sections (Fig. |2.8| ) suggest that diffractively produced states with large mass are relevant at 
large Q 2 . Therefore at large Q 2 shadowing probes the coherent interaction of quark-gluon 
configurations with large invariant mass. This is complementary to fixed target experi- 
ments at FNAL (E665) and CERN (NMC) where the coherent interaction of low mass 



vector mesons played a dominant role (see Sections |5.1.4| and 



In Section |5.2| we have argued that to leading order in 1/Q 2 shadowing is dominated by 



the interaction of large-size hadronic fluctuations of the exchanged photon. This suggests 
a weak energy dependence of nuclear shadowing. However, at very small x together with 
very large Q 2 , the steadily growing number of partons in the photon- nucleon system makes 
quark-gluon configurations of the photon interact like ordinary hadrons, even if they 
have small transverse size proportional to 1/Q 2 (see e.g. Eq. Q5.34| )). As a consequence 
one expects a more rapid energy dependence of shadowing as compared to the case of 
small Q 2 . On similar grounds the diffractive leptoproduction cross section at x <C 0.1 and 
Q 2 3> 1 GeV 2 should rise more strongly than suggested by Regge phenomenology. Possible 
indications for this behavior have been found at HERA |50| , |52fl , signaling the onset of a 
new kinematic regime with a complex interplay between soft (large size) and hard (small 
size) partonic components of the interacting photon. A systematic investigation of strong 
interaction dynamics in this kinematic region is a major challenge. Here electron-nucleus 
collider experiments could give important new insights ||254| . 



8.3 Physics of high parton densities 



An elementary QCD treatment of radiative corrections in nucleon and nuclear structure 
functions, or equivalently, in quark and gluon distributions is provided by the DGLAP 
evolution equations (see Section |2.4| ). In leading logarithmic approximation one sums 
over a s lnQ 2 /A 2 terms, each of which represents a quark radiating a gluon or a gluon 
splitting into a qq or gluon pair. A contribution (a s lnQ 2 /A 2 ) n is associated with the 
radiation from n partons in a physical gauge. Due to radiation partons loose momentum. 
Therefore DGLAP evolution leads to a rise of quark and gluon distribution functions at 
small x: strongly increasing numbers of partons carry smaller and smaller fractions of the 
total momentum. For example, consider gluons which dominate the dynamics of parton 
distributions at small x. Suppose that we are given an initial gluon distribution, prior 
to radiative QCD corrections, at a momentum scale Qq. This initial distribution can be 
non-singular (finite) at x — > 0. Now turn on QCD radiation and DGLAP evolution. The 
resulting asymptotic behavior of the gluon distribution function at x <C 1 and large Q 2 is 
(see e.g. 




xg(x, Q 2 ) ~ exp { In HV*) ) ■ (8-4) 
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One observes a strong rise of xg(x,Q 2 ) at i <C 1. This increase of the gluon density 
with decreasing x can however not continue indefinitely. At very small x the density of 
gluons becomes so large that they interact with each other reducing their density through 
annihilation. Thus one enters a new regime of very high parton densities where standard 
methods of perturbation theory are inappropriate despite small values for the coupling 
a s , and resummation techniques to infinite order have to be applied (for references see 
e.g. [ p4]p5|p6| , |i37| , |138| , |139|| ) . 



For further illustration let us view the scattering process in a frame where the target 
momentum is very large, P = \P\ — ► oo. A measurement at specific values of Q 2 and 
x probes partons over a longitudinal distance Az ~ l/(Px) and with a transverse size 
Ab ~ 1/ \[Q 1 . Thus at small values of x and moderate Q 2 parton wave functions overlap, 
leading to high parton densities, so that the rate of parton-parton annihilation processes 
increases. This rate involves the probability of finding at least two partons per unit area 
in the nucleon |[289|J29C1|| . In this respect parton recombination can be regarded as a "non- 



linear" correction to radiation processes described by standard DGLAP evolution, where 
parton distribution functions enter linearly (see e.g. Eq. ( J2.39) )). 



The existing HERA data on free nucleon structure functions at Q 2 > 1 GeV 2 and x > 
10~ 4 do not show a clear sign for the need of additional parton fusion corrections in the 
evolution equations (see e.g. ||). However, phenomena related to high parton densities 
at small x should be magnified in nuclei since nuclear parton densities are enhanced. For 
example, if shadowing is ignored and the nuclear gluon distribution is assumed to be the 
sum of the gluon distributions of individual nucleons, gx(x,Q 2 ) = Ag^(x,Q 2 ), then the 



ratio of gluon densities in a nucleus and a nucleon per (transverse) area is |291|: 



9AM 2 ) I g N (x,Q 2 ) _ Ar| _ Q 5 ^ (g g) 



ttR 2 a I nr^ R\ 



where « 0.8 fm and Ra ~ tqA 1 ^ fa 1.1 fm A 1 / 3 have been used for nucleon and nuclear 
radii, respectively. Effects beyond DGLAP evolution should therefore be amplified in 
nuclei and set in already at a larger value of x as compared to free nucleons. 



In Ref. [P90|J292|| corrections to DGLAP evolution for nuclear parton distributions have 
been calculated in the leading logarithmic approximation taking into account corrections 
due to parton-parton recombination. The nuclear gluon distribution is written in two 
parts: xgA{x,Q 2 ) = xAg^s(x,Q 2 ) + 5 ( xgA{x, Q 2 ))- The first term is associated with 
independent nucleons and evolves accordingly. The second term is the correction of interest 
here and describes the interaction of gluons from different nucleons. Its evolution due to 
gluon-gluon recombination reads ||290| , |292| : 



Q 1 



d 5(xg A (x,Q 2 )) 81 A 1 / 3 



dQ 2 



A 



16 Q 2 r 2 



o? s {Q 2 ) 



dx' 



x 



x'g^(x',Q 2 ) 



(8.6) 
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Clearly gluon-gluon recombination is enhanced in nuclei. One finds ||293| , 294 l that the ar- 
range where non-linear effects ( |8.6| ) become significant, differs for heavy nuclei and free 
nucleons by more than two orders of magnitude, assuming xg^(x) ~ x~ 



-0.2 



Investigations of unitarity constraints in hard two-body amplitudes ||166|| also suggest that 
non-linear effects at high parton densities can be significant in nuclei at Q 2 ~ 10 GeV 2 
and x ~ 10~ 4 . 



While the kinematic bounds for the applicability of "normal" DGLAP evolution are quite 
well defined, the dynamical mechanisms responsible for slowing down the rapid increase 
of parton distributions at large Q 2 and very small x are not yet clear. One is entering a 
new domain of QCD, dealing with partonic systems of high density, which presents new 
challenges. 
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